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FOREWORD 


This  report  vas  prepared  by  Rocketdyne,  a  Division 
of  North  American  Aviation,  Inc. ,  Canoga  Park, 
California,  on  Air  Force  Contract  AF04(6ll)-8515 
under  Task  No.  3850306  of  Project  No.  3850,  "Total 
Impulse  Measuring  System  for  Solid-Propellant  Rocket 
Engine  (Research)."  Contract  AF04(6ll)-8515  consists 
of  a  program  for  the  analysis  and  design  (Phase  l) , 
fabrication  and  testing  (Phase  II),  and  installation 
and  testing  (Phase  III),  of  an  accurate  (O.l^)  solid- 
propellant  total  impulse  measurement  system  for 
Edvards  Air  Force  Base.  This  report  is  submitted  to 
present  the  development  status  of  the  system,  and 
constitutes  the  Phase  I  report  in  the  over-all  program. 
It  vas  prepared  by  the  Research  Instrumentation  Group 
of  the  Rocketdyne  Research  Department. 
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ABSTRACT 


Discussion  and  results  of  the  three  main  considerations 
conprising  the  initial  (Phase  l)  effort  are  presented: 

(i)  study  and  discussion  of  the  use  site  and  existing 
measurement  system  at  lAlTJ,  with  recommendations  for 
the  installation  of  the  Rocketdyne  system;  (2)  general 
and  applied  mathematical  analyses  of  spring-mass  total 
impulse  measurement  systems;  and  (3)  mechanical  and 
electrical  design  of  the  Rocketdyne  system  to  be  con¬ 
structed.  Specific  results  of  these  separate  but 
related  program  efforts  are,  respectively:  (l)  con¬ 
currence  with  the  general  approach  taken  in  the 
existing  FIAFB  system,  but  disagreement  with  some  specific 
design,  installation,  and  operational  features  of  the 
system;  (2)  new  and  original  methematical  analyses  whose 
results  justify  the  Rocketdyne  approach  to  the 
measurement  problem;  and  (3)  new  and  original  mechanical 
and  electrical  component  designs  wtiicb  set  a  new 
standard  in  accuracy  and  performance,  producing  an 
over-all  error  breakdown  commensurate  with  program 
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The  final  results  of  effort  in  the  initial  (Phase  l)  portion  of  the  con¬ 
tract  were  twofold.  First,  production  of  a  detailed  design  for  the  total 
impulse  measuring;  system  with  the  desired  measurement  capabilities  was 
accomplished.  The  results  of  all  associated  analyses  of  this  system  indi¬ 
cate  that  its  performance,  in  anticipated  solid-propellent  studies,  will 
produce  measurements  with  total  errors  not  exceeding  O.l/,  Second,  the 
analytical  development  undertaken  in  support  of  the  design  effort  produced 
general  quantitative  results  useful  in  any  design  of  spring-mass  systems 
for  transient  thrust  and  total  impulse  measurement  systems.  The  general 
trapezoidal  excitation  of  both  one-  and  two-degree-of-f reedom  systems, 
with  viscous  damping,  w'as  ini^estigated ,  and  the  methods  used  to  obtain 
the  errors  expressions  and  tabulations  are  applicable  to  comparatively 
short,  as  well  as  long  engine  operating  durations.  Consequently,  the 
mathematical  results  obtained  w-ith  this  program  are  directly  applicable, 
for  example,  to  measurements  associated  with  pulsed-mode  operations  of 
space  engines,  /s  a  result,  analj’^ical  and  design  experience  was  obtained 
that  appears  useful  for  future  applications  in  d\*namic  force  and  impulse 
measurements . 

The  following  specific  subjects  are  discussed  and/or  analyzed: 

1.  Existing  thrust  integration  system  at  EA.FB; 

2.  Recommendations  for  installation  cf  Rocketdyne-designed 
system  at  EAFB; 

3.  General  mathematical  discussion  of  the  relationship  betw'een 
the  integrated  input  and  integrated  output  for  certain  linear 
dvTiamic  systems ; 

U.  Analytical  dynamics  of  a  spring-mass  total  impulse  measurement 
system  (one  degree-of-f reedom) :  theor\^  and  application ; 

3.  Dynamic  response  of  Rocketd\*Tie  thrust  and  total  impulse  measure¬ 
ment  system  (two  degree-of-f reedom) :  theory  and  application; 
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0. 


Theory  and  appl iration  of  paral lelogram  coraprobion-type 
flexure  Hunpension  systems; 

?•  Detailed  description  of  the  design  of  mechanix'^al  and  elecliical 
system  and  components; 

8.  Determination  of  (Duiplete  spr in^^-mass  system; 

9.  Description  of  the  design  of  addi ti  onal  instruraeiitati  on 
developed  for  operational  use  with  the  total  impulse  measurc'- 
ment  system. 
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EDWARDS  AIR  FORCE  BASE  THRUST  INTEGRATION  SYSTEM 


INTRODUCTION 

Discussions  and  use-site  inspections  with  EAFB  personnel  have  resulted 
in  agreement  regarding  the  general  measurement  system  approach  previously 
taken  by  EAFB  in  the  design  and  construction  of  equipment  and  facilities 
for  solid-propellant  evaluation  studies  (BATES).  A  total  impulse  measur¬ 
ing  system  consisting  of  a  precision  strain-gage  load  cell  supplying  an 
analog  signal  to  a  voltage-controlled  oscillator,  to  produce  a  proportional 
pulse  repetition  rate  whose  output  is  cumulatively  registered  by  an  elec¬ 
tronic  counter,  is  regarded  as  basically  sound  and  accurate  because  of  the 
continued  development  such  equipment  has  received  in  recent  years.  How¬ 
ever,  very  substantial  improvement  has  been  immediately  realized  in  this 
method  through  concentrated  design  and  analytical  efforts  directed  at  com¬ 
ponent  and  system  refinement.  As  a  result,  specific  design,  installation, 
and  operational  features  of  the  existing  EAFB  system  appear  undesirable 
and  incompatible  with  the  measurement  objectives  in  this  contract.  A 
two-part  discussion  follows  of  practical  instrumentation  engineering  con¬ 
siderations  and  their  application  to  the  EAFB  system. 

PART  I:  GENERAL  DISCUSSION 

Error  Factors 

The  design  of  a  high-precision  recording  system  for  low-level  signals 
requires  a  knowledge  of  several  factors  which  are  not  readily  apparent. 

The  three  error-producing  elements  which  are  commonly  neglected  are: 

(l)  contact  resistance,  (2)  leakage,  and  (3)  thermal  emf. 

Contact  Resistance.  Contac t  resistance  is  present  whenever  a  coup  1 i ng  is 
introduced:  connector,  terminal  strip,  program  board,  switch,  relay,  etc. 

Variations  in  contact  resistance  are  usually  caused  by  dirt  and  by  ’’dry 
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circuit”  CDiilttcts.  (l^ry  circuits  occur  when  the  current  level  is  too  low 
to  remove  oxidation  and  ’powdering”  of  contact  surfaces  )  Resistance  can 
be  controlled  by  good  maintenance  and  by  choosing  relays  and  switches 
which  provide  good  wiping  arli  »n  between  chemically  stable  contact  materi¬ 
als.  In  most  cases,  contact  lesi stance  is  nut  a  serious  problem,  but  it 
can  become  serious  if  it  is  disregarded. 

Leakage .  Leakage  occurs  wherever  a  signal  lead  wire  is  exposed,  and  it 
IS  present  at  those  devices  named  in  the  foregoing  paragraph  on  contact 
resistance  (Note:  This  is  often  called  "leakage  resistance,  '  which 
makes  it  an  inverted  term.  It  is  better  to  use  "leakage”  alone  os  a  noun 
than  t  If  chance  the  confusion  caused  by  "leakage"  and  "resistance"  being 
(ff  njiposite  inclination.)  Although  leakage  can  be  caused  by  sloppy  tech¬ 
niques  (dirt,  oil,  finger  tracks,  etc  ),  it  is  roost  often  the  result  of 
atmospheric  humidity  and  is  consequently  variable.  Its  most  annoying 
characteristic,  h<»wever,  is  that  it  is  very  difficult  to  measure. 

The  reason  for  this  measurement  difficulty  is  that  if  a  res  is  ton ce -me as - 
uring  device  is  placed  across  two  terminals  where  moisture  is  suspected, 
the  voltage  from  the  device  causes  electrolysis  of  the  m(»isture.  Gas  is 
formed  at  one  terminal  (at  least),  and  the  measurement  shows  the  resist¬ 
ance  of  the  gas  film.  The  thickness  of  the  gas  film  is  determined  by  the 
amount  of  voltage  across  the  t terminals.  This  effect  is  shown  in  Table  1 
which  gives  the  leakage  measured  between  four  j>airs  of  terminals  on  a 
mi»istened  fanning  (terminal)  strip.  The  "low-voltage  bridge"  was  as.sembled 
in  the  laboratt^ry  and  designed  to  apply  approximately  the  same  voltage 
across  the  terminals  as  occurs  from  a  transducer  signal.  It  should  be 
noted  that  a  device  like  the  Voltohmyst,  with  its  i. 5-volt  source,  does 
ni)t  introduce  sufficient  error  to  be  serious,  but  that  the  high  voltage 
source  con  give  an  extremely  erroneous  result. 

The  most  serious  problem  caused  by  leakage  occurs  when  a  bridge  network 
IS  used  as  a  transducer.  I.eakage  will  reduce  the  sensitivity,  and  if  it 
19  unsymme ti 1  cal ,  it  will  change  the  balance  point.  Very  little  leakage 
18  necessai^y;  for  example ,  using  a  350 -ohm,  3  volt  bridge,  2.5  megohms 

across  one  leg  represents  IJt  of  unbalance. 
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TABLE  1 


Type  of  Device 

Megger  (Biddle) 

Ohmineter  (Triplett) 

Voltohmyst  (RCA) 

Lov-Voltage  Bridge 
(Lab  Assembly) 


UEAKAGE  MEASUREMENTS 

Approximate 
Voltage  From  Source 

150 

45 

1.5 

0.01 


Measured 
Resistance,  ohms 

500,000 

100,000 

6,500 

3,200 


Thermal  Emf ,  Thermal  emf  is  produced  wherever  dissimilar  metals  are  in 
contact.  Fortunately,  these  dissimilar  metal  contacts  come  in  pairs 
and  are  opposite  in  direction,  so  that  if  the  two  contacts  of  the  pair 
are  at  ijie  same  temperature,  the  thermal  emf's  will  cancel.  For  example, 
a  copper  lead  wire  is  connected  to  one  brass  lug  of  a  switch.  The  other 
brass  lug  naturally  connects  to  the  copper  wire  continuation.  Basically 
(disregarding  the  other  metals  in  the  switch) ,  this  becomes  a  copper- 
brass  -brass  -copper  sequence,  and  the  emf  generated  from  copper  to  brass 
is  equal  and  opposite  to  that  generated  from  brass  to  copper,  provided 
that  the  junction  temperatures  are  equal.  This  is  the  function  of 
Becquerel's  Law  for  Intermediate  Metals,  which  says  that  a  dissimilar 
metal  may  be  inserted  in  a  line  provided  that  the  junctions  are  at  the 
same  temperature. 

However,  the  insistance  on  no  AT  (temperature  difference)  between  junc¬ 
tions  should  be  noted.  The  emf  generated  by  a  t/c  (thermocouple)  system 
is  a(AT),  in  which  OC  is  the  difference  in  potential  between  the  metals 
in  the  thermal  emf  series.  When  a  t/c  is  used  for  measurinj^  temperature, 

OC  should  be  as  high  as  is  practical,  but  in  any  interconnections  carrying 
low  emf,  OC  should  be  low  enough  so  that  a(AT)  will  be  ne^^ligible,  or 
care  should  be  taken  to  keep  AT  low.  At  one  time,  constantan  was  widely 
used  for  winding  resistors,  and  these  resistors  created  a  severe  problem. 
ydhen  connected  into  a  copper  line,  two  copper  constantan  t  c  were  formed 
and  a  slight  AT  between  the  ends  of  the  resistor  produced  considerable 
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voltage,  because  for  copper-cons tantan ,  (X  is  about  20  microvolts  per  F 
of  At.  Resistors  are  now  usually  made  of  manganin,  evanohm,  etc.,  which 
have  a  low  value  of  Jtt  with  copper.  Regardless,  design  should  consider 
that  a  series  of  elements  (switches,  resistors,  relays,  etc.)  in  a  low- 
level  signal  line  can  accumulate  erough  emf  errors  to  produce  appreciable 
signal  error.  Whenever  possible,  such  items  eis  span  and  balance  controls 
should  be  put  in  the  high-voltage  (power)  lines  to  a  bridge  transducer  and 
not  in  the  low-level  (signal)  lines. 

As  an  example,  tests  for  thermal  emf  through  the  normally  closed  contacts 
of  a  certain  relay  showed  that  it  produced  one  microvolt  per  F  of  change 
in  ambient  temperature.  Vhen  the  coil  was  activated,  the  heat  from  the 
coil  was  sufficient  to  produce  8  microvolts  through  the  closed  contacts 
after  it  had  stabilized.  This  is  a  particularly  severe  case,  but  all  the 
relays  tested  showed  some  thermal  effect. 

The  zero  balance  control  can  be  an  awkward  source  of  thermal  emf’s.  This 
control  is  often  necessary  to  bridge  circuits,  but  is  awkw'ard  because  the 
control  is  usually  a  slidewire  (VDR)  which  is  a  mass  of  dissimilar  materi¬ 
als:  brass  terminals,  resistance  wire,  and  contact  slider  (usually  a 

mixture  of  graphite  and  metal).  There  is  one  way  that  the  control  wiring 
can  be  designed  to  minimize  its  effect.  First,  the  VDR  should  be  connected 
between  the  power  input  terminals ,  not  between  the  s ignal  terminals .  Sec¬ 
ond,  it  should  use  a  limiting  resistor  of  as  high  a  resistance  as  possible. 
In  Fig.  1,  R^  is  the  transducer  bridge  resistance,  R^  is  the  VDR  (zero  con¬ 
trol)  and  R^  is  the  limiting  resistor.  kTien  wired  in  the  fashion  shown, 


the  emf’s  developed  at  the  three  VDR  terminals  will  either  be  in  the  power 
input  circuit  in  which  microvolts  are  negligible  compared  to  volts,  or 
will  be  isolated  by  the  ratio  of  R^  to  R,  .  For  this  reason,  R^  should  be 
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as  large  as  possible.  It  can  be  calculated  by  using  the  majcinium  unbalance 

condition  which  occurs  when  the  R  slider  is  at  either  end,  thus  putting 

s 

across  one  leg  of  For  a  350-ohm,  3  nrv'/volt  transducer,  500,000  ohms 

for  will  provide  more  than  ±5^  zero  control,  which  is  adequate  because 
commercial  transducers  are  (or  better)  zero  balance.  Any  transducer 

which  is  very  far  beyond  its  zero  balance  specification  has  something  basic¬ 
ally  wrong,  and  should  be  removed  from  service  for  repair.  If  is 
500,000,  then  the  effect  of  the  thermal  emf  from  the  zero  control  will  be 
reduced  by  the  factor  350/500,000  or  0.07^,  and  can  be  neglected. 

The  worst  single  source  of  thermal  emf  error  is  the  strain-gage  transducer. 
By  their  nature,  the  strain  wires  and  the  copper  lead  wires  form  t/c's,  but 
these  are  partially  compensated  for  by  Becquerel's  Law.  The  gage  wires, 
however,  have  a  finite  dimension  and  the  t/c  jiinctions  are  displaced  far 
enough  so  that  the  signal  terminals  can  be  at  different  temperatures  (i.e., 
At  is  not  zero),  and  this  may  be  complicated  by  the  presence  of  internal 
zero  and  modulus  resistors.  If  the  transducer  can  be  located  in  a  protec¬ 
tive  temperature-controlled  box  (e.g.,  pressure  transducers)  the  problem 
is  not  serious,  but  a  load  cell  measuring  force  from  ein  operating  thrust 
chamber  is  quite  vulnerable.  Radiated  heat  from  the  flame  can  reflect 
from  surrounding  elements,  convected  heat  can  be  carried  by  swirling  air, 
and  conducted  heat  can  be  transmitted  from  the  chamber  through  the  supports. 
The  most  aggravating  characteristic  of  this  problem  is  that  it  varies  inde¬ 
terminably  during  the  test  and  causes  errors  which  cannot  be  corrected. 

There  is  a  simple  test  to  determine  the  presence  of  thermal  emf.  It  requires 
that  a  motor  be  fired  without  taking  data.  For  exaiiq)le,  if  a  load  cell  cir¬ 
cuit  is  in  question,  the  motor  is  mounted  in  preparation  for  firing  (the 
longest  practical  duration  should  be  used).  A  sensitive  zero-center 
recorder  (e.g.,  -1  to  +1  millivolt) is  connected  to  the  signal  terminals, 
and  the  power  is  removed  from  the  transducer  input.  The  system  is  now  set 
up  to  record  any  undesirable  emf's.  The  recorder  input  terminals  are  then 
shorted  together  and  either  the  recorder  is  adjusted  to  zero  or  the  pen 
mark  is  noted  as  zero  input.  Then  the  short  is  removed.  Any  resulting 
movement  of  the  pen  will  show  initial  thermal  emf.  With  the  recorder  run¬ 
ning,  10  minutes  should  be  allowed  to  record  prerun  thermals,  then  the 
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motor  should  be  run,  waiting  10  to  15  minutes  afterward  to  allow  stabili¬ 
zation.  The  record  will  show  the  thermal  emf’s:  prerun,  run,  and  post¬ 
run.  This  test  will  be  indicative  of  the  presence  of  thermal  emf’s,  and 
should  not  be  used  for  correction  (thermal  errors  are  highly  nonpredictable) , 
but  should  be  used  to  determine  whether  redesign  is  necesseiry. 


System  Complexity 

Except  for  the  force  transducer,  the  error  elements  discussed  above  are 
caused  by  complexity;  complexity  is  in  turn  usually  caused  by  efforts  to 
provide  flexibility.  In  attempting  to  produce  a  measuring  system,  it  is 
advisable  to  remember  that  flexibility  and  precision  are  not  compatible, 
especially  in  low-voltage  circuits.  Ideally,  a  strain-gage  bridge  should 
be  coupled  to  its  conditioning  equipment  with  four  well-insulated  wires 
with  copper-to-copper  solder  connections  at  each  end.  This  is,  of  course, 
impossible  because  it  is  usually  necessary  to  have  zero-balance  controls 
and  simulation  devices.  However,  it  is  possible  to  eliminate  all  unneces¬ 
sary  interconnections  such  as  fanning  strips,  connectors,  relays,  and 
switches,  thereby  reducing  to  a  minimum  the  possibility  of  thermal  emf  and 
leakage  errors. 


PART  II:  SPECIFIC  APPLICATION  TO  EXISTING 
BATES  SYSTEM 


Thrust  Integration  System 

The  complete  data  system  as  installed  at  EAFB  is  impossible  to  analyze  for 
the  probability  of  error  because  there  is  no  way  of  determining  the  effects 
of  the  various  environments  to  which  it  is  exposed,  especially  heat,  humid¬ 
ity,  and  cleanliness.  However,  the  various  areas  in  which  error  could 
occur  can  be  discussed,  including  an  estimate  of  the  amount  of  error  which 
can  be  produced  in  each. 
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Calibration ,  A  mechanical  loading  calibrator  is  incorporated  in  the  system, 
No  detailed  study  of  this  unit  was  made,  but  in  general,  it  should  be  ade¬ 
quate  within  the  limits  imposed  by  the  load  cell.  Any  calibrator,  however, 
should  be  periodically  checked  for  accuracy  against  some  standard  (a  proving 
ring  should  be  suitable  for  this  case).  A  device  which  uses  mechanical 
linkages  is  susceptible  to  variation  with  use.  This  unit  should  be  accurate 
within  0.25^. 


Load  Cell .  This  unit  is  the  key  element  in  the  system,  eind  unfortunately  is 
also  the  most  susceptible  to  short-term  errors.  In  severe  cases,  the  load 
cell  can  be  as  much  as  ^  in  error  for  short  periods.  A  test  to  determine 
its  thermal  susceptibility  was  described  in  Part  I,  in  the  discussion  of 
thermal  emf. 


Interconnections .  Conversation  with  operating  personnel  indicated  the  pres¬ 
ence  of  four  or  five  connections.  These  could  be  reduced  to  two.  Thermal 
emf  errors  are  small  in  these  cases  (less  than  0.1^  in  any  case),  but  they 
cam  accumulate.  The  most  serious  problem  with  interconnections  is  the  pos¬ 
sible  presence  of  leakage  between  uninsulated  terminals,  and  this  is  often 
uncontrollable  because  temperature  fluctuations  cause  a  connector  housing 
to  "breathe"  humid  air,  thereby  depositing  moisture.  Errors  of  3^  have 
been  found  in  "waterproof"  connectors  which  were  not  properly  tested. 


Bridge  Completion  Unit.  Apparently,  the  equipment  furnished  by  B  &  F  was 
originally  designed  for  flexibility  to  avoid  the  extra  cost  of  an  individual 
design  for  each  problem.  It  is  usable  with  one-,  two-,  and  four-arm  bridges, 
and  is  supplied  with  the  necessary  boards  for  completing  each  bridge. 
Although  these  boards  are  within  a  rack  in  the  recording  center,  there  is 
still  a  possibility  of  some  leakage  error.  Although  the  possibility  is 
small,  such  flexibility  should  be  avoided. 

Zero  Balance  Control.  This  control  is  a  VDR  properly  located  between  the 
power  leads  with  a  limiting  resistor  to  isolate  the  control  from  the  bridge. 
However,  to  moke  this  unit  usable  with  bridges  of  different  resistances  and 
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sens itivites ,  the  isolating  resistor  is  5000  ohms.  This  value  vill  give 
reasonable  isolation,  but  the  control  becomes  extremely  sensitive.  For  a 
350-ohm,  3-nrv/volt  bridge,  this  control  gives  ±500^  variation  in  zero.  A 
resistor  of  100,000  to  250,000  ohms  would  be  adequate.  The  error  from  this 
unit  is  probably  negligible  unless  the  control  knob  was  accidentally 
touched  or  vibrated  after  calibration.  Wiper  motion  between  wires  could 
cause  appreciable  error  because  of  its  sensitivity. 

Resistance  Simulation  (*^R"Cal^0  Switch.  There  is  probably  no  error  in  this 
device,  but  two  factors  should  be  considered.  It  is  a  ”dry"  circuit  (i.e., 
low  current)  and  contact  resistance  should  be  monitored,  and  any  thermal 
gradients  between  the  switch  terminals  should  be  avoided.  Wiping  action  of 
the  contacts  should  keep  the  resistance  low,  and  there  does  not  seem  to  be 
any  heat-producing  element  near  the  switch. 


Voltage  Calibration .  The  switch  to  supply  calibrated  voltage  to  the  VCO 
input  is  subject  to  the  comments  of  the  foregoing  paragraph. 


Data  Recording  and  Reduction 

The  manner  in  which  the  data  are  recorded  and  reduced  can  affect  accuracy 
as  much  as  any  other  single  factor.  Figure  2  shows  the  beginning  and  end 
of  a  thrust-time  curve.  (The  steady-state  period  is  omitted.)  As  the 
system  is  operated,  a  gating  signal  is  used  to  turn  the  recorder  ’^on^’  when 
the  level  reaches  10^  of  design,  and  ’*off'*  when  it  returns  to  10^  at  cutoff. 
In  this  way,  the  force  indicated  from  the  igniter  is  not  recorded.  However, 
this  introduces  several  chances  for  further  error.  The  idea  is  predicated 
on  the  assumption  that  ignition  and  burnout  characteristics  are  the  same 
for  all  grains. 

With  ideal  conditions,  the  grain  surface  will  ignite  uniformly,  and  thrust 
will  build  up  rapidly,  leaving  only  a  small  area  (lined  shading  in  Fig.  2) 
before  the  gate  ”on’^  time;  but  if  ignition  spreads  slowly,  a  larger  thrust¬ 
time  area  is  lost  (dotted  shading).  Burnout  is  even  more  critical,  because 
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no  tvo  grains  will  burn  to  the  case  at  all  points  in  the  same  manner.  The 
lined  shading  again  shows  an  ideal  case,  but  because  there  are  always  small 
pockets  left  as  burnout  approaches,  it  would  be  expected  that  considerable 
energy  is  left  even  though  the  thrust  is  below  the  lOjt  level  as  shown  by 
the  dotted  area.  Thus,  a  biased  error  is  produced  which  is  variable  and 
consequently  uncorrectable. 

Another  error  which  can  occur  is  the  activation  level  of  the  gate.  This  is 
controlled  by  another  system  (transducer  plus  conditioning  equipment),  and 
its  use  could  easily  add  as  much  error  as  is  introduced  by  the  igniter.  The 
burnout  level  is  particularly  doubtful  because  the  heat  from  the  burning 
grain  can  produce  nonuniform  temperature  distribution  in  the  pressure  trans¬ 
ducer,  and  this  can  produce  not  only  thermal  emf's  but  will  not  permit  the 
temperature  compensation  to  operate  properly.  The  latter  is  predicated  on 
uniform  temperatures  throughout  the  sensing  element. 
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A  more  controllable  method  of  making  these  measurements  vould  be  to  test 
each  igniter  type  for  the  amount  of  impulse  it  adds  and  use  that  value  as  a 
correction.  The  total  impulse  could  then  be  recorded  and  the  igniter  impulse 
subtracted.  This  would  give  much  better  precision,  because  the  variable 
errors  would  be  reduced  as  much  as  possible  and  one  biased  but  measureable 
error  would  be  added.  The  correction  for  igniter  thrust  need  not  be  highly 
accurate  because  it  is  a  small  part  of  the  total. 

Without  analyzing  a  large  number  of  high-speed  analog  records  of  firings,  it 
is  impossible  to  determine  the  amount  and  variability  of  error  introduced  by 
gating  the  recorder.  However,  the  error  will  be  greater  for  small  grains 
■ftian  for  large  ones,  thereby  producing  an  untraceable  variation  in  specific 
impulse  (ig)-  Furthermore,  1^  is  always  in  error  in  the  negative  direction, 
which  prejudices  performance  with  respect  to  design;  i.e.  ,  the  system  can 
never  record  all  of  the  available  performance. 

An  element  which  might  have  been  a  factor  in  the  decision  to  use  the  lOjt  gat¬ 
ing  level  is  a-c  noise.  The  voltage-controlled  pulse  rate  generator  has  a 
characteristic  which  must  not  be  neglected;  it  does  not  respond  to  negative 
voltages.  Consequently,  if  pure  ac  is  fed  to  it,  it  will  count  only  the 
positive  portion  and  will  not  subtract  counts  for  the  negative,  so  that  the 
result  looks  like  a  positive  dc.  This  characteristic  means  that  if  the  com¬ 
plete  curve  is  integrated  and  if  a-c  noise  is  on  the  signal  line  at  the 
beginning  and  end  of  the  run  (near  zero  dc) ,  it  will  produce  an  error. 

There  are  three  different  ways  that  this  error  can  be  eliminated  without 
introducing  further  errors,  but  such  a  design  is  beyond  the  scope  of  this 
discussion. 
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INSTALLATION  REQUIREMENTS  OF  ROCKETDYNE  SYSTEM 


CONCRETE  TEST  PAD 

A  concrete  test  pad  of  sufficient  size,  strength,  and  flatness  is  necessary 
for  the  proper  support  of  the  Rocketdyne-designed  system.  The  complete  test 
stand  has  over-all  dimensions  of  approximately  20  feet  in  length,  5  feet  in 
vidth,  and  5  feet  in  height,  and  is  composed  of  two  in-line  sections  having 
approximately  13-  and  7-foot  lengths,  respectively.  The  sections  weigh 
about  18,000  and  4000  pounds,  respectively,  and  four  lifting  eye-bolts  are 
supplied  with  each  section.  All  points  on  the  sixrface  of  the  supporting 
section  of  concrete  should  be  within  l/4  inch  of  a  horizontal  surface  so 
that  alignment  and  leveling  problems  of  the  system  are  not  too  severe.  The 
design  and  configuration  of  the  two  existing  test  stand  mounting  beams,  24 
inches  apart  and  imbedded  in  the  concrete  of  the  BATES  test  pad,  are  con¬ 
sidered  satisfactory  for  the  Rocketd^me  system.  Consequently,  the  system  was 
designed  taking  into  consideration  an  identical  type  of  tiedown.  A  cross- 
sectional  area  of  an  imbedded  BATES  tiedown  beam  is  shown  in  Fig,  3. 

ELECTRICAL  POWER  AND  WIRING  REQUIREMENTS 

Satisfactory  operation  of  the  hydraulic  calibrator  and  the  strain-gage  load 
cell  during  calibration  and  measurement  operations  requires  an  adequate  quan¬ 
tity,  size,  and  configuration  of  wiring,  and  a  suitable  type  of  electrical 
power,  A  minimum  of  12  No.  18  wires  in  two  shielded  cables  of  6  wires  each 
are  required  for  operation  of  the  load  cell,  and  25  No.  18  wires  and  1  No.  8 
wire  are  required  for  operation  of  the  Ruska  hydraulic  calibrator. 

Electrical  power  required  at  the  test  stand,  from  an  EAFB  source,  is  10 
amperes  of  unregulated  24-vdc.  Regulated  d-c  power  for  the  operation  of 
the  load  cell  and  calibrator  will  be  supplied  from  rack-mounted  sources 
located  in  the  test  control  room  enclosure  supplied  by  Rocketdyno.  The*  12 
No.  18  wires  in  two  shielded  cables  of  6  wires  each  required  for  the  load 
coll  should  be  in  a  separate  conduit  from  all  other  wiring,  so  that  all 
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STEEL  BEAM 


Figure  3*  Cross  Section  of  Test  Stand  Tie-Dov7i  Beam 
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possible  electrical  interference  be  avoided.  All  instrumentation  virin^^, 
both  for  load  cell  and  calibrator,  should  pass  through  a  waterproof  junction 
box  located  in  the  ground  about  25  feet  from  the  test  pad.  The  junctions  in 
this  box  of  the  load  cell  wiring  should  be  the  only  breaks  in  the  load  cell 
cables  between  the  load  cell  connections  at  the  test  stand,  and  the  chassis 
connections  on  the  rack-mounted  units  in  the  control  room.  This  single  break 
in  each  load  cell  cable,  although  undesirable  from  a  measurement  accuracy 
standpoint,  is  needed  for  practical  rewiring  purposes  in  the  event  of  a  motor 
failure  resulting  in  the  destruction  of  wiring  in  the  vicinity  of  the  test 
pad.  All  cable  connectors  should  be  of  the  Bendix  waterproof  type.  Both  of 
the  cable  connections  at  the  load  cell  should  use  Bendix  socket  connectors 
No.  10-726l4-6S,  and  the  junction  box  connectors  in  each  load  cell  cable 
should  consist  of  Bendix  pin  connector  No.  10-lb9H^j-0I^  and  socket  con¬ 
nector  No.  10-72614-6S.  The  Ruska  calibrator  control  wiring  connectors  in 
the  waterproof  junction  box  should  consist  of  Bendix  pin  connector 
No.  10-72128-12P  and  socket  connector  No.  10-72628-125  for  the  25  No.  18 
wires,  and  pin  No.  10-72115-3P  and  socket  No.  10-72615-35  for  the  single 
No.  8  wire . 

GRAVITY  MEASUREMENT 

Accurate  total  impulse  measurements,  by  the  RocketdjTie  system,  require  cali¬ 
bration  of  the  measuring  system  by  the  hydraulic  calibrator.  Because  very 
accurate  dead  weights  will  be  used  to  produce  the  precise  hydraulic  pres¬ 
sures  utilized  to  obtain  the  calibrating  forces,  a  knowledge  of  the  specific 
value  of  the  acceleration  due  to  gravity,  at  the  test  site,  is  required  so 
that  the  weights  can  be  properly  fabricated.  In  accordance  with  these  con¬ 
siderations,  Rocketdyne  requested  and  received  from  EAFB  the  local  g  value, 
to  six  significant  figures,  at  Test  Pad  No.  5,  Solid  Motor  Test  Area, 

Rocket  Propulsion  Laboratory,  The  measurement  wais  performed  on  18  January 

2 

1963  and  was  determined  to  be  979*^77  cm,  sec  .  This  value  of  g  was  utilized 
in  the  design  of  the  calibrating  weights  by  the  Ruska  Instrument  Company, 
manuf actiirers  of  the  hydraulic  cal ibration  sys tern. 
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INTRODUCTORY  DISCUSSION 


Mathematical  analysis  is  an  essential  part  of  any  adequate  effort  to  design 
and  build  a  total  impulse  measurement  system  possessing  the  accuracy  required 
in  the  contractural  effort.  Without  this  analysis,  all  design  and  construc¬ 
tion  effort  would  be  based  solely  upon  qualitative  judgment  resulting  from  a 
combination  of  previous  trial-and-error  experience  and  guess-work.  Its 
application  ranges  from  its  use  to  determine  both  general  and  specific  guide¬ 
lines  during  the  system  design  phase  to  its  use  in  an  accurate  and  detailed 
evaluation  of  the  system  after  its  construction.  The  objective  of  all  mathe¬ 
matical  analysis  undertaken  in  this  contractual  effort  was  the  mathematical 
characterization  of  the  system  for  the  purpose  of  a  quantitative  determina¬ 
tion  of  its  behavior  in  the  performance  of  its  measuring  function.  The  anal¬ 
yses  presented  in  this  section  of  the  report  consist  not  only  of  static  equi¬ 
librium  considerations  needed  in  motor  suspension  flexure  system  design  and 
constructional  alignment  operations,  but  also  of  dynamic  response  determi¬ 
nations  of  the  mechanical  system  to  transient  force  excitations  similar  to 
thrust  outputs  encountered  in  rocket  motor  propellant  tests.  Accordingly, 
the  range  and  type  of  analyses  reflect  the  variety  of  inquiries  directed,  at 
various  times,  to  the  design  aspects  of  the  over-all  constructional  problem, 
and  the  sequence  of  the  following  mathematical  efforts  corresponds  to  their 
actual  chronological  occurrence. 

The  first  mathematical  effort  was  concerned  with  a  general  treatment  of  the 
relationchip  between  the  actual  and  measured  impulse  for  certain  linear 
dynamic  systems,  including  a  preliminary  discussion  of  the  impulse  errors 
incurred  by  less-than-inf inite  integrations  of  indicated  displacements.  No 
numerical  examples  were  undertaken  in  this  initial  effort  because  its  pur¬ 
pose  was  a  general  discussion  and  elaboration  of  the  basic  idea  underlying 
the  measuring  function.  The  second  mathematical  effort  was  concerned 
exclusively  with  trapezoidally  shaped  force-time  excitations  of  a  one 
spring-mass  system  viscously  damped.  Its  purpose  was  to  obtain  some  numer¬ 
ical  results  with  a  system-excitation  combination  which  approximated  an 
actual  thrust-measuring  system  and  rocket  motor.  The  results  of  this  under¬ 
taking  provided  the  initial  quantitative  indication  of  the  adequacy  of  the 
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mechanical  system  design.  The  third  mathematical  effort  vas  concerned  with 
a  more  accurate  and  detailed  study  of  the  actual  measurement  system,  con¬ 
sisting  of  a  series-connected  massive  concrete  block,  a  thrust-measuring 
load  cell,  and  the  motor  and  mount  assembly.  The  representation  of  the 
system  in  this  study  was  the  most  advanced  model  undertaken  in  the  Phase  I 
effort.  Its  purpose  was  the  determination  of  the  effects  on  the  impulse 
measurement  error  of  such  characteristics  as  concrete -block  compliance , 
variable  motor  mass  (propellant  burning),  nonviscous  damping,  and  nonlinear 
force -displacement  cal ibrations  because  of  load-cell  aging  and/or  imposed 
bending  moments.  The  results  of  this  advanced  study  provided  the  greater 
part  of  the  justification  of  the  design  employed  in  the  designed  system. 

The  last  of  the  mathematical  developments  contained  in  this  report  was 
undertaken  to  provide  an  analytical  approach  to  the  flexure  system  design 
of  the  motor  suspension.  The  result  of  this  effort  was  the  determination 
of  a  suitable  combination  of  flexure  system  and  load  cell  mechanical  stiff¬ 
nesses  that  effectively  isolated  the  load  cell  from  extraneous  forces  and 
moments  that  would  degrade  its  basically  linear  performance.  The  applica¬ 
tion  of  the  results  to  the  flexure  design  thereby  minimized  the  problem  of 
nonlinearities  in  the  system  calibration. 
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ANALYTICAL  DYNAMICS:  GENERAL  EXPRESSIONS 


INTRODUCTION  AND  SUMMARY 

Some  elementary  observations  are  presented  regarding  the  relationship  between 
the  integrated  input  and  integrated  output  of  certain  simple  types  of  linear 
dynamic  systems.  An  application  of  these  considerations  is  the  solid-propel¬ 
lant  total  impulse  measurement  system  designed  by  research  instrumentation 
personnel.  The  material  presented  constitutes  an  elaboration  of  some  com¬ 
ments  made  in  Ref.  1  on  the  relationship  between  the  total  applied  impulse 
and  total  measured  impulse  of  such  a  system. 

Three  basic  mathematical  models  are  considered  which  represent  increasing 
levels  of  generality:  a  simple  model  which  is  appropriate  for  initial  design 
studies  of  the  total  impulse  measurement  system,  and  two  more  general  models 
which  may  be  needed  in  subsequent  cmalyaes. 


Model  1 


Model  1,  a  damped  spring-mass  system,  is  shown  schematically  in  Fig.  4  and  is 
described  by  the  second-order  differential  equation 

m  ^  +kx(t)  =  f(t),t^0,  (l) 

dt  dt 

where  the  input  force  f(t)  =  0  for  t  <  0,  and  the  displacement  x(t)  of  the 
mass  m  is  subject  to  the  initial  conditions 


dt 

x(t) 


(2) 
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Model  2 


Model  2  is  a  system  described  by  an  nth  order  linear  constant-coefficient 
differential  equation: 


+  a^x(t)  =  f  (t)  ,  t  2:  0,  (3) 


vhere 

/  0,  f(t)  =  0  for  t  <  0,  and  the  displacement  x(t)  is  subject 
to  the  initial  conditions 
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Model  3 


Model  3  is  a  geoeral  linear  time-invaixant  stable  (physically  realizable) 
dynamic  systero  vith  impulse  response  (or  veighting)  function  W^(t)  and  fre¬ 
quency  response  function  Yj(w) ,  subject  to  zero  initial  conditions.  This 
model  includes  as  special  cases  systems  such  as  the  above,  vhose  input  force 
and  resulting  displacement  are  related  by  ordinary  linear  differential 
equations  vith  constant  coefficients.  Variable  coefficient  systems  are 
excluded  because  they  are  not  time-invariant. 


For  each  model,  vhat  is  of  interest  is  the  relationship  between  the  inte¬ 
grated  actual  input  force, 

=  J  £(t)dt  (!>) 

o 


and  the  integrated  displacement  indicated  by  the  spring, 

J  x(t)dt  (6) 

o 

For  the  case  of  an  idealized  total  impulse  measurement  system,  the  inte¬ 
grated  input  represents  the  total  impulse  generated  by  the  rocket  ujotur 
(the  ** total  applied  impulse")  and  the  integrated  output  is  a  measure  of 
the  total  impulse  measured  by  the  system  (the  ’’total  measured  impulse”). 
The  indicated  limits  of  integration  in  Eq.  5  and  6  extend  over  all  time 
0  ^  t  ®,  although  in  practice,  the  applied  force  f(t)  vanishes  for 
t  >  T,  where  T  is  the  duration  of  the  run,  and  x(t)  is  essentially  zero 
after  some  later  time. 


For  the  three  models  under  consideration,  the  following  is  shown: 

Model  1:  =  K,l.  +  K, (me,  +  be.) 

M  1  A  1^2  1 

Model  2: 

M  2  A  2  n  D  n-1  n-1  1  1 

Model  3: 


(7) 
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where  the  quantities  K^,  and  Kj  appearing  in  Eq,  7  are  system  constants 
given  by 


K, 


k 


K 


2 


K 


3 


W(t)dt 


(8) 


These  results  are  valid  provided  that  the  input  f(t)  is  a  sufficiently  veil- 
behaved  function  emd  each  of  the  systems  is  stable  (Ref .  2) .  The  conditions 
required  of  the  input  will  undoubtedly  always  be  satisfied  by  thrust 
functions . 


RESULTS 

Before  deriving  the  results  stated  in  Eq.  7  and  8,  some  alternative  formu¬ 
las  are  presented  for  the  constants  K^,  Kg,  and  and  some  of  their  prac¬ 
tical  implications  are  discussed.  It  should  be  noted  first  that  under 
zero  initial  conditions  the  terms  enclosed  with  parentheses  in  Eq.  7  drop 
out,  and  for  each  of  the  three  systems  the  relationship  between  the  total 
measured  impulse  and  total  applied  impulse  is  of  the  form 

I„/I^  -  t  (9) 

where  K  is  a  constant  determined  by  the  system.  Equation  9  implies  that 
whenever  the  system  is  excited  under  zero  initial  conditions  by  an  input 
f(t),  the  ratio  of  the  total  measured  impulse  to  the  total  applied  impulse 
is  a  constant  which  does  not  depend  on  the  specific  form  of  the  input.  More¬ 
over,  this  relationship  indicates  one  way  in  which  the  constant  K  may,  in 
principle,  be  determined  experimentally.  If  possible,  an  input  with  known 
total  impulse  I^  would  be  applied  (the  pointvise  behavior  of  the  input  need 
not  be  known),  Ij^  would  be  measured,  and  the  ratio  would  be  computed. 

Such  an  experiment,  which  may  be  described  as  a  total  impulse  calibration 
procedure,  would  be  repeated  with  sufficient  frequency  to  obtain  a 
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statistically  reliable  estimate  of  the  constant  K.  Measurement  analysis 
of  the  ’’constant"  K  could  then  bo  performed  in  the  usual  manner. 


A  second  way  (as  an  alternative  to  Eq.  8)  can  now  be  considered,  in  which 
the  constants  and  in  Eq.  7  may  be  interpreted.  This  will  lead 

to  another  characterization  which  has  possible  practical  applications.  The 
three  systems  under  consideration  possess  frequency  response  functions 
^,(^41  and  y,(t4  which  determine  the  effect  of  the  systems  on  the 

1  ^  j 

various  frequency  components  of  the  inputs.  For  models  1  and  2,  the  fre¬ 
quency  response  functions  and  may  be  expressed  explicitly  in 

terms  of  the  coefficients  in  the  differential  equations  describing  the  sys¬ 
tems.  In  fact,  the  following  is  true  (Ref,  2): 


Y  i - 

^  m(ito)^  -e  b(ic4  k 


From  Eq.  8  and  10,  it  follows  that 


Vi(o) 

Ygfo) 


Ine  constant  may  be  written  in  a  similar  form.  For  if  Y^(c4  denotes 
the  frequency  response  function  of  system  3i  then  in  terms  of  the  system 
impulse  response  function  W^(t), 

Y^(i4  -  I  w^(t)e""‘^^dt,  (] 

0 

and  evaluating  Eq.  12  at  zero  frequency  (u)  =  O)  gives 


Y3(o) 


CO 

I  VI  (t), 


Hence,  from  Eq.  8  and  13, 


K3  =  Y^(0)  (14) 

Each  of  the  constants  appearing  in  Eq.  7  is  therefore  of  the  form 

K  =  Y(0)  ,  (15) 

where  Y(o)  is  the  system  frequency  response  function  Y(co)  evaluated  at 
zero  frequency. 

The  above  discussion  also  indicates  a  third  way  of  describing  the  constants 
Kg,  and  K^.  All  three  systems  possess  impulse  response  fun  ctions  which 
are  related  to  their  frequency  response  functions  in  the  manner 

Y(u;)  =  f  dt  (16) 

0 

That  is,  the  frequency  response  function  is  the  Fourier  transform  of  the 
impulse  response  function.  For  cc  =  0  in  Eq,  l6 

CD 

y(o)  =  J  W(t)dt  (17') 

0 

Thus,  from  Eq.  15  and  17  it  can  be  seen  that  each  of  the  constajats  in 
Eq.  7  equals  the  total  time  integral  of  the  system's  impulse  response 

function, 
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where  W(t)  is  the  system  impulse  response  function.  Allowing  T  to  become 
arbitrarily  large  in  Eq.  19  gives 

lim  U(t)  =  \  W(t)dt  .  (20) 

r  -  CD 

From  Eq.  18  and  20  it  follows  that  the  constants  K^,  Kg,  and  may  be 
described  es  the  steady-state  values  of  the  responses  of  the  systems  to 
unit  step  inputs.  This  lost  characterization  of  the  constants  appearing  in 
Eq.  7  indicates  another  way  in  which  they  may,  in  principle,  be  determined 
experimentally.  A  step  input  would  be  applied  to  the  system  (subject  to 
zero  initial  conditions,  for  example)  and  the  input  would  be  maintained  at 
a  constant  level  until  the  measured  response  became  essentially  constant. 

If  the  step  has  height  A  1 ,  then  the  value  of  the  steady-state  response 
is  to  be  divided  by  A  to  obtain  the  constant  K.  Again,  such  a  procedure 
would  be  repeated  sufficiently  to  form  a  statistically  reliable  estimate  of 
the  constant  K  and  its  variability. 


It  should  be  noted  in  passing  that  if  it  is  not  possible  to  apply  a  true 
step  input  to  the  system  but  only  a  ramp  function 


0  for  t  0 

J  atforO^t^t 

o 

A  for  t  ^  t 
L  0 


then  the  above  procedure  (measuring  the  steady-state  value  of  the  response 
as  an  estimate  of  the  constant  K)  may  still  be  applied.  A  step  input  is 
preferable,  however,  because  it  can  be  shown  that  the  experimental  error 
arising  from  not  being  able  to  pass  to  the  limit  T  ®,but  instead  taking 
as  an  estimate  of  K  the  response  x(t^)  at  some  large  time  smaller 

for  a  step  input  than  for  a  corresponding  ramp  input.  The  difference  in 
the  errors  for  the  two  types  of  input  becomes  small  as  increases,  how¬ 
ever,  and  therefore  should  provide  no  practical  difficulties  in  the  appli¬ 
cation  of  this  method. 
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miVATION  OF  RESULTS 


In  this  discussion,  the  results  expressed  in  Eq.  7  ftJQd  8  of  the  introduction 
are  demonstrated.  Model  5  is  used  first  because  the  result  for  this  case 
is  almost  a  matter  of  definition  and  also  is  used  in  the  other  two  cases. 


Model  3 


A  linear  time-invariant  stable  dynamic  system  is  characterized  by  its 
impulse  response  function  W(t),  which  represents  the  response  of  the  system 
to  a  unit-impulse  applied  at  time  t  =  0.  For  such  a  system,  the  response 
x(t)  to  an  input  f(t),  with  zero  initial  conditions ^  is  given  in  terms  of 
the  impulse  response  function  W(t)  by 

x(t)  =  J  W(T)f(t  ~  T)dT  (21) 

0 


Here  it  is  assumed  for  simplicity  that  the  system  is  physically  realizable 
(the  response  depends  only  on  the  past  of  the  input  and  not  on  its  future) 
and  that  the  input  vanishes  for  t  <  0.  It  is  also  assumed  that  the  system 
possesses  a  frequency'  response  function 


y(uJ)  =  r  w(t)e  dt 

•^0 


and  that  the  input  f(t)  has  a  Fourier  transform^ 


j  f(t)e“^ 

*'  A 


-itet 


dt 


(22) 


♦For  example,  a  function  f  i  defined  for  t  ^  0,  possesses  a  Fourier  trans¬ 
form  if  is  continuous  (except  possibly  for  a  finite  number  of  jumps) 

and  if  Jol'Wl  dt  exists.  To  avoid  complications,  it  is  assumed  from 
this  point  that  these  conditions  ore  met  by  the  inputs.  Note  that  the 
latter  requirement  follows  automatically  from  the  first  if  f(t)  vanishes 
for  all  t  greater  than  some  T  >  0,  which  will  certainly  be  true  for  all 
inputs  to  the  total  impulse  measurement  system. 
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Then,  if  F  (uj)  denotes  the  Fourier  transform  of  the  output,  the  frequency 

X 

domain  relationship  (corresponding  to  £q.  2l)  is 


Nov,  if  ,  y(w),  and  Fj(u^  are  evaluated  at  to  «  0, 


F^(0) 

y(o) 

Ff(0) 


J  3c(t)d 

i; 

00 

J  f(t) 


dt 


dt 


(23) 


‘A 


(24) 


Hftcce,  from  Eq.  23  and  24 

00 

\  W(t)dt  ]l. 


M 


which  is  the  desired  result.  It  can  be  observed  again  that 


J  w{t)dt 


y(o) 


(25) 


Model  1 

Considering  now  the  linear  second-order  constant  coefficient  differential 
equation,  1,  with  initial  conditions  Eq.  2;  every  solution  x(t)  of  Eq,  1 
for  t  ^  0  is  of  the  form 


x(t)  =  h(t)  +  y(t) 

where  h(t)  is  the  solution  of  the  homogeneous  equation 

n,  ^  b  +  kh(t)  =  0,  t  ^  0 


(26) 


dt 


(27) 
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satisfying  the  initial  conditions  dh(t)/dt  *  Cg  and  h(t)  -  at  t  =  0, 
and  where  y(t)  is  a  particuleir  solution  of  the  nonhomogene ous  equation 
(Eq.  l)  with  zero  initial  conditions  for  y(t) .  Here  it  is  assumed  that  the 
input  is  sufficiently  well-behaved  so  that  it  possesses  a  Fourier  transform 
given  by  Eq.  22,  Corresponding  to  the  decomposition  (Eq.  26),  the  total 
measured  impulse  may  be  split  into  two  parts: 


■  ^ 


where 


j'hW. 


is  the  contribution  to  arising  from  h(t)  and 

00 

^Hy  -  J  (50) 

is  the  contribution  to  arising  from  y(t).  Although  for  the  total  impulse 
measurement  system  it  seems  appropriate  to  assume  zero  initial  conditions 
(in  which  case  h(t)  is  identically  zero  and  =  l^j  )»  the  expression 
(Eq.  7)  for  which  is  obtained  i  valid  in  the  general  case. 

Considering  first  the  contribution  1^^  to  which  is  associated  with  the 
given  initial  conditions,  the  solution  h(t)  of  Eq.  27  may  easily  be 
expressed  in  explicit  form  in  terms  of  the  parameters  m,  b,  and  k  (for 
example,  in  Ref.  2),  but  it  is  not  necessary  to  deal  with  such  explicit 
solutions  to  derive  the  result  (Eq.  7)-  It  need  only  be  shown  that  by 
virtue  of  the  assumption  that  the  system  is  stable  (a  condition  ensured  by 
the  assumption  b/m  >  O),  the  solution  h^t)  as  well  as  its  first  derivative 
dh(t)/dt  approaches  zero  as  t  becomes  arbitrarily  large. 


With  this  in  mind,  Eq.  27  is  integrated  over  the  interval  0  ^  t  ^  T  to 


obtain 


“  J  ^  dt  +  b  J  dt  +  k  J  h(t)dt  =  0 


0  dt 


’0  dt 
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and  therefore 


(nih'(r)  +  b.b(T)  )  -  (mh'(o)  +  bb(0)  )  +  k  J  h(t)dt  «  0 

0 

where  dh(t)/dt  has  been  denoted  by  h'(t).  Thus 

J  h(t)dt  =  "  (inh'(O)  +  bh(0)  )  -  (mh'(T)  +  bh(r)  )  (31) 

Now,  because  and  h(T)  approach  zero  as  T  ^  Eq.  31  becomes,  for 

arbitrarily  large  T, 

h(t)dt  =  i  (mcg  +  bCj)  (32) 

where  the  initial  values  Cg  and  c^^  have  been  substituted  for  h’(o)  and 
b(o).  As  anticipated,  =  0  for  zero  initial  conditions  (c^  =  Cj  =  O). 

Considering  the  contribution  to  arising  from  y(t) ,  because  system  1 
is  a  linear  time -invariant  stable  dynamic  system  with  frequency  response 
function 


klj^  -  (mCg  +  bCj) 


(36) 


Before  continuing  to  model  2  consider^  briefly^  the  errors  which  result 
from  not  being  able  to  perform  the  integration 


r  ® 

x(t)dt  -  lyi  (3 

over  the  infinite  range  0  ^  t  i  The  situation  will  be  restricted  to 
where  f  (t)  vanishes  outside  of  an  interval  0  ^  t  T  (for  total  impulse 
measurement  T  is  the  full  duration  of  the  run).  Similar  remarks  will 
apply  to  the  more  general  case  of  model  2. 

Supposing  the  response  x(t)  is  known  only  over  the  time  interval 

,  where  t  >  T,  then  what  is  of  interest  is  the  error 
o’  o  ’ 

*M  ■  ‘m  -  1° 

r^o 

made  in  approximating  by  x(t)dt  and  also  the  error  made  in 

M  j  0  r^o  ^ 

determining  from  Eq.  36  when  x(t)dt  is  used  in  place  of  The 

latter  error  ia  given  by 


and  ia  easily  seen  by  subtracting 


j*t 

k  x(t)dt  -  (mCg  +  bCj) 


frost  the  right  side  of  £q«  36, 


from  Eq.  37 


I  ana 

'V  “  Jo  *^*)‘*^  “  Jo  *<^-^‘**  '  Jt 

O 


Because  f(t)  =  0  for  t  ^  x(t)  is  solution  of  the  homogeneous 
equation 

m  — 1  b  4  kx(i)  -  0,  t  ^  t  (41) 

dt"  dt  ® 

The  problem  of  detennining  6  as  given  by  Eq.  40  from  the  differential 
equation  (Eq.  4l)  is  entirely  similar  to  the  problem  of  determining 
from  Eq.  27.  Here  the  initial  time  t  =  0  is  replaced  by  and  the 
initial  conditions  dh(t)/dt  =  H(t)  =  c^  at  t  =  0  are  replaced  by 
dx(t)/dt  =  x(t)  =  x(t^)  at  t  =  t^. 

Reasoning  as  before,  the  expression 

i«  obtained  Instead  of  that  given  by  Eq.  32.  Also,  from  Eq.  39  and  42, 


6.  =  mx'(t  )  +  bx(t  ) 

A  '  o  'o' 


It  can  be  observed  further  that  by  virtue  of  Eq.  3^,  38t  «wd  42 

P  t 

I.  ^  k  1  ®  x(t)di  +  (nac’(t  )  +  bx(t  )  )  -  (ibx’(O)  +  bx 

A  V  Q  0  0 


(0)  ) 


or,  if  zero  initial  conditions  x'(0)  =  x(0)  =  0  are  assumed. 


r 

^A  "  Jo  i(t)dt  +  (mx*(t^)  +  l»x(t^)  ) 


The  quantities  |  x(t)dt  and  x(t  )  which  appear  in  £q.  42,  43i  and  43 
^(3  o 

are  directly  measurable  from  recorded  data,  and  presumably  the  parameters 
m,  b,  and  k  would  also  be  known.  If  x^(t^)  could  be  estimated  from  given 
data,  or  if  *'(t )  were  known  to  be  nearly  zero,  then  in  principle  6^, 

6^,  and  could  be  determiued  from  Eq.  42,  43,  and  45.  It  should  be 
noted  that  under  the  assumptions  made  in  this  section,  x' and 
will  approach  zero  for  t^  . 
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Model  2 


The  proof  of  the  result  (eq.  ?)  for  system  2  may  be  carried  out  in  a  manner 
entirely  similar  to  that  for  system  1.  The  solution  x(t)  of  Eq.  vith 
initial  conditions  (Eq.  k) ,  is  represented  as  the  sum  of  solution  h(t)  of 
the  homogeneous  equation  (vith  initial  conditions,  (Eq.  4)  )  and  the  solu¬ 
tion  y(t)  of  the  original  nonhomogeneous  equation  (Eq.  3)  having  zero 
initial  conditions.  Again,  splitting  into  tvo  parts. 


The  result  for  model  3  is  applied  to  conclude  that 


^(O)  \ 


where  Y2(0)  is  the  frequency  response  function  of  system  2  evaluated  at 
W=  0.  Because,  as  noted  earlier, 


Yo(w)  = 


^(iW)”  +  a^_j(ia))"“^  +  ...  +  aj(ia})  + 


Then,  Y  (o)  =  l/a  and  therefore 


=  —  I, 

My  a  A 

^  0 


To  determine  the  homogeneous  equation 


i^h^"\t)  +  a^_jh^“"^^(t)  -t-  ...  +  a^h(t)  = 


is  integrated  over  the  interval  0  ^  t  ^  T  to  obtain 

r  Y  1 

J  h(t)dt  =  —  (a^c^  +  a^_^  c^_j  +  ...  +  ajCj) 
0 


.  +  a 
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n.  If  this 


where  =  d^hCt^df^  and  c  .  =  h^'^“^^(o);  j  =  1,  2,  .  . ., 

J 

system  is  stable  in  the  dynamic  sense,*  then  as  T  ®  the  terms  within  the 
second  pair  of  parentheses  approach  zero,  and 


IvfK  =  f  h(t)dt  =  —  (a  c 
Mh  a  '  n  n 


+  a  ,  c  , 
n-1  n-1 


+  a 


(47) 


Thus ,  from  Eq 


46  and  47  the  desired  result  is  obtained: 

1  T  1 

I.  + 

a  A  a 
0  0 


lac  +  a  -  c  - 
n  n  n-1  n-1 


ajCj) 


For  this  purpose  it  is  sufficient  to  assume  that  the  (generally  complex) 
roots  of  a  +  ...+a^X+a  =  0  have  negative  real  parts. 

For  model  ?,  which  corresponds  to  the  case  n  =  2  in  model  2,  the  assumption 
b/g  ^  ^  which  was  made  is  equivalent  to  the  assumption  that  the  roots  of 
mX  ^  bX  -f*  k  =  0  have  negative  real  parts. 
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ANALYTICAL  DYNAMICS: 


SPECIFIC  LXPIIFISSIO.VS 


FREE  VIBRATION  ^vTTB  VISCOUS  DAMPING 


With  reference  to  Fig.  5,  a  spring- suspended  mass,  m,  is  considered  free  to 
vibrate  vertically  in  a  medium  offering  a  viscous  resistance,  i.e.,  a  resis¬ 
ting  force  proportional  to  the  first  power  of  the  velocity.  Such  resistance 
is  usually  termed  viscous  damping.  The  mass  of  the  spring  itself  is  con¬ 
sidered  small  compared  with  the  suspended  mass,  so  that  the  discussion  is 
concerned  with  the  motion  of  a  single  particle.  In  this  case,  if  the  mass, 
m,  be  displaced  any  distance,  x,  (within  the  elastic  range  of  the  spring) 
from  its  position  of  equilibriiun  and  then  released,  it  will  begin  to  move 
under  the  action  of  a  resultant  force. 

X  =  -kac  -  ci  (^8) 

where 

k  =  spring  constant  =  force  necessary  to  produce  unit  extension  of  the  sprin 
c  =  coefficient  of  damping  =  magnitude  of  the  viscous  resistance  at  unit 
velocity. 

Newton’s  second  law  of  motion  states  that  when  a  force  acts  upon  a  given 
particle,  that  force  produces  an  acceleration  of  the  particle  in  the  direction 
of  the  force  and  proportional  to  its  magnitude.  For  the  force  expressed  by 
Eq.48y  the  equation  of  motion  becomes: 


mx  =  -kx  -  ci 

Dividing  through  by  m  and  introducing  the  notations 


^  «*  u?  and  ~  *•  2n 
m  m 


the  equation  may  be  written  in  the  fora 
X*  ^  2nx  +  <J?X  =  0 


(49) 

(50) 

(51) 
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EQUILIBRIUM 
POSITION 
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X 

1 

Fi^zure  '5.  Simple  Spring-Mass  System,  Free  Vibration 

Thus,  for  free  vibrations  with  viscous  damping  a  linear  differential  equation 
with  constant  coefficients  is  obtained. 

As  a  trial  solution  of  Eq.  31 9  take 


X  =»  Ce 


rt 


(52) 


where 

r  a  constant. 


Substituting  this  trial  function  for  %  in  Eq,31  results  in  the  auxiliaty 


+  2nr  +  uF  =  0 


(5-^) 


This  oquation  dotennincs  two  values  of  r  for  which  expression (‘52)  can  satisfy 
Eq.  51: 


(54) 


Hence,  the  general  solution  of  Eq, 51  becomes: 


where 


A  and  B  =  arbiiraiy  constants  of  integration. 

To  attach  any  physical  significance  to  this  solution,  two  distinct  cases 

r? — ? 

must  be  di3tin(]:ui3hcd,  depending  on  whether  the  radical  Jn  -  ^  in 
EqvK\tioii  ia  real  or  imaginary,  in  other  words  whether  n>coor  n  <  U>. 

Lqvuition  ^0  imiicales  that  this  depends  on  the  relative  magnitudes  of  the 
daipping  coefficient,  c  and  the  spring  constant,  k.  Generally  speaking,  a 
large  coefficient  of  damping  and  a  small  spring  constant  will  result  in 
real  values  of  r^  and  while,  for  the  reverse  of  these  conditions,  r^ 
and  will  be  complex  numbers. 


In  connection  with  this  solution,  it  should  be  noted  that  both  r^  and  Tg 

are  negative,  being  numerically  the  larger.  Thus  the  displacement  x 

has  the  some  sign  os  x^  and  approaches  zero  as  a  limit  when  t  becomes  infinitely 


"^7 


large.  The  motion  is  not  a  vibration  at  all  but  simply  one  in  which  the 
suspended  mass,  after  its  initial  displacement,  gradually  creeps  back  toward 
the  equilibrium  position  but  takes  theoretically  infinite  time  to  get  there. 
It  results  from  the  fact  that  the  damping  coefficient  is  too  large  compared 
with  the  spring  constant  and  is  sometimes  of  practical  interest  in  con¬ 
nection  with  certain  types  of  electrical  measuring  instruments  such  i  the 
galvanometer.  In  the  case  where  n  =  W,  aperiodic  motion  is  also  o^^iincd, 
and  the  corresponding  value  of  the  damping  coefficient  c  =  2nm  becomes 


'crit 


2  (dm 


and  is  called  critical  damping. 


Case  2,  n  <  (0 


More  often,  the  case  is  encountered  where  n  <  CO  so  that  the  roots  r^  and 
are  complex.  In  discussing  this  case,  it  will  be  convenient  to  change  the 
form  of  Eq.'35  to  indicate  more  clearly  its  physical  significance.  Remembering 
that  the  concern  now  is  only  with  the  case  where  n  <  (o,  let 

(J^  -  n^  =  (60) 

Then  expressions  in  Hq.  5^*  become 


rj  =  -n  +  iW^ 


Tj  .  -n  -  iW, 


where 


i  =  ./-I 


and  Eq.55  nay  be  written  in  the  form 


-nt  /.  iW  t  -it*J,tv 

X  =  e  (Ac  1  +  Be  1  ) 
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Using  the  knovn  relations 


e  *  cos  z  +  1  sin  z 


-iz 

e  3  cos  z  -  1  sin  z 


Eq. 55,  in  turn,  may  be  written  in  the  form 


=  0  (Cj^  cos  Wjt  +  Cg  sin  Wj^t) 


where 


and  Cg  =  new  arbitrary  constants. 


To  evaluate  these  constants,  assume  now  the  initial  conditions 


X  =  Xq  ,  X  =  Xq,  when  t  =  0 


Then  substituting  Eq.  65  into  Eq.  6^,  together  with  its  first 
derivative  with  respect  to  time,  there  is  obtained 

Xq  nx^j 

Cl  =  Xq  and  Cg  - 

and  the  solution  becomes 


^Xq  cos  Wjt  +  C^O  ^  ^0^  sin  <*^tj 


The  general  solution  (64)  can  be  represented  in  another  form  by  using  the 


relation 


cos  (e  -<p)  3  COS0  cos^  +  sin0  ein(p 


to  obtain 


X  =  Ae~°^  cos  (Wjt  -  a) 
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vhtre 


A  -  ♦  C 


<X  >  tan 


2  TT 

2 

-1  ^2  foras  of  the  Arbitrary  constants « 


Using  the  values  of  and  from  expressions  (6o)  these  new  constants  become 


and 


+  nx. 


\ 


a  =  tan 


■if  *0  A 

1 


(70) 


(71) 


A  displacement-time  curve  plotted  from  Eq.69  is  seen  to  be  vibratory  in 
nature  and  represents  so-callcd  damped  free  vibrations.  Equation  (67)  is 
therefore  the  expression  for  the  vibratory  effects  of  any  initial  displacement 
or  initial  velocity  x^,  when  t  ==  0. 


FORCED  VimUTION,  GFSUViL  DISTURBING  FCSICE 


l*/hile  damped  free  vibrations  arc  important,  there  are  many  ruses  for  which 
thz  diiferentiai  equation  of  motion(5l)  of  a  spring-suspended  laass  iu  Fig.  6a 
takes  the  more  general  forxM 


V  +  2ai  -♦"uFx~Q  =  f(t) 


(?2) 


where 


q  - 


the  disturbing  force  per  unit  of  suspended  mass, 


In  dealing  with  this  equation,  a  procedure  aomewhat  different  from  that  used 
previously  will  be  followed.  Referring  to  Fig.  6b,  it  is  f^ssumed  that  the 
disturbing  force,  q,  per  unit  of  mass  is 'given  by  a  curve  AB,  At  any  instant, 
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Figur^^  6a.  Spring-Mass  System, 
Forced  Vibration 


t*,  consider  one  elemental  impulse  q  dt’  as  represented  by  the  shaded  strip  in 
the  diagram.  This  one  impulse  imparts  to  each  unit  of  rtiass  an  instantaneous 
increase  in  velocity  dx  =  q  dt*,  regardless  of  what  other  foices,  such  as  the 
spring  force,  may  be  acting  on  it,  and  regardless  of  its  displacement  and 
velocity  at  the  instant  t*.  Treating  this  increment  of  velocity  as  if  it 
were  an  initial  velocity  (at  the  instant  t’)  and  using  Eq.  ^7  »  ii-  is  concluded 
that  the  corresponding  dispiacoment  of  the  spring -suspended  mass  at  any 
later  time,  t,  will  be 


dx 


(71) 


Since  each  impulse  q  dt*  botw^een  t*  «  0  and  1*  =  t  has  a  like  effect,  there  itt 
obtained,  as  a  result  of  the  continuous  action  of  the  disturbing  force  q,  the 
total  di8placemea\u 


i  r  -n(t-t') 

=  —  !  qe  ' 


id. 


&rn 


<it' 


) 


This  expression  still  docs  not  include  the  effect  of  any  initial  displacement 
or  initial  velocity  wJien  i  ~  0.  These  effects,  however,  are  exactly 


^li 


those  represented  by  Eq*67  ;  hence,  for  a  complete  solution  of  tq*  72*  there 
is  obtained 


At  first  glance  it  night  appear  that  the  first  tera  in  this  expression 
represented  free  vibratiorxs  and  the  laot  term  forced  vibrations,  but  this 
is  not  quite  the  case.  Actually  the  first  ternK  represents  only  the  effect 
of  initial  displacement  and  initial  velocity;  the  last  torn  represents  the 
complete  effect  of  the  disturbing  force  q  f(t’).  The  disturbing  force 
produces,  on  its  ora  account,  both  free  and  forced  vibrations  and  all  this 
together  is  accounted  for  by  the  integral  in  Eq.  75  ,  i,e.,  by  the  solution 
(74),  For  this  reason,  Eq.  74  and  75  arc  especially  useful  in  studying  the 
early  effects  of  a  disturbing  force  before  damping  has  had  timj  to  dissipate 
the  free  vibrations. 

For  convenience  in  further  calculations,  jt  will  be  helpful  to  transform  the 
integi/al  in  Eq.  75  as  follows:  Using  the  relation 

sJin  ("^4^t  ~  Wjt’  -  cos  bJ^t  sin  W^t*  (7^^) 


The  complete  solution  is  given  as 


X  ^  e 


-nt 


[xg  cos  wt  +C  —  )ain  W^t] 


“■Tit  /  t  A  A 

e  8 1  n  Cil  t  r  t 


to 


1:  f 

Jo 


cos  dt' 


/.I  * 

<■  COB  W.t 


'»t  ft 

—  qe”  sin  dt* 


(77) 


Application  of  this  expression  to  rectangular  and  trapezoidal  force-tiise 
inputs  will  DOW  be  Ktade. 


42 


APPLICATION  OF  METHOD  TO  RECTANGULAR  AND  TPjVPEZOIDAL  FOPXE-TIHE  INPUTS 


RECTANGULAR  FCStCE-TTME  INPUT 

Application  of  Eq*  77  for  a  rectangular  force-tLne  input  to  the  spring-aaaa 
system  provides  a  relatively  simple  demonstration  ol  the  mathematical  method 
and  the  results  obtained  by  its  use«  Althougii  stepped  inputs  are  not 
encountered  in  actual  rocket  engine  performance  evaluati?>n8 ,  they  are  em¬ 
ployed  for  test  stand  dynamic  evaluation  tests  by  the  use  of  fracture  links 
and  snapped  wires.  Figure  7  shows  the  shape  of  the  input  force-time  curve. 
The  constant-level  force  duration  is  arbitrarily  taken  to  be  great  enough  to 
dissipate  essentially  all  oscillation  of  the  system  produced  during  the  up- 
step  by  the  damping  (always  present  in  some  degree)  by  the  time  the  down-step 
occurs  at  the  input  force  termination.  Separate  evaluations  of  Eq.  77  for 
the  two  steps  arc  made  because  of  the  difference  in  initial  conditions  and 
load  existing  at  the  instants  the  two  steps  occur. 


0 


Figure  ?•  lorce-Timi*  Jlec tangular 


Case  ] .  Stepped  Loading 


0  . 


Initial  conditions  nt  t  0  are  =  0,  =  0,  and  q  =  —  is  constant  during 

the  loaded  phase.  The  subsequent  displacement,  x,  of  the  mass,  m,  using 
Eq.  77  is 


X  = 


F  c""* 

raWj 


[“"  "i‘  Jo 


^  nt'  r*'  lit'  1 

e  cos  Wjt'dt'  -  cos  tdjt  J  e  sin  dt'j 


from  which  is  obtained 

-nt 


(n  ) 


(78) 


(79) 
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Multiplying  this  displacement  equation  by  the  spring  constaEit,  k,  and 
utilizing  the  relations 


lu  -  F  , 


u?  .  IL 


'  _2  ,.,2 


vhere 


F  =  force  indicated  by  the  spring 


there  re  ults 


F  =  Fq  (1 


^  e  sin  w  t 

I  ^ 


By  integration  of  Eq.  81  for  the  duration  of  the  applied  step  loading  an 
expression  is  obtained  relating  the  impulses  of  the  actual  input  force 
and  the  force  measured  by  the  spring. 

Jo**' ■  fflljo*  ^  lo**'"*  ““  “i‘  -  lo  “•  “i*  “‘J 


Since  the  actual  impulse  in  this  case  is  using  the 

notation  for  measured  impulse,  =  Jg  Fdt,  provides  the  expression 


n  j 


o  rk 

/  A  A 


Equation  83  indicates  that  the  difference  between  the  actual  and  measured 
impulses  (error);  for  the  duration  of  the  stepped  loading,  is  a  fimction 
of  the  force  level,  the  damping,  the  damped  nacui'al  frequency  mid  the  force 
duration  time. 

Case  2,  St(p  >ped  Italoading 


Taking  t  >■  0  at  t^e  instant  of  system  unloading  the  initial  conditions  are 
taken  to  be  Xg  *  «  0  and  q  «  0,  Although  a  zero  initial  velocity, 

at  the  stated  displacementi  would  theoretically  require  an  infinite  time  to 
obtain,  by  damping,  subsequent  to  an  upstep,  it  is  en  irely  realistic  and 
practical  to  regard  it  obtained  in  a  finite  time*  Fox  example,  examination 
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of  the  derivative  of  Eq.  79  ahovs  that  i  ^  10  in/sec  in  2  sccondi  after 
a  10,000-po\md  stepped  load  is  applied  to  a  175  cps  systen  damped  only 
3  percent  of  critical  damping. 

The  subsequent  displacement,  x,  of  mass,  m,  again  using  Eq.  77  in  which 
these  values  arc  substituted,  is  seen  to  be: 

(84) 


and  since  =»  0  for  periods  subsequent  to  the  stepped  imloading,  Eq.  85 
represents  the  error  incurred  during  the  time  t  after  the  down-step.  By 
inspection  of  Eq.  83  and  83  it  is  seen  that  the  two  impulse  errors  asso¬ 
ciated  with  a  rectangular  force-time  input  are  equal  in  magnitude  and 
opposite  in  sign,  for  equal  integration  times  subsequent  to  the  two  steps, 
and  result  in  a  zero  net  overall  measurement  error  by  the  system  due  to 
mechanical  response.  The  only  restrictions  incorporated  in  the  foregoing 
application  were  (l)  a  constant-level  force  diiration  long  cnox'gh  toebtain  essen¬ 
tially  static  equilibrium  by  the  time  of  the  do\*Ti-step  and,  (2)  an  equal 
integration  time  after  the  down-step.  Of  course,  equal-time  integrations 
for  stepped  loading  and  unloading  of  the  system  and  the  statement  ■  0 
at  the  time  of  the  down-step  are  equivalent  to  stating  that  the  system  is 
once  again  in  static  equilibrium  by  the  end  of  the  integration  ^jerformed 
after  imloading.  The  purpose  of  the  foregoing  exercise  was  to  examine 
the  relation  between  system  design-excitation  and  measxirement  errors, 
through  a  breakdown  of  the  over-all  integration  into  separate  integrations. 


-nt 


l^cos  ^  sin 


Aga 


in  utilizing  the  relations  kz  =  F,  =  Jq  Fdt  there  results 
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TRAPEZOIDAL  FORCE^TIME  INPUT 


Tlie  U9«  of  Eq«  77  for  the  more  gener&l  caae  of  trapezoidal  force-time  iiipute 
consists  of  its  successive  use  vitb  each  of  the  four  separate  phases  of  the 
system  excitation,  identified  by  reference  to  Fig,  8.  Phase  numbering  is 
based  on  s>maetry  considerations  rather  than  on  chronological  order.  The 
force-time  input  excitation  consists  of  tvo  ramps,  one  up  and  one  dovn, 
separated  by  an  intermediate  constant-level  force  whose  duration  is  again 
great  enough  to  obtain  essentially  static  equilibrium  at  the  time  of  occiur- 
rence  of  the  dow-romp.  Force-time  inputs  of  this  type  provide  useful 
analytical  approximations  to  rocket  motor  thrust-time  curves. 


Phase  1:  Startup  Ramp 

Initial  conditions  are  ^  ’  0,  and  the  excitation  force  is  Q  «  Kt, 

where  K  is  the  slope  of  the  ramp,  during  the  time  period  0  i  t  ^  tj.  On 
substitution  of  these  expressions  into  tq,  77,  the  displacement  x  of  the 
mass  m  during  phase  1  is  given  by: 


which  results  in 


-nt  W  t  e 
ae _ _  1 

01  m  Ln^'  + 


2nWje"*^  -  (n^  -  -in  Wjt  -  2nCiJ^  co. 


H  - 


t  ^ 

(n  + 


Wj  ) 


(87) 
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On  multiplying  this  equation  by  the  spring  constant,  k,  and  utilizing  the 
relations  given  in  Eq.  80  there  is  obtained  the  expression: 


2Kn(e'*°  cos  t  -  l)  K(n^  -  e~°'’  sin 

F  =  Kt  +  - -5 — - -  +  - r — - 5—2 - - 

n^  +  +  t*^^) 


2\  -nt 


(88) 


The  expression  relating  the  impulses  of  the  applied  and  measured 
startup  forces  is  obtained  by  integration  of  Eq.  88  to  time  that  in. 


^1 

f  Fdt  =  k[  tdt  ^  ^  !.  dt  ^  ^ 


‘1 

-  dt  +  - ~~T - ST- e  sin 


2Kn 


2  ,o2  !0 

n  +  hi 


2Kn 


2  ,,2 

n  +  hi 


f  ^  -nt 

Jo  * 


dt 


(89) 


The  integrated  result,  utilizing  the  expressions  for  measured  and 

Kti2 

Fdt  and  I,  = 


^1 

actual  impulse  during  phase  1,  respectively  =J 
provides  the  error  expression:  ^ 


2Knt, 


^1  ’  'a  - 


2  2  2 


jr -  «  I  .Id  o^t,  , 


(3n^  -  CO*  -  3d^  .  (9*') 

Eq,  90  shows  that  the  impulse  error  incurred  during  the  startup  phase  is 
a  function  of  the  slope  of  the  startup  ramp,  the  damping,  the 
damped  natural  frequency,  and  the  startup  duration  time. 


Pliase  2,  Cutoff  Ramp 

Ktj 

Initial  conditions  arc  *  ^0  *  ^  essentially,  and  the  excitation  force 

is  Q  •  T —  (t  -  t)  during  the  time  period  0  ^  t  ^  t  ,  As  before,  the  state- 
2  2 

ment  x^  »  0  is  considered  entirely  realistic  for  the  contemplated  ayatea  and 
its  excitation.  It  will  be  demonstrated  later  that  no  measureable  error  is 
inctirred  by  stating  «  0  at  this  point  of  the  contemplated  trapezoidal  excit¬ 
ation  of  the  planned  total  impulse  measurement  system. 


Oil  nubatituiioii  of  ihoao  cxprcsBiona  in  Kq,  77,  (he  diaplacomcnt  x  of  th© 
mnsi  m  during  phase  2  is  given  by: 


cos  CO, t  f  ~  sin 


KtjC  sin 


CO^m 


cos  W^t'dt' 


Kfc.e""*'  sin  OLt  -t  ,  , 

e  ain  --  — -  Jq^  W^t'clt'  -*■ 


^irntg 


K  e  cos  to  t  pt  , 

- 1 - - i-  t  c  sin  tO,t’dlt' 

tOj^B  ^2  •*  0  ^ 


(91) 


vhich  results  in: 


Kt  e""*  r  n  1 

[cos  tOjt  .  sin  to^ij .  ■-ia^T“L;;27^ 


Kt.e  n  sin  (Oj^t  f  tO^^  cos  tO^t 


2  ,.2 
n  +  W. 


to  to 
1 


nt 


2nt0^e”*’  “  (n"  -  tO^^)  sin  tO^t  -  OntO^  cos  tO^t 

^  “"2””  2.2  J 


(92) 


On  multiplying  this  cqiintion  by  the  spring  constant  k  and  again  utilizing  th® 
relations  given  in  Fq.  R('  there  is  obtained  the  expression: 


Kt, 


hr 

t„  L 


Kt  + 


2Kn(e 


-nt 


coa  a^t  -  1)  K(n^  (0^^)  e*'"*’  sin 

- - -  .  - - J 


2 

n  + 


n' 


ta^(„"  4 


(93) 


The  expression  relating  the  iiopulses  of  the  applied  and  measured  cutoff 
forces  is  obtained  by  integration  of  Eq.  93  to  time  t^,  i.e., 


The  integrated  result,  utilizing  the  expressions  for  measured  and  actual 
impulse  during  phase  2,  respectively 

r^2 

Im  =  Jo  F  dt  and  , 


provides  the  error  expression: 
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Coapariftrn  of  th«  two  orror  on8  for  phaaos  1  and  2,  Flq.  9'^  and  95 

ghova  thAc  complet<*  cancellation  of  the  two  errors  occurs  for  aymmetric  J 
startup  i.atoff  romps  (i,o,,  tj  ^  t^^).  With  naymmotrical  rarapa,  (i.e,, 
^  incomplete  cancellation,  in  general,  results  in  n  net  impulse 

measurement  error  d\iring  these  two  phaaen  of  the  system  excitation.  Simi¬ 
lar  results  are  runv  (»btained  for  phases  *5  and  4, 

Phase  3  3  f  an  t-ljt*vel  V v 


Initial  conditionH  of  this  phnae  are  obtained  from  the  displacemeri t  ('xpreHsion 
(l’!q.  B7)nn(l  its  derivative,  at  the  instant  of  terminaticn  of  plmne  H .  Tiie 
result  is: 


The  excitation  for(M‘  is  Q  Ktj  during  the  timt*  period  0  ^  i  ^  t^.  On  suh- 
Htitution  of  these  express i orif^  into  Mq.77  >  mul  tipi  J  eut  ion  of  the  riv^ulting 
equation  by  the  spring  cnriMleni,  K,  and  usi*  of  the  expiessions  In  Hq.  BO.  the 
expressiori  for  the  measured  force,  I'\  du*'ing  phase  *5  Is; 


,  ^  ri(., 

t.ie  (’OS  CU,  T  \  I - - 

•  1 


a,,(  ,  tJ 


( ‘>V) 


where  C.  and  (1  are,  respectively,  the  brac'/efed  terms  in  Kq.  9((  for  and 
ijj.  Thp  impulae  of  the  measured  force  F  during  phatt^o  "5  i.i  obtained  by 

the  integration  of  Kq.  fo  time  t,.  Tin*  roHulting  ('»:pr«*Ha  i»>ii ,  utilizing 

'  *  1  «  ' 

the  relations  I..  L,  ^  F  di  and  I.  K  t,  t_  for  tJu*  measuri'd  and  actual 

M  J  0  A  11 

impulses,  respectively,  provides  th©  error  expreasioii  for  phase  3,  ai 
shown  f nl I  owl ng. 


h-h-  Im 

=  - ^ - ^^[n(n^  -  3tu.“)(c  ^sin  +  e""^5  sin  W^t  )  + 

W^(n  +Wj  )“ 

CUj(3n^  -  Ul^^)(i‘  cos  -•-  c  cos 

a;ij(Wj^  -  3n^)  ^  *"3^  cos  «j(tj  ^  t^)  - 

n(3Wj^  -  n^)  0  ^  ^3)gi,i  (tj  *  t^)  -t- 

(CC^^  -  3n^)]  (91 

» 

Equation  98  shows  that  the  impulse  error  incurred  during  the  constant-level 

phase  is  a  function  of  the  system-excitation  parameters  previously  stated 
♦ 

for  phases  1  and  2  plus  the  constant-level  duration  time  t^. 

Phase  4 .  Post— Cutoff  Period 


Initial  conditions  of  this  phase  are  obtained  from  the  displacement 
expression  feq.  92)and  its  derivative,  at  the  instant  of  termination  of 
phase  2.  The  result  is: 


)  -nto  .  - 

=  "IT  ^t2(n2  +  "  - 


It 

Kt,  ,-e-"^2 


cos  ^1^2  1  "I 


e  cos 

.  (  2  ,,2n  J 

t2(n  "  ) 


The  excitation  force  during  the  post-cutoff  period  is,  of  course ,  zero , 

On  substitution  of  these  expressions  into  Fq.77,  multiplication  of  the 
resulting  equation  by  tlio  spring  constant,  k,  and  use  of  the  expressions  in 
Eq,  80,  the  expression  for  the  measured  force,  F,  during  phase  h  is: 


C,  +  nC. 


-  .  'T  IIV^^ 

F  =  Ktj[c  e""''  cos  - i) 


-nt 

e  sin  CC, 


(100) 


51 


where  and  are,  respectively,  the  bracketed  quantities  in  Eq.  99  for 


’■o  *0- 


The  impulse  of  the  measured  force,  F,  during  phase  li  is  obtained  by 

the  integration  of  Eq.  100  to  time  t,  .  The  resulting  expression,  utilizing 

r  ^4  * 

tile  relations  1.,  =  F  dt,  I.  =  0,  for  the  measured  and  actual 

M  J  0  A 

impulses,  respectively,  provides  the  error  expression  for  phase  h: 


-  Ivf 
A  M 


|^n(n^  -  3^^)  (e  sin  ^  sin  Cd^t^^)  + 


i0j(3n‘'  -  (c  cos  ^ 

Cdj(cOj^^  -  3n^)  e  *^^^2  cos  ^(^2 


n(3aj^  -  n  ) 


cos  t^(t2  +  t^) 


tUj  ({d^^  -  3n^)J 


(101) 


The  similarity  of  the  error  expressions  for  phases  1  and  2,  and  for  phases 

3  and  4f suggests  their  joint  consideration,  and  useful  simplification  is 

obtained  by  expressing  the  damping  parameter,  n,  and  the  damped  natural 

frequency  in  terms  of  the  damping  ratio  h  =  — ^ —  and  undamped  natural 

^crit 

frequency  With  reference  to  Eq.  50  and  59>  the  following  expressions 
are  obtained: 


(102) 
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The  use  of  these  relations  in  the  sum  of  Eq.90  and  95  esults  in  the  error 
expression  for  phases  1  and  2: 


h  (5  -  411**)  sin  Cdtj  ^1-h^  + 

1)  cos  CO  tj^  Jl  -  h^  -  l)J  - 

1^-  h  (5  -  4h^)  -  h^  + 

Jl  -  ('ih^  -  1)  (e“"‘'^2  cos  -  l)]}  (lO^) 


Because  the  actual  impulse  for  phases  1  and  2  is  given  by 

I,,  -  =  1/2  Kt  ^  +  1/2  K  t  t  =  ^^(t,  -  t  ) 
Al-2  ^  1  '  1  2  2  M  2’ 


^^1*2 


2  V 


the  fractional  error  for  phases  1  and  2  is 
E  *  E„  2(E  *  E,) 

^  X  ^  J  *. 


1,2 


Al-2 


K  h‘2  Ct: "  0 


and  on  using  the  notations 


1 

2 


Vi  - 

h(5  -  Ah^) 

^1  -  h2 


(104) 


(10-3) 


(106) 
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the  fractional  error  for  phases  1  and  2  is  expressed  by 


-COhti 


- — ^ -  {n_  sin  CClI.t  -  11  cos  WiLt.  -  l) - 

o  ,ti  ,  i  A  3  11 

urtit2(-ji  +  1) 

[n^e  ^^1^2  ~  ^**^2  cos  ^  -  Oj]"  (lO?) 


Similarly,  for  phases  3  and  h 


(-..(3 --.=)( 


-oihti  .  ,,,  f.  72  -a-htT  .  ,  . 

e  sinWtj,yi-h  +  e  -^sinCCt^ 


^1  -  h^  -  l)(c~^*'^l  cos  (jCt^Jl  -  h^  +  e  ^^*5  cos  Ut^Jl  -  h^) 

/l  -  h^  (/ih^  -  1)  c"“  ''^^1  ^  ^3)  cos  W{tj  i^)  J\  -  + 

h(3  -  ^h*")  e  ^  *'^^1  ^3^  sin  oXt^  t_)  ./I  -  h~  -  (^di^  -  l) 


-~  1^-  li  (3  -  /jh“)( 


2w  -Whto 

)  ( O  SI 


sin  W  ,yi  -  h*"  e  sin  CO  -  h'^) 


71  -  ir  (^:r  -  ^  ^  _  ,2) 

,/l  -  h^  {kh^  -  1)  e"^  ^^^2  ^  H)  cos  ^  Ji  _  ,j2  ^ 
h(3  -  4h^)  ^  H)  sin  coCt^  +  t^^)  Jl  -  h^  - 


/i  _  n2  .  „-Wht;. 


yTTIr  (Mr  -  1)  ]  } 


(108) 


5^ 


Because  the  actual  impulse  for  phases  3  h  is  gii'cn  by 


'a3--.  “  ‘iS 


(100) 


tic  fractional  error  for  phases  3  and  k  is 


E  +  E  +  E^^ 


3,4  ■  I 


A3-4 


K  tjt3 


(110) 


rid  on  again  using  the  notations  in  Eq.l06  the  fractional  error  for  phases 
'  and  h  is  expressed  by 

.  =  - -  |h_  sin  CO  H.t.  +  e  ^^^3  sin  COH.t  -  e~^  ^^^1^^3^cos  COHi  (t.  +  t_) 

13 

”3  cos  +  e"“'^^3  cos  H^t^  -  e~^^  ^^l‘"^3)cos  C0H^(tj^+t^)  -  l}- 

t^  si"  ^  sin  W  H^t^  -  e"^  ^^^2'*-i^4)co8  tOH^Ct^+t^))  - 

cos  COHjtg  +  CO  Hj^t^  -  cos  C0Hj^(t2+t^) 

(111) 
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DISCUSSION 


A  discussion  of  the  over-all  measurement  error  is  appropriate  at  this 
point.  7^  begin,  it  should  be  stated  that  a  methematical  integration  of 
an  indicated  force-time  curve  performed  to  an  infinite  time  vill  ghov  that 
no  error  occurs  in  the  measurement  of  total  impulse  of  any  kind  of 
excitation  of  a  spring-mass  system.  The  purpose  of  this  analysis,  however, 
is  to  obtain  the  expression  for  the  error  that  is  incurred  for  finite 
integrations  of  indicated  forces  by  a  spring-mass  system  trapezoidally 
excited.  It  is  therefore  reasonable  to  examine  the  sum  of  Eq.  103  ^nd 
108  (total  error)  to  determine  if  it  does  vanish  for  an  infinite  post¬ 
cutoff  integration  time,  t^.  Accordingly,  for  t^  «  oo, 


^  -W  ht 
‘  e 


3 


sin 


*  t  ) 


' - 9  A  ^  ^ 

(t^  -  t^)  ,/l  -  J-  11^  ^  cos 


h(t.  +  t  )  —  -s  -• 

-^cos  (tj  *  t  )  .Jl  -  y  J 


(112) 


Inspection  of  Eq.  112  for  to tal  measurement  error  shows  that  -•  0  as 


tl  +  t^  —  t^.  The  accuracy  of  Eq.  112  for  the  contemplated  measurement 

system  and  a  very  rapid  transient  excitation  is  demonstrated  by  the  fact 
.  -95 

that  E^  ^  1.3  X  10  ^Mb  sec  for  t^^  =  for  a  175-cps  system,  damped 
10^  critical,  subjected  to  a  startup  ramp  of  10,000  pounds  in  5  milli¬ 


seconds  followed  by  a  constant-level  force  duration  of  2  seconds. 


Therefore ,  for  all  practical  purposes ,  Eq.  103  and  108  regarded  either 
separately,  or  jointly  for  total  error,  quantitatively  show  the  relation¬ 
ship  between  measurement  system-design  parameters  and  total  impulse 
measurement  errors  incurred  by  finite  times  of  integration  for  general 
trapezoidal  excitations .  The  only  1 imi tat ion  imposed  on  the  entire 
mathematical  development  of  these  expressions  is  the  statement  of  static 
equilibrium  at  the  time  of  the  excitation  downramp,  and  the  insignificance 
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of  the  error  resulting  from  this  assumption  for  the  planned  system 
and  its  excitation  vas  demonstrated. 

Because  percentage  errors  are  usually  more  desirable  than  absolute 
errors,  Eq.  107  and  111  arc  of  more  interest  for  measurement  purposes. 
Values  of  error  obtained  by  the  use  of  these  equations,  when  multiplied 
by  100,  vill  be  in  the  desired  form. 

A  few  statements  regarding  the  actual  use  of  such  a  measurement  system 
are  appropriate.  Examination  of  the  error  expressions  shows  that  the 
actual  operation  of  the  system  described  will  be  governed  by  the  level 
of  error  that  can  be  tolerated.  The  percentage  of  error  incurred  in  the 
total  measurement,  as  shown  by  the  sum  of  Eq.  107  and  111,  is  a  function 
of  the  system^s  natural  frequency,  CO,  and  damping  ratio,  h,  the 
excitation  symmetry  ratio,  tj^/tg,  and  the  individual  values  of  t ^ ,  tg, 
t^ ,  and  t^^.  All  of  these  design  and  excitation  parameters,  except  t^ , 
are  beyond  control  of  the  system's  operator,  once  a  measurement  system 
is  constructed  and  put  in  use.  Consequently,  the  duration  of  the  post¬ 
cutoff  integration  time  will  be  the  controlling  operational  factor  in 
holding  measurement  errors  below  a  desired  level. 

Any  decision  regarding  the  amount  of  post-cutoff  integration  time,  t^^ , 
is  consequently  dependent  upon  the  sliape  of  the  input  excitation  curve, 
if  errors  are  to  be  held  to  the  desired  level  with  the  constructed 
system. 
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ANALYTICAL  DYNAMICS  OF  THE  ROCKETDYNE 
THRUST  AND  TOTAL  IMPULSE  MEASUREMENT  SYSTEM 

INTRODUCTION 

In  studies  associated  with  the  design  of  the  thrust  total  impulse  measure 
ment  system  for  use  at  EAFB,  it  was  desirable  to  develop  a  mathematical 
model  which  more  adequately  describes  the  d3mamic  response  characteristics 
of  this  thrust  measurement  system.  The  development  and  analysis  of  such  a 
model  is  described  herein. 

Because  the  system  which  was  investigated  possessed  nonlinear  character¬ 
istics,  it  was  necessary  to  adopt  a  numerical  approach  for  the  solution  of 
the  resulting  equations  of  motion.  However,  if  several  appropriate  simpli- 
fying  assumptions  are  made,  it  is  possible  to  solve  the  equations  in  closed 
form  for  a  special  case  which  is  outlined  in  the  following  analysis. 


MATHEMATICAL  MODEL 

In  systems  of  this  type,  it  is  customary  to  describe  the  mathematical  model 
by  means  of  a  network  of  springs  and  dashpots.  In  particular,  the  thrust 
measurement  system  may  be  represented  by  two  springs  and  dashpots  as  shown 
in  Fig.  9. 


Figure  9  •  Two— Degree-of-Freedom  System 
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In  Ftg,  ,  thr  ronrrpi^  block  and  load  roll  havp  been  replaced  by  equiva-* 
lent  Bpringfl  and  da^hpotH,  The  nomenclature  ia  na  follows: 

M  maa  H  o f  t ho  c  on  c ro  to  block 

r 

k  ^  apririg  conaiant  for  the  block  evaluated  on  the  boHia  of 

ahear  deflection  produced  by  a  thniat 

m  ^  masH  of  entire  motor  and  mount  aRaembly 

k  *  spring  con'itunt  for  load  cell 

b,c  damping  coef f i(‘ ientn  (including  windage  efferta)  of  the 

ce 1 1  und  h  lock,  reaper  live  ly 

F(t)  '  applied  thrust 

the  d i Hplacementfl  of  the  block  and  load  cell,  respectively, 
from  their  equilibrium  poRitiona 

It  ahoiild  be  noted  here  that  M  ,  k  ,  and  k  are  uaaumed  to  be  constants,^ 
whereas  m,  b,  and  c  are  variable  quantities  given  in  the  form 

m  ni(t ) 

b  b(x„) 

c  .  r(ij) 

DIF1<™NT1AL  KQUA'flONS  OK  MOTION  FOU  'IIIF  SYHTI'^M 

hy  liigraiige'H  method,  beginning  with 
potential  energies  of  the  system  as 


-  *i)'’  (in) 


!4«>cJ.inn  of  tluM  d i hoiihm  ion . 


The  equatioriH  of  motion  uro  derived 
the  expresHlons  for  th«  klnotir  and 
given  by 


T 

V 


¥  *  I  *  ¥  “*  2 


Variable  k  is  trow  tod  in  tho  final 


00 


Ths?  reaiHting  forcen  may  bo  aBHumed  in  the  form 


fj  -  c(x^)s^ 

^2  b(xo)*o 

Thus,  in  termn  of  the  quantities  T,  V,  and  f^,  Lagrange's  equations  of 
motion  taay  be  vritten  as 


±_ 

dt 


(f ,) 


av  „ 

IT  +  '^  • 
1 


1,2 


(115) 


where  =  externally  applied  forces.  With  the  aid  of  Eq.  113  and  114, 
the  coupled  equations  of  motion  for  the  system  are 

”■  (mig)  +  k(x2  ~  Xj)  +  ^>(*2)^2  “  “  ^"(^)  (^^^) 

There  are  two  sets  of  initial  conditions  associated  with  the  system  (Eq.  II6) 
to  be  considered  in  the  present  analysis: 


1.  x^(0)  «  0 

Xj(0)  0 

2.  Xj(0)  .  0 
Xj^(O)  r.  Uj 


*2(0)  -  0 

x^(0)  -  0 
■* 

x^,(0)  0 

’<„(<>)  %  (117) 


EVAUJATIQN  QJ'  SPRING  CONSTANT 


Th^  concrete  block  in  assumed  to  act  a»  a  beam  built  in  at  one  end  with  a 
cuicentrated  load  acting  at  the  free  end  nhovn  aa  follows. 


p 


By  definition,  the  spring  constant  for  the  beam  is 

k  -  P/6  (118) 

where  d  ia  the  static  deflection  due  to  a  load  P.  Because  the  cross- 
sectional  dimensions  of  the  block  are  of  the  order  of  magnitude  of  the 
length  of  the  block,  the  deflection  due  to  the  effect  of  shear  must  be 
Included  in  addition  to  that  due  to  bending, 


From  ]^£,  5i  the  following  expression  is  obtained  for 


mi 


[l  h  0.71  =  0.10  (hA)^ 


(119) 


where 


E.  Young's  modulus  for  reinforced  concrete 


t 

I 


length  of  beam 

moment  of  inertin  of  c rosH-sec t ion  about  axis  00 
(oh)(ll/L>)- 


eh^ ! 12 


Thus,  the  equivalent  spring  constant  for  the  concrete  block  is  given  by 


k 

c 


p/6 

-  m  1  a'5  1  < 


^  (lA)E^(h/f  ) 


().7l(h/t)^  O.lo(hA)^ 

e/  [  1  ♦  0.7l(hA)''''  -  0.1()(hA)^ 


(120) 


A  eompu<.f*r  program  huH  b«»eri  written  to  oolv«^  the  coupled  oyRtcMo  (Kq,  llh) 
for  a  wide  variety  of  f\inrtionn  m(t),  and  c(x).  Parameter 

studies  and  impulHe  error  eBtimateR  uaing  thia  prograna  will  be  deorribed 
subsequently. 


SPECIAL  CASE:  ONi]  DEGREE  OF  ITlEtlDOM,  VISCOUS 
DAMPING,  AND  TIME-VARYING  MASS 

If  J,he  OHMuraption  iw  nmde  that  f.he  conrietr  block  io  «  »  "Htiff"  that  itM 
motion  cun  be  neglected,  then  the  syetem  of  equations  (Eq.  k)  reduces  to 
a  Mingle-degree-of-f reedom  system  given  by 


dt 


(mx„)  +  kx,, 

£3 


"  f'(t) 


(121) 


In  addition,  if  the  danping  factor  »  b  »  constant,  and  m(t)  is 

assumed  in  the  form  m  »  at,  t  a  >  0,  the  above  equation  is 

linearlaied  and  closed'^f orm  solutions  are  possible. 


Considering^  therefore,  the  lineari^.ed  equation 


(■  -  at)y  4  (b  -  a)A„  +  kx„  -  -  F(t) 


(122) 


write  B  ■'  m  “  at 
o 


Thus,  Kq.  122  reduces  to 


s  x;;  f  a.x;,  »  a„x„ 

e.  1  h.  «> 


a^F( 8 ) 


(121) 


where  primes  denote  differentiation  with  respect  to  o, 

S.J  r  (a-b)/a 
a,,  k/a*" 


F(«) 


-1/a- 
m  -R 

H  . ) 


Equation  123  oay  be  identified  with  an  inhomogeneous  Bessel's  equation  of 
a  special  kind.  The  homogeneous  solution  (x2^H 

=  [AJ^()3Vi)  +  (124) 

where  and  are  respectively  Bessel  functions  of  first  and  second  kind 
of  order  V. 

V  =  b/a 

R  =  (2/a)^/k 

Applying  the  method  of  variation  of  constancs,  the  general  solution  to 
Eq.  123  is: 

s 

x^(s)  =  s^/2®[AJ^(i9V^)  +  BY^(^V"s)  -  yr  a^JjR^)  J  F(y)y~^'^2^f^ 

m 

I'®  /  ^ 

(5Vy)dy  +  rr  a  Y  (/3Vs)  1  F(y)y  ^  (j3Vy)dy  (l25) 

)  U  <■>  V  J 

m 

o 

In  terms  of  the  time  t,  this  may  he  written  as 


(126) 

The  constants  A  and  B  are  now  to  he  determined  from  the  initirl  conditions 
stated  earlier  in  Eq,  117. 

Case  1 

X2(®)  =  0  =  X2(0) 
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It  is  readily  shown  for  these  conditions  that  A  =  B  =  0,  and  thus  the  solu¬ 
tion  is  given  by 


*2^^)  =  ^  I  (m^-ar)"^'^^^ (t) (j9\/m^-aT)dT 

0 

-  Y^{$^/m^T)  J^(m^-aT)"^/^^’(T)jj^(flV^^T)dT]  (l27) 

Case  2 

*2(0)  =  0  *2^°)  =  "2 

These  conditions  yield 


Solving  for  A  and  B  obtains 

A  .  (,T  Y^_j(3VS;) 

B  =  -  (it  (129) 

Thus,  X2(t)  is  given  by 

+  1Ta^a(ro^-at)^T^^“|^J^(flV™Q~®l-)  J  (iiip-aT)“**^^“F(T)Y^(BVni^-aT)dr 


65 


The  solution  for  this  special  case  should  provide  some  insight  in  inter¬ 
preting  the  results  of  the  numerical  treatment  of  the  general  case 
(Eq.  116). 


MODIFICATION  OF  EQUATIONS  OF  MOTION 
FOR  VARIABLE  SPRING  CHARACTERISTICS 

If  the  load  cell  spring  exhibits  nonlinearities,  the  equations  of  motion 
derived  earlier  must  be  modified  by  assuming  first  that  the  force-displace¬ 
ment  relation  for  the  spring  is  given  in  the  form 

f  =  kx  +  a 


where  k  and  Of  are  constants  to  be  determined  from  information  given  by  the 
lead  cell  manufacturer  or  from  direct  static  calibration. 


V*1 


Thus,  the  potential  energy  stored  in  the  load  cell  spring  now  becomes 

> 

=  J 

«  ‘'o 


f  dx 


With  this  modification  of  the  potential  energy,  the  revised  equations  of 
motion  are: 

Mc’xi  +  c(xj)xj  +  (k^-^k)xj  -  kx^  -  a(x2-Xj)^  =  0 

(mx2)  +  b(x2)x2  +  k(x2-Xj)  +  a(x2-x^)^  =  -  F(t)  (l32) 

The  coupled  equations  (Eq.  132)  represent  a  tvo-degree-of-f reedom  system, 
with  the  specif ied  general  type  of  nonlinear  damping ,  the  load  cell  spring 
of  which  has  the  specified  general  type  of  nonlinear  force-displacement 
characteristic. 
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DESCRIPTION  OF  COMPUTER  PROGRAM 


Introduction 

An  IBM  7094  computer  program  has  been  written  which  numerically  integrates 
the  preceding  equations  of  motion  (Eq,  132),  determines  the  cumulative 
impulse,  and  estimates  the  error  in  the  latter.  The  program  includes  a  main 
program  and  four  subprograms,  all  written  in  FORTRAN.  In  operation,  the  main 
program  simply  calls  the  principal  subprogram,  SOLVDE,  once  for  each  set  of 
data.  A  brief  indication  of  the  operations  performed  within  each  subprogram 
is  given  below. 


SOLVDE 


SOLVDE  is  a  subrout ine-type  subprogram  which  serves  the  dual  purpose  of  an 
executive  routine  and  an  integration  routine.  It  is  executive  in  that  it 
’’reads  in”  the  input  data  and  control  cards,  prints  out  headings  and  results 
according  to  the  control  cards,  and  calls  the  other  subprograms.  In  solving 
the  differential  equations,  it  uses  Runge-Kutta  type  integration  formulas. 

A  Runge-Kutta  technique  was  selected  for  convenience  of  starting  the  solution 
and  because  of  its  natural  stability  property  (such  toc’niques  are  always 
stable  for  a  sufficiently  small  step  size).  Because  various  forms  of  the 
integration  formulas  are  given  in  the  literature  (Ref.  6  and  7),  the  for¬ 
mulas  are  written  as  follows,  in  vector  notation,  with  no  further  discussion. 
The  coupled  equations  of  motion  are  written  as  one  equation: 

X  =  G  (t,  X,  x)  (133) 

where  x,  x,  V,  and  G  =  two-dimensional  vectors.  The  given  initial  condition 
on  the  displacement  is: 

x(t^)  =  x^  {13'i) 
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The  step-wise  integration  is  performed  accordingly  as 


X  ,  =  X  +hx  +-^(A  +  B+C) 
n+1  n  n  6  '  ' 

X  -  =  X  +  (A  +  2B  +  2C  D) 
n+1  no'  ' 


where 


A 

B 

C 

D 


hG  (t  .  X  ,  X  ) 
'  n  n  n'^ 


hG  (t  +  "^  ,  X  +  ^  X  ,  X  +  ^) 
'  n  2  n  2  n  n  2 


h 

2 

h 


^  ^  /  A  1*  M  •  -D  \ 

hG  (t  +  TT  ,  X  +  TT  X  +  7-  A,  X  +  -pr) 
'n  2’n  2n  4’n  2' 

hG  (t  +  h,  X  +  h  X  +  ^  B,  X  +  C) 

'  n  ’  n  n  2  ’  n  ' 


and  h  =  the  step  size. 


(135) 


(136) 


RTS  ID 


RTSID  is  a  function-type  subprogram  which  is  called  by  SOLVDE  to  evaluate 

the  right  side  of  Eq.  129.  It  has  arguments  i,  t,  x, ,  x  ,  x, ,  and  x  ,  and 

"t  h  i.  ^  X  ^ 

hence  evaluates  the  i  component  (i  =  1,2)  of  G  at  the  specified  values  of 

time  and  of  both  components  of  x  and  x.  The  proper  evaluation  requires 

information  regarding?  the  time  variation  of  the  system  parameters  and  input 

excitation.  For  the  system-excitation  combination  under  consideration,  the 

mathematical  expressions  for  time  dependence  were: 


1.  Mass  of  entire  motor  and  mount  assembly: 

m  -  at  for  t  ^  t^ 

»(t) »  » 


at«  for  t  -  t 

J3 


B 


where 


m 

c 

a 

t 


B 


initial  mass 

burning  rate  of  propellant 
cutoff  time 
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2.  Thrust  excitation  (trapezoidal): 


where  p,  t^,  tg,  and  t 


3 


t  <  t 

tj  ^  t  ^  t 


for  tg  <  t  <  t^ 
for  t^  =  t 
inputs  to  the  program. 


If  a  more  general  thrust  function  is  desired,  it  can  be  put  into  th**  pro¬ 
gram  in  tabular  form.  A  modified  version  of  RTSID  has  been  written  which 
evaluates  F(t)  by  linearly  interpolating  this  table.  In  all  other  respects 
the  two  versions  of  RTSID  are  identical.  The  trapezoid  in  the  above  ligure 
used  by  the  modified  version  by  inputting  the  table  given  as  follows: 


t 


0 


0 

P 

P 

o 

0 


where  must  be  at  least  as  larp;e  as  Nh;  N  is  the  number  of  inte^^ration 
steps  of  size  h. 
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Damping 


The  damping  funrtionn  and  c(ij)  nnHor.iated  rp.mnpr  ti  ve  1  y  with  the 

iS  I 

niotor-load  cell  mmhinat  ion  and  the  concrete  hi  ork ,  were  exproMNod  in  the 
following  form; 

b  .  (x  .  )  A  .  tank  (B.  x  .  )  ♦  (I . 

1 ^  w  1  'll'  1 

()  I,'.’)  (ir/) 

Thin  exprewHion  im  versatile  in  that,  either  a  fonntnnt  damping  roe  f  f  i  c  i  en  t , 
c.,  can  be  lined  (A.  o),  or  the  function  can  become  nrhitnirilv  (  Iomc  io  ii 

tttep  function  with  a  dincont  i  nu  i  ty  at  Xj  (). 


MPI/LS 


IMPULS  Ih  a  subroutine  type  of  subprogram  which  computes  the  cumulative 
impulse  to  time  T,  mathematically  exprcHHed  by 


I(T) 


1 


s  j  ( t )  j  t  a 


X 


(i) 


1 

)(i(. 


(I  VI) 


and  which  also  eatis^tes  the  error  incurred,  expressed  by 

Gd 


e(t) 


T 


(k  |x,,(0  =  xj(0|  '  »L«o(0  “  5^,(0!  )  ««<  (no) 


where  T  ■'  and  k  and  fX  are  defined  by  fhe  f orre-d  1  Hplaresent  relaflnri 


kx  +  6?  X 


(1/40) 


It  ahoul'l  be  emphaa i zed  that  the  kind  of  error  ideiiiilfied  Ih  tiie  error 
incurred  only  by  failure  to  (oathemat  icn  1  ly  lnte^rate  to  infinity  to  obtain 
the  raaximuiB  value  of  l(t).  A  twttliematical  integration  Inherently  <!omp»'n- 
aatea  for  the  nonlinear  f  orce-d  i  ap  Ineesent  cliarao  ter  i  h  t  i  c  of  the  aystem.  An 
electronic  integration  of  an  indirated  f  firee-t  ime  plot.,  by  meane  of  only  a 
VO  1  tagp--f:ontrn  1  led  oacillutor  and  a  count<>i',  eaiinot  neitarate  the  error’s 
incurred  by  de|Wirture  from  c.a  1  i  brat  i  tin  linearity  from  the  erroiH  iiieiured 
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hy  prpmiure  ierminaiiou  of  iut  i  (m .  Additional  infomstion  would 

nrcf'rtortr i  ly  hnvr  <  o  hr  nuppliod  with  axi  rlortronir  irit^gi’fttion  of  a 
poHMnuain^  a  nnnlinrar  f  orcr  o  i  r^p  lacorarnf  rlmrari  erintir  >  if  ar-rurat#?  meaa- 
urrmenta  of  Inial  impular  tiro  to  l>r  arhirvril.  C!o?>h  idorwM  r  ronipl  i  ra  t  i  on  of 
aqu  i  pmoxif.  would  hr  involved  witJi  tho  iieieMMary  automatir  elertronir  eompen-- 
Nation.  Hence,  all  poHeihle  effort  wan  made  in  toe  hv*-; f.em  Hpai^n  to  obtain 
II  load  cell  and  vn  I  ta^e -control  1  ed  oaiillator  wHh  Mie  hiftheMt  po.*^aible 
decree  of  linearity,  in  or<lrr  to  avoid  thia  type  of  uyntefn  complex  itY»  To 
obtain  information  ^e^a^(iin^  the  error  aaaociated  with  electronic  inte^ 
gration  of  a  thniat  signal  obtained  from  a  avatem  exhibiting  a  imnllneftr 
cn  1  1  hra  1 1  (HI ,  the  computer  i)ro,9/'ttiii  (Ucludcd  i  in  eg  rat  i  (inn  of  the  expreonion 
for  neveral  vmluea  of  dojmrt.ore  from  a  linear  oharact@ri  ntic,  aper  1  f  i  ra  1  ly : 


I 


& 


«,(t) 


dt 


(141) 


in  whlnh  asveral  valuer  of  &  nro  utilised.  The  prastleal  usgfulnoBB  of 
this  kind  oi  nompuintion  rtsulta  from  a  jud^ugnt  of  th@  Nllowabl#  trror  of 
ihia  typO)  hgforo  n  logd^co'll  rgplnooBent  or  gygttn  siechanioal  r«iiiigiia@nt 
la  roquirod. 


Additionikl  rg^rki  g^rg  ^pproprisio  lo  thg  dlgousilon  of  thg  "norMl"  typg 
of  urror,  obt»in@d  by  tgrviniiion  of  thg  intg^,rntion. 


Thg  intogrgnd  of  Eq.  HB  osoiilnigg  about  egro,  for  t  ^  with  th#  trea 
of  osoh  lobg  SMlIor  In  mgnitud#  than  thg  prtegdini  ont  and  gvontunlly 
tonding  to  itro.  For  ililg  rgeson  I(T)  con  hm  trtmt#d  oo  s  oonvtrggni 
eltgmatlng  ggriui  if  iho  vatugg  of  thg  intggretion  liait,  T,  aro  roquirtd 
to  eorrggpond  to  stro  valuos  of  thg  integrand,  Thg  Msiwus  error  tgret  of 
thg  setHee  1b  the  arte  A(t)  of  the  neit  lobe,  so  iliat  thr  iapulee  grror 
(Eq.  119)  ia  related  a@  follows: 


k(t) 


a(t) 


(142) 


The  output  of  the  IMl^JLS  subprograa  i  a  a  table  of  the  va  1  vies  of  t  imp  T, 
corresponding  to  the  /.ern  valiips  of  the  Integrand ,  and  the  r  or  respond  ing 
valupp  of  i(t),  A(t),  and  HK)  A(T)/I(T).  The  latter  value  le  the  percent 
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of  ^rtor  in  th^  iotnl  impulf^f  asnociaied  with  tprrainfttion  of  t.hp 

integration-  The  irape/oidnl  nile  in  uaed  for  the  integration.  Heitoe,  if 
the  aiep  ai/.e  uaed  ia  too  largi^  the  computed  value  a  of  A(T)  may  not  bo 
monotomical ly  decrenoing  in  nmgnitude.  However,  if  the  integration  in  SOLVDR 
has  been  atahle,  the  roriiputed  valuei^  of  A(t)  will  manifenr  a  general 
decreaaing  trend. 

According  to  previoua  optiima,  determined  by  input  parametera,  SOLVDK  can 
call  the  aiibroutine  CHTfflIT  to  give  cathode  ray  tube  displays  of  the  deflec¬ 
tion  of  the  load  cell  and,  if  doMired^  of  tlin*  deflection  of  ttie  concrete 
block.  Each  display  la  phot ogi'aphod  flo  that  the  results  are  preserved  on 
9-^inch  by  0-inth  photogrnphp.  CRTflSJT  determinefl  the  scaling,  density  of 
the  pointa,  etc,,  according  to  input  parametarm.  The  jioints  are  connected 
with  llnea  uaing  library  tape  subprograms. 

AI»PL1CAT10N  Of  rW:ORY  TO  TW: 

R(K3KETDYNE  HYSTl^ 

Numerous  nombinat  ions  of  design  and  operational  |>aramaiors  have  been  pro*^ 
oeaaed  in  the  computer  program  to  make  a  preliminary  determination  of  the 
accuracy  of  the  measurement  of  total  impulne.  The  various  parameter  values 
used  in  these  determinations  represent  the  best  posnlble  eatitoateai  prior 
to  the  actual  construction  and  testing  of  the  measurement  ayatemi  of  theae 
variables.  The  fiarameter  values  listed  in  (hise  lA  below  wrre  regarded  as 
nominal  values  and  moat  of  the  computer  "runs"  for  case  1  consisted  of 
varying  one  of  these  parameters  while  holding  the  others  at  the  nominal 
values.  To  conserve  computing  time,  t^  was  held  at  two  seconds  for  most 
of  the  runs  because  the  worst  percent  errors  are  associated  with  the  short*- 
eat  excitation  times,  other  parameters  being  held  constant.  Input  excita**> 
lion  times  of  H  and  10  seconds  w*ere  also  computed  with  appropriate  corre» 
spondiiig  values  of  t,,  t,^,  and  propellant  burning  rate  (a).  The  imrameter 
values  end  assoc  i  a  ted  impu  1  se  error  resu  Its  11s  ted  1  n  ('ase  1 H  arc  res  tr  i  r  ted 
to  those  of  linear  (vihcoum)  damping,  linear  load-cell  displacement  chnra- 
ieiisticH,  and  variable  mol.or  mass  during  excitatituu  Case  i.nritniUN 


v&lu@3  aaaociated  with  nonlinear  system  calibrafci  ofib  and  nonlinear  damping, 
together  with  the  aaaocinted  impulsp  meanureroent  errors . 


Case  1:  Linear  (viscous)  damping  and  linear  aystem  calibration 
A,  Nominal  values  of  ))arameterH  (Run  No .  1 ) 


M 

c 

k 

c 


maBH  of  concrete  block ,  3^7  «  1  ugn  (  1 1 ,  *100  pounds  ) 

spring  constant  of  concrete  block,  1,02x10^  Ib/ft 
(50  times  stiffer  than  load  ce  11) 


a 


m 


at. 


P 

t. 


x(0),x(0)  .. 


conaiftni  damping  of  oonorctp  block, 

1,H12  X  10'^  111  (damping  ratio  «*  0,175) 

constant  damping  of  load  cell,  1.2  x  lo’  lb  aec/ft 

spring  constant  of  load  coll,  2.04  x  10'  Ib/ft 

initial  mass  of  motor  and  mount  assembly^ 

15.55  slugs  (500  pounds) 

terminal  mass  of  motor  and  mount  assembly, 

13.35  slugs  (430  pounds) 

amplitude  of  applied  thrust,  10  pounds 

rise  time  of  applied  thrust,  0.005  second 

cutoff  of  applied  thrust,  1.95  seconds 

total  time  of  applied  thrust,  2  seconds 
initial  conditions,  0 


The  results  obtained  with  the  above  nominal  values  of  parameters  are  given 
in  the  acoompanying  plots  (Fig.  10  and  II )  of  load  cell  and  concrete  block 
displacement  and  in  Table  2.  An  ciTor  of  only  0.00117^  is  Incurred  by 
terminating  the  integration  concurrently  with  motor  burnout  (2  seconds). 


B.  Other  values  of  parameters  (all  those  not  listed  were  held 

at  the  nominal  vnlueB):  all  listed  errors  are  those  incurred  by 
tertninating  integration  at  burnout  ('J  Meconds), 

Run  No .  2 

fj 

roiicrfe*i4'  hlnck  damping.  1  20H  x  10  lb  hcc  ^f  (, 

(30^  I  CMH  t.luiii  tuimiitnl) 


V 


Impnl BP 

error  (ai  2“Berondfl  itHrgral.ion 

time);  -0.001 IR^ 

Run  No. 

''i  ^ 

eonrreto  block  damping,  2  ^*10  n 
(lO^  more  than  normal) 

1  0  *  1  h  see /f  t. 

Run  No. 

A 

k 

conrro  f<'  h  1  !»rk  ppr  j  ng  re^UR  Lani , 
(oTily  T  timea  atiffor  lluMi  load 

<1.12  X  10'  Ih/fl. 
rel  I) 

10  1 h  see/ f  < 

"i  " 

concrete  block  fiamping,  h  x 

(damping  ratio  0.  I?*”)) 

Trap\il  H'* 

error  (at  2-BecnndR  int.egra t.ion 

iiitie):  -O.OOnO'K 

Run  No. 

k 

cnnrreie  block  spring  rnnat.ant  ^ 
(200  iimoa  at.iffer  than  load  re 

>4.08  %  10*^  ib/n. 

n)^ 

0,  . 

cone  re  in  block  damping^  1. 1>2^  % 
(damping  ratio  ^  0,175) 

10^  Ih  see/ ft 

Impuln? 

error  (at  g«aecoridB  integration 

time):  -O.OOlORli 

Hun  No . 

load  cell  damping,  1.020  x  10^ 
(damping  ratio  »  O.O^O) 

Ib  ser/ft 

Impulse 

error  (at.  S^aeconda  Integrat.ion 

time);  -O.OOllOi^ 

Him  No, 

J 

‘g  " 

1 

^  oad  re  1 1  damping,  1.171  ^  IR 
(damping  ratio  ^  O.O^O) 

11)  ser/f  t, 

Impuloo 

error  (at  2«aerondM  integration 

time):  -O.OOIK)^ 

Hun  No. 

k  e 

I  oad  cell  apr  i  ng  t'ona  t  an  t  ,  1  ,  H 

s  10^  Ih/rt, 

V  « 

2 

load  (ell  damping,  1.185  x  10^ 
(damping  ratio  «  0.015) 

Ih  HIM'/ ft. 

Tmpiil^p 

error  (at  2-BecondB  int.egrat.lou 

time):  -0.00117“^ 

Run  No. 

0 

p  aeapl  i  t  of  applioil  t^hrusi.^  1000  pounds 


Impulso  rrriM'  (at  2-m<mmiiuIh  i  ti  i.ogrit  <  i  on  limi»): 


0.001 17^ 


ilun  Nfc.  10 

t|  rise  time  of  applied  thrust^  0.02*)  seconds 

tg  "  cutoff  of  applied  thrust,  U.92^  aecondn 

-  total  time  of  applied  thrvist,  ^.000  aer.onde 
a  -  propellant  burning  rate,  sluga/ser 

Impulse  error  (n(,  'j-^econda  integration  time):  0.00031^ 

Run  No.  11 

tj  «•  rlor  time  of  applied  thrust,  O.O'JO  aerorids 
tg  r'»itoff  time  of  applied  thruot,  '1.90  seconds 

t^  total  i.ime  of  applied  tlirust,  10,00  seconds 

a  w  propellant  buniing  rate,  0.2222  BlugB/sec 
Impulse  error  (at  I0»se('on(ls  Integration  tine);  -O.OOOOb^ 


Cftis  2:  Nonlinear  damping  andA>r  nonlinear  syetem  calibration 
Run  No,  12  (nonlinear  load  cell  damping) 

Af,  faaximun  change  in  load  cell  damping,  lb  sec'^fi 

Bg  rapidity  of  change  in  load  cell  damping,  10  sec^ft 

(damping  ratio  is  0.011  'O.OlO) 

Impulse  error  (at  2-seconds  integration  time):,  -0,00111, ^ 

Hun  No,  11  (nonlinear  f orae>Hligplacenent  relation) 
Nonlinearity  coefficient,  “7.1  %  10^^*  Ib'^  cu  ft 

Impulse  error  (at  L’-seronds  integration  time);  -O.OHH^ 

(error  caused  by  using  linear  calibration) 

Run  No,  111  (nonlinear  force-displanenent  equation) 

Nonlinear  coefficient,  “11  «  10^'“  Ib^cu  ft 

Impulse  error  (al,  g=soconds  integrall«>n  lime):  -l.HZ^ 

(error  caused  by  using  linear  roMbraiion) 

llun  No.  11  (both  damping  and  forcc-illsplncomenl  ,  iiotilinosr) 

Ag  and  same  as  in  Run  No,  IE 
Hase  as  Hun  No.  Pi 

Impulse  error  sane  sh  In  l(\in  No.  Pt 
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XriWE  C  SECOMDS'I 


(133d )  lN3W3DVldSia 


TABI£  2 


TOTAL  IMPULSE  MEASUREMENT  DATA 
(CASE  lA:  NOMINAL) 


Integration 

Time,  seconds 

Cumulative 
Impulse,  lb  sec 

Error, 
lb  sec 

Error, 

2.0000 

19725.2 

-0.23 

-0.00117 

2.0025 

19725.0 

0.20 

0.00103 

2.0050 

19725.2 

-0.18 

-0.00090 

2.0075 

19725.0 

0.15 

0.00078 

2.0100 

19725.2 

-0.13 

-0.00067 

2.0125 

19725.0 

0.13 

0.00064 

2.0155 

19725.2 

-0.10 

-0.00051 

2.0180 

19725.1 

0.089 

0.00045 

2.0205 

19725.2 

-0.078 

-0.00039 

2.0250 

19725.1 

0.067 

0.00034 

2.0255 

19725.1 

-0.063 

-0.00032 

2.0285 

19725.1 

0.050 

0.00026 

2.0510 

19725.1 

-0.045 

-0.00023 

2.0555 

19725.1 

0.039 

0.00020 

2.0560 

19725.1 

-0.034 

-0.00017 

2.0585 

19725.1 

0.032 

0.00016 

2.0M5 

19725.1 

-0.025 

-0.00013 

2.04^0 

19725.1 

0.022 

0.00011 

2.0465 

19725.1 

-0.020 

-0.00010 

2.0490 

19725 . 1 

0.017 

0.00009 

ro 

0 

19725.1 

-0.015 

-0.00007 

2.0540 

19725.1 

0.014 

0.00007 

2.0570 

19725.1 

-0.011 

-0.00006 

2.0595 

19725.1 

.0098 

0.00005 

2.0620 

19725.1 

—  .0086 

-0.00004 

2.0645 

19725.1 

.0074 

0.00004 

2 . 0670 

19725.1 

-  .0070 

-0.00004 

2.0700 

19725.1 

.0055 

0.00003 

2.0725 

19725.1 

-  .0049 

-0.00002 

2.0750 

19725.1 

.0043 

0.00002 
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TABLE  2 


(Continued ) 


Integration 

Time,  seconds 

Cumulative 
Impulse,  lb  sec 

Error, 
lb  sec 

Error, 

1^ 

2.077-3 

19725.1 

-.0037 

-0.00002 

2.0800 

19725.1 

.0032 

0.00002 

2.0825 

19725.1 

-.0030 

-0.00002 

2.0855 

19725.1 

.0025 

0.00001 

2.0880 

19725.1 

-.0022 

-0.00001 

2.0905 

19725.1 

.0019 

0.00001 

2.0930 

1S725.1 

-.0016 

-0.00001 

2.0955 

19725.1 

.0015 

0.00001 

2.0985 

19725.1 

-.0012 

-0.00001 

2.1010 

19725.1 

.0011 

0.00001 

2.1035 

19725.1 

-.0095 

-0.00000 

2.1060 

19725.1 

-.0082 

0.00000 

2.1085 

19725.1 

-.0077 

-0.00000 
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THEORY  AND  APPLICATION  OF 
PARALLELOGRAM  COMPRESSION-TYPE  FLEXURE 
SUSPENSION  SYSTEMS 


INTRODUCTION 

Analysis  of  compression-type  flexure  suspension  systems  vas  made  primarily 
for  design  clarification  purposes  so  that  a  suitable  suspension  system 
could  be  provided  for  the  BATES  motor  in  the  Rocketdyne-designed  test 
stand.  The  need  for  such  analysis  vas  established  by  the  dual  design 
objectives  of  high  system  dynamic  response  to  transient  thrust  forces, 
and  a  high  degree  of  linearity  in  the  measurement  system  indications  of 
static  forces.  These  two  objectives  resulted  in  the  proposed  use  of  a 
load  cell  having  no  compliant  in-line  flexures  or  ball  joints  with  a 
motor  suspension  system  whose  operational  deflections  resulted  in  suit¬ 
ably  small  load  cell  bending  stresses.  The  objectives  were  accomplished 
through  the  use  of  a  flexure  system  whose  static  translational  and  ro¬ 
tational  stiffnesses  were  suitably  related  to  those  of  the  load  cell. 

The  design  criterion  vas  based  on  the  load  cell  electrical  response  to 
the  bending  moments  imposed  on  it  by  mechanical  system  distortions  re¬ 
sulting  from  motor  thrusts  and  weight  changes,  calibration  forces,  and 
thermal  expansions.  Accordingly,  mathematical  expressions  of  the  flexure 
system  translational  and  rotational  stiffnesses  as  functions  of  flexure 
component  and  system  dimensions  were  derived.  The  complementary  stiff¬ 
nesses  of  the  proposed  load  cell  and  i ts  associated  electrical  responses 
were  obtained  from  its  designer,  the  Baldw’in-Lima-IIamil ton  Corporation. 

The  load  cell-flexure  system  combination  vas  then  examined  to  determine 
the  flexure  system  dimensions  which  produced  a  force  and  moment  distri1}u- 
tion  yielding  the  nece*-  •.  v  low'  level  of  load  cell  bending  moment  associ¬ 
ated  with  the  needed  low  level  of  measurement  error.  The  results  of 
this  process  consequently  yielded  the  proposed  flexure  system  design, 
and  the  steps  in  this  process  will  now  be  discussed. 
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DERIVATION  OF  FLEXURE  S^'STEM  STIFFNESS  EQUATIONS 


With  reference  to  the  figure  belov  and  Fig.  12,  the  following  equation 
in  which  E  =  Young’s  modulus  and  6  =  deflection,  have  been  developed 
for  a  flexure-strut  combination. 


Vertical  Stiffness 

-  Ft 


Ver ti cal  deflection: 


table 
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Vertical  Deflection 


Because  two  flexure-strut  combinations  are  incorporated  in  the  parallelo¬ 
gram  suspension  system,  the  vertical  stiffness  of  the  flexure  system  is: 


(145) 
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In  the  case  of  a  composite  beam  constructed  of  two  relatively  short, 
rather  thick,  flexure  sections  separated  by  a  long,  thicker  strut,  all 
junctions  of  the  sections  are  builtin  and  the  displacement  of  the  movable 
support  is  considered  to  be  translation  only  (no  rotation).  The  curvature 
of  the  entire  flexure-strut  combination  is  therefore  symmetrical  about 
a  line  through  the  center  inflection  point  C,  parallel  to  the  planes  of 
the  tvo  supports.  Consequently,  the  displacement  of  the  movable  support, 
associated  vith  the  bending  force  F„,  is  twice  that  of  the  center  point 

o 

C.  Considering  only  bending  of  the  flexures,  the  vertical  reaction  at 

the  fixed  support  is  R,  =  F„  and  the  couple  exerted  by  the  wall  is  M,  = 

IB  t 

-Fg  (L  +  The  bending  moment  at  any  cross  section  mm  of  the  flexure 

adjacent  to  the  fixed  end  is  M  =  Mj  +  Rx  .  Using  the  well-known  differ¬ 
ential  equation  of  the  deflection  curve  of  a  bent  beam, 

El.  ^  =M  (li*6) 

^  dx 


there 


is  obtained, 


El  ^ 
dx® 


by  substitution  for  M: 

=  -Fg  (1  -  2<.j)  *  FgX 


(l'>7) 


where  1^  is  the  moment  of  inertia  of  the  flexure  section  about  a  horizontal 
transverse  axis  through  its  section.  A  first  integration  gives  the  slope 
equation  for  the  flexure  section: 


dy 

dx 


-Fg  (L  .  2<.,)  ,  .Fg  , 


(US) 


and  because  dy/dx  =  o  at  x  =  o,  it  follows  that  Cj  =  o.  A  second  inte¬ 
gration  gives  the  bending  deflection  equation  for  the  flexure  section: 

EIj  y  =  -Fg  (L  .  2tj)  I?  .  Fg  .  Cg  (U9) 

and  because  y  =  o  at  x  =  o,  it  follows  that  Cg  =  o.  The  slope  and  de¬ 
flection  of  the  bent  flexure  section  at  the  junction  point  P  at  the  end 
of  the  strut  are  found  from  the  above  equations  by  substituting 

A  • 

Because  for  small  angles  tan  ?  «  9,  it  follows  that 
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slope 
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deflection 


6 

P 


6 

P 


I 

6  +  6" 

P  P 


-  F, 


B 


( 


"  6 


It’ 


El 


1 


s 
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The  deflection  of  the  strut  in  bending  is  obtained  in  similar  manner. 
The  vertical  reaction  at  the  end  of  the  strut  (point  P)  is  =  F„  and 
the  couple  exerted  by  the  flexure  is  =  -F  (L  +  'f-, ) .  The  bending 
moment  at  any  cross  section  nn  to  the  left  of  the  center  point  C  is 
M  =  Mg  +  RgX,  vhere  x  is  nov  the  distance  of  section  nn  from  the  end 
point  P.  Again,  using  the  differential  equation 


=  M 


dx 


(153) 


there  is  obtained,  by  substitution  for  M: 

dx 


(154) 


where  Ig  =  the  moment  of  inertia  of  the  strut  section  about  a  horizontal 
transverse  axis  through  its  centroid.  A  first  integration  gives  the 
slope  equation  for  the  Sv-rut  section; 


El 


2 


dx 


-Fb  (l  .  *  Fj  +  Cj 


and  because 


dx 


e 

p 


at  X  =  0,  it  follows 


that 

k  Elg  I 

AE  J 

s 


(155) 
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and  therefore: 


El, 


^2  dx  ^  V 

r  Ip  -,  2  '' 

[  17  '^'1  ■"  2  S  )  ^ 


X  +  F  - 

^  'b  2 


AE 


(156) 


A  second  integration  gives  the  deflection  equation  for  the  strut  section: 


EIj  y  .  -Fg  (L  -  S)  -  Fj  ^  - 


r  ^2 


"bl  17  S  *  i 


k  EI2 

AE  J 


X  +  C, 


and  because  y  =  6  at  x  =  o,  it  follows  that 
P 


<^2  =  [r 


+ 


k  Elg  ■f'j 

“ae 


(157) 


and  therefore: 


El^y  =  .F„  (Ut,)  ^  F„  -  F„  ^  (1,4^  *  |  x 


-''b  * 


1 


k  Elg  -, 

AE  J 


(158) 


Because  the  strut  thickness- to-length  ratio  is  small  (<  l/20),  deflection 

due  to  shear  is  negligible  and  the  deflection  of  the  bent  flexure-strut 

combination  at  the  center  point  C  is  found  from  the  above  equation  by 

substituting  x  =  L.  The  result  is: 

2 


r  5  ,3  '■%  '2  ,  7, ,2  5  ,3,  '"^'2-^1  1 

.F„  lk8^  ~  ^  TT  ^  “2“  4  ^^1  ^  6  ^1  ^  ~AE  J 

B  1  s 


El, 


(159) 
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The  total  deflection,  6  of  the  flexure-strut  combination  in  the  horiz- 

B 


ontal  longitudinal  direction,  is,  by  symmetry  considerations,  twice  that 
of  point  C,  and  the  stiffness  of  the  combination  is,  numerically: 


k* 

L 


"b 

% 


El, 


24^'^  li  2  3  ^  AE 

1  s 


(160) 


Because  two  flexure-strut  combinations  are  incorporated  in  the  paral¬ 
lelogram  suspension  system,  the  horizontal  longitudinal  stiffness  of  the 
flexure  system  is  given  as  follows,  w^here 


I,  = 


wt" 


WT^ 


1  12  *  ^2  12  ^  ^  wt ,  and  k  2 


have  been  substituted; 


k,  = 


2  EWT' 


5,3  ,,2  .  r  2  .  7  ,  »  2 

2  E  +  3L  "Ej  +  12  — ^  I^L  "Ej  +  g  L  •6^ 


3  EWT^  I 


1 


3  1  J  E  wt 

s 
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Horizontal  Transverse  Stiffness 


-  X 
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STRUT 


In  the  case  of  a  beam  which  has  a  large  thickness  at  all  parts  compared 
to  its  length,  so  that  the  deflection  is  due  both  to  bending  and  to 
shearing,  the  composite  beam  consists  of  relatively  short,  deep  flexures 
separated  by  a  long,  deep  strut.  All  junctions  are  builtin  and  the  dis¬ 
placement  of  the  movable  support  is  again  considered  to  be  one  of  trans¬ 
lation  (no  rotation).  The  curvature  of  the  whole  flexure-strut  combina¬ 
tion  is  therefore  symmetrical  about  a  line  through  the  center  inflection 
point  M,  parallel  to  the  planes  of  the  two  supports.  The  deflection  of 
the  flexure-strut  combination  is  due  in  part  to  bending,  6’’  ,  and  in  part 
to  shear,  6^  ,  of  both  the  flexures  and  the  strut.  The  deflection  due 
to  bending  can  be  stated  by  using  the  expressions  already  obtained, 
with  regard  to  the  proper  expressions  for  the  moments  of  inertia  1^^  and 
Ig  which  are  now 

I  =  i 

^12 

T  (162) 

2  12  * 


and  the  expression  for  the  total  deflection  of  a  flexure-strut  combina¬ 
tion  due  to  bending  alone; 


The  deflection  of  the  flexure-strut  combination  due  to  shear  is  obtained 

bv  the  use  of  the  express i on  dv, /dx  =  kF  /AE  with  both  the  flexures 

^  '1^  s^  s 

and  the  strut  where; 


dx 

k 


slope  of  the  shear  deflection  curve  of  the  section  considered, 

numerical  factor  by  which  the  average  shearing  stress  must 
be  multiplied  in  order  to  obtain  the  shearing  stress  at  the 
centroid  of  the  cross  section,  (k  =  3/2  for  a  rectangular 
section) 
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=  shearing  force 

A  •=  cross-sectional  area  of  the  section  considered 

E  =  shear  modulus  of  the  material 

s 


For  a  single  flexure, 


^  dv,  V  ®  1  o  F 

(  1  ^  _  si  _  2  s 

V  dx  yi  vt  E 

1  s 


and  for  the  strut, 

r  ^ 

V  dx  ^9  L 


WTE 


(164) 


(16-) 


Therefore,  for  the  tvo-f lexure-strut  combination,  the  total  deflection 
due  to  shear  is: 


6" 

9 


wt  E 


6" 
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3F. 


2  E 
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—  F 
2  5 


WTE 


L 
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Thus,  the  total  deflection  of  the  combination,  due  to  both  bending  and 
shear ,  is : 
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Because  two  flexure-strut  coubinations  are  incorporated  in  the  parallelo¬ 
gram  suspension  system,  the  horizontal  transverse  stiffness  of  the  system 
is: 


- 


2F 

s 

s 


2ETV' 


2^ 


12^  (A,  . 


3  1  ’ 


3ETW^ 


2E 


L  ^ 
Vvt  '*■  WT^ 


(168) 


Rotational  Stiffness  Aoout  Horizontal 
Transverse  Axis 


Ro ta tional  moments  about  a  horizontal  transverse  axis  displace  the 
flexure  system  in  tension  and  compression,  and  because  tvo  flexure-strut 
combinations  are  equally  stressed  in  opposite  directions,  the  suspension 
system  stiffness  expression  is: 


Rota tional  Stiffness  About  Vertical  Axis 

A  transverse  couple  about  a  vertical  axis  displaces  the  fl exure  sys tem 
in  shear -bend  i  n^iC,  ®nd  because  tvo  fl  exure- s  trut  combinations  are  equally 
stressed  in  opposite  directions,  the  suspension  system  stiffness 
expression  is: 


r 

v 


F  D 
s 


F  D 


2fi, 


(170) 
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T 

V 


E  Vvt  ^  VT^ 


(170  Conti nijod) 


Flotational  Stiffness  About  Horizontal 
Longi t  nd inal  Axis 


This  stiffness  expression  wa?  not  developed  because  moments  about  (his 
axis  are  not  expected  to  be  encountered  in  apprevi?ble  amounts,  and 
also  because  the  loe*f  cell  response  to  torsio  al  moments  is  nepli^jible. 

The  flexure  sys  tem  stiffnesses,  expressed  by  the  ; reced in^r  equations, 
are  related  to  the  corresponding  load  cell  stiffnesses  by  the  design 
requirements  based  on  measurement  erro^^s .  T!ie  complementary  B-L-H 
proposed  load  cell  stiffnesses  are  given  belov: 


1. 

2. 

3. 


Axial  stiffness:  = 

1  •  b  / 

X  10^  lb  in. 

Transverse 

s  t  i  f  f nesses : 

k  ’ 

V 

=  k^'  =  10^  lb  in. 

Ro  tationa 1 

s  ti f  fnesses  ; 

^r’ 

^  r  ’  =  T  '  =  2")  X  10'’  lb  in.  Vad 
I.  V 

The  electrical  responses  of  the  locid  cell  to  transverse  loads  and 
bending  moments  are  as  follows; 

1.  A  700-pound  transverse  loa<i,  or  a  moment  of  2700  Ib^.  in.  ,  causi^s 
a  maximum  sensitivity  change  of  0.17^^  of  full  scale. 

2.  Cell  responses  to  torsional  {twisting)  moments  are  negligible. 


DETERMINATION  OF  FimTlE  SVSTIM  STIFIM-^ 

REQUIRFMENTS 

The  design  of  the  flexure  system  is  largely  determined  by  the  anticipated 
degree  of  departure  from  a  perfect  thrust  calibration  and  measurement 
system.  An  ideal  system  is  one  in  which  no -error  calibration  forces  are 
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impoflpd  oxactly  on  tho  axis  of  a  noncompliant  force^meaflurinj?  transducer 
which  responds  onlv  to  axial  forces,  connected  to  a  suspension  system  in 
vhich  the  motor  thrust  axis  is  perfectly  aligned  with  the  ideal  transducer. 
With  such  a  system  no  errors  would  result  from  misaligned  calibration 
links  or  motor  thrust  axis,  and  the  load-measuring  tran-^ducer  would 
identically  indicate  equal  calibration  and  thrust  forces.  The  linearity 
of  the  transducer  calibration  would  be  solely  dependent  upon  the  trans¬ 
ducer  linearity  itself.  However,  actual  calibration  and  measurement 
sy s  terns  employed  for  rocket  motor  thrus  t  de termina  tions  mus  t  necessarily 
depart  somewhat  from  the  ideal  system  described,  and  the  purpose  of 
this  discussion  is  to  identify  the  nature  of  the  departures  and  quanti¬ 
tatively  describe  their  effects  upon  the  calibration  and  measurement 
system  accuracies. 


Identification  of  Causes  of  Bending  Moments 

It  should  be  stated,  first,  that  an  accurate  emplo\'ment  of  a  load  cell 
is  obtained  by  merely  imposing  identical  load  conditions  upon  it  during 
calibrations  and  measurements.  With  only  an  accuracy  requirement, 
the  presence  or  absence  of  transverse  forces,  along  with  couple  bending 
and/or  twisting  moments  on  the  load  cell  is  not  material,  provided  that 
the  transducer  loading  conditions  are  duplicated  in  the  calibration 
and  operation  of  the  mencurement  system.  However,  tlie  calibrations 
of  such  systems  are,  in  general,  nonlinear.  Accuracy  plus  linearity 
in  the  calibration  and  use  of  a  sys  t  em  requires  cons iderably  more 
attention  to  system  design  detail.  All  direct-reading  electromechanical 
transducers  require  axial  corap 1 iance  in  some  degree  for  their  opera t ion , 
and  they  are  sens i ti ve  in  some  degree  to  transverse  loads  and  couple 
benciing  moments.  The  cause  of  such  loads  and  couple  moments  is  not 
only  fabrication  and  ass emb ly  tolerances,  but  also  operational  variables 
such  as  shifting  motor  thrust  axes,  decreasing  motor  weight  (burning 
solid  propellant)  and  differential  thermal  expansions.  Because  these 
loads  and  moments  result  in  calibration  nonl inear i tics  and  measurement 
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errors,  it  is  tiseful  to  list  and  quantitatively  discuss  tlieir  sotirces 
as  fol lows: 

1.  I^endin^^  moment  produced  by  flexure  system  axial  stiffness 

2 .  Tons  t  rue  t i ona I ly  misa 1 i  ^lod  calibration  and/ or  motor  thrust 

forces 

a.  Horizontally  and  or  vertically  offset  forces  parallel  to 
load  cell  axis 

b ,  f^orces  angularly  inclined  to  load  cell  axis  in  the  horizon¬ 
tal  and/or  vertical  planes 

3.  Motor  test  and  environmental  conditions 

a.  Decreasing  motor  weight  (burning  propellant) 

b.  Shifting  motor  thrust  axis  (unusual  propellant  combustion) 

c.  Vertical  and/or  horizontal  differential  thermal  expansions 

Ten  or  more  sources  of  load  cell  bending  moment  liave  been  identified, 
and  it  is  useful  to  regard  them  as  equally  probable.  As  already  stated, 
2^00  lb  in.  of  bending  moment  on  tlie  load  cell  produces  a  possible 

change  in  load  cell  sensitivity.  Load  cell  sensitivity  changes 
are  directly  related  to  calibration  and  measurement  nonlinearities  and, 
therefore,  measurement  errors.  Because  it  is  desirable  to  limit  the 
to  f,a  1  possible  error ,  due  to  load  cell  non  1  inear  i  t  i  es  ,  to  values  of 
O.Ol'?  or  less,  it  is  appropriate  to  require  that  the  maximum  transverse 
load  and  the  maximum  amount  of  bending  moment  necessarily  supplied  by 
tlie  load  cell  with  each  single  source  be  limited  to  10  pounds  and 
'lO  lb  in.  respectively.  In  the  following  discussion,  the  individual 
sources  are  examined  and  the  flexure  system  stiffness  requirements 
are  determined  on  the  basis  of  these  limitations. 


Flexure  System  S ti f fness  Requi remen ts 

Flexure  system  stiffness  requirements  will  be  determined  on  the  basis  of 
errors  produced  by  the  most  unfavorable  circumstances  of  system 
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construction  and  operation.  Considering  first  the  forces  ana  moments 
which  can  occur  in  a  vertical  plane  containing  the  load  cell  and  motor 

axis,  in  the  figure  below,  the  concrete  abutment  is  considered  fixed,  and  the 
initial  motor  weight  is  simply  supported  on  the  flexure  system  beneath 
with  no  associated  restraining  forces  or  moments  supplied  by  the  load 
cell.  If  a  weight  W’  is  subtracted  from  the  original  weight,  the  changes 
in  flexure  force  are  Fj  and  Fg  respectively  and  the  associated  transverse 
load  and  pure  moment  supplied  by  the  load  cell  are  R  and  M,  respectively. 


At  equilibrium, 

W*  =  Fj  +  Fg  +  R 

RS  MFj  -  Fg)  D/2  ^  M  (171) 

A  70-pound  vcight.  change  due  to  burned  propellant  should  be  considered, 
letting  S  =  2^  inches  and  D  =  "50  inches.  Because  both  a  transverse 
load  and  pure  moment  are  supplied  by  the  load  cell  with  this  single 
source,  the  maximum  permissible  values  of  R  and  M  should  be  10  pounds 
and  50  lb  in.,  respectively. 

Then: 


70  =  Fj  ^  Fg  +  10 

240  =  (Fj  -  Fg)  D/2  +  50  (l72) 
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Therefore,  ^2  *  flexure  system  vertical  stiffness  is 

seen  to  be  about  six  times  the  load  cell  transverse  stiffness.  In  the 
same  manner,  (Fj  -  Fg)  D/2  =  190  and  the  flexure  system  rotational 
stiffness  about  a  horizontal  transverse  axis  is  seen  to  be  about  four 
times  the  load  cell  bending  stiffness. 

In  the  next  case,  (belov),  a  perfectly  aligned  10,000-pound  calibrating 
force  compressing  the  load  cell  and  deflecting  the  motor  flexure  sus¬ 
pension  system  in  the  thrust  direction  is  considered.  Denoting  by  f  the 
backvard  force  exerted  by  the  flexure  system,  it  is  seen  that  fh  =  M  =  100 
and  for  h  =  15  inches  and  M  =  100  lb  in.,  f  *  6  pounds.  Therefore,  the 
axial  stiffness  of  the  load  cell  should  be  at  least  1700  times  the 
flexure  system  horizontal  longitudinal  stiffness. 


Next,  a  motor  thrust  force  of  10,000  potiiids  vertically  displaced  a 
distance  d  from  the  load  cell  axis  is  roTisidered  (  belov  ): 


MOTOR 


At  equilibrium: 


Fj  -  Fg  -  R  (173) 

and  the  moment  equation  about  eg  is: 

10,000  d  =  (Fj  +  Fg)  D/2  +  M  +  RS  ( I?'*) 

Assuming  that  d  =  l/4  inch  (abnormally  large  misalignment)  and 
S  =  24  inches, 

2500  =  (Fj  +  Fg)  D/2  +  M  +  24  R  (l75) 

So  that  the  transverse  force  R  and  the  bending  moment  M  supplied  by 
the  load  cell  are  10  pounds  and  50  lb  in.,  respectiv’ely,  the  flexure 
system  rotational  stiffness  about  a  horizontal  transverse  axis  should 
be  4^  times  stiffer  than  the  load  cell  bending  stiffness,  and  the 
flexure  system  vertical  stiffness  should  be  15  times  the  load  cell 
transverse  stiffness.  By  consideration  of  a  l/4-inch,  horizontally 
displaced,  10,000-pound-thrust  force,  the  same  conclusions  can  be 
reached  regarding  the  flexure  system  rotational  stiffness  about  a 
vertical  axis  and  the  transverse  longitudinal  stiffness. 

On  the  basis  of  the  foregoing  considerations  of  possible  misalignments 
that  can  occur  in  motor  test  operations,  the  following  design  criteria 
are  given  for  the  flexure  system  stiffnesses: 

1.  Vertical  stiffness:  about  15  times  the  load  cell  transverse 
stiffness  (same  for  horizontal  transverse  stiffness) 

2.  Horizontal  longitudinal  stiffness:  about  I/17OO  of  the 
load  cell  axial  stiffness 

3.  Rotational  stiffness  about  a  horizontal  transverse  axis: 
about  45  times  the  load  cell  bending  stiffness  (same  for 
rotational  stiffness  about  a  vertical  axis  for  transverse 

couples) 
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APPLICATION  OF  THEORY 


Flexure  System  Design 

The  application  of  the  foregoing  flexure  sysfcm  design  criteria  vas 
made  to  the  Rocketdyne  design  of  a  test  stand  using  the  BATT5  motor 
for  propellant  studies.  The  practical  uso  of  the  flexure  stiffness 
expressions  sugges ts  that  sui table  al lovances  be  made  for  fabrication 
equipment  limitations.  Accordingly,  it  vas  reasonable  to  inquire  if 
flexure  dimensions  that  were  reasonable  from  geometrical  and  machining 
standpoints  were  adequate  to  provide  the  desired  flexure  system  stiff¬ 
nesses.  Therefore,  a  determination  of  flexure  stiffnesses  vas  made 
on  the  basis  of  the  following  dimensions; 

1.  Flexure  component  dimensions,  inches: 

=  1/2 
V  =  16 
t  =  1/8 

2.  Strut  dimotisions  and  spacinc,  inrhos: 

L  =  n 

W  =  16 
T  =  l/h 
P  =  30 

3.  Material;  steel,  E  =  30  x  10^’  pp  i  and  E^  =  12  x  10^’  ps  i 

The  flexure  system  stiffnesses,  alrmg  with  t!ie  ratios  of  flexure  s  vs  tern 
to  load  cell  stiffness  that  were  obtained  are  listed  below: 
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Vertical  Stiffness. 


k 

V 


2E 

wt  VT 


13.3  X  10^ 


11/ in. 


k 

Ratio:  =  13.3 


(176) 


Horizontal  Longitudinal  Stiffness. 


k,  = 


2EWT" 


|l’  +  3L^^-j  +  12~  + 


E  vt 
8 


k  =  1000  Ib/in. 
L 


Ratio: 


1 

1670 


(177) 


Rotational  Stiffness  About  a  Horizontal  Transverse  Axis. 


= 


2^1  I 

WT^ 


=  3  X  10”  lb  in. /rad 


Ratio:  =  120 

T 


(178) 
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Rotational  Stiffness  About  a  Vertical  Axis. 
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Ratio:  -  25.6 

V 
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Horizontal  Transverse  Stiffness. 


^^1  = 


2ETV 


.3 


2^ 
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2  3 


3  -t 
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J\  /  I  \ 

1  '  2E  ^vt  Vt' 


=  3  X  lo”  Ib/in. 


_ 

Ratio:  — ,  4  3 


(  180) 


Although  this  last  ratio  is  only  5,  compared  to  the  requirement  of  15,  it 
is  considered  adequate  in  view  (►f  the  magnitude  (  1^ inch)  of  the 
assumed  thrust  misalignment. 

Because  the  obtained  ratios  of  flexure  system  stiffness  to  load  cell 
stiffness  compare  very  favorably  with  the  desired  ratios,  determined 
on  the  basis  of  measurement  errors,  the  stated  flexure  dimensions 
are  used  in  the  suspension  system  design. 


Flexure  Bending  Stress  and  Column  Strength 

Because  the  flexure  system  supports  the  motor  and  motor  mount  assembly, 
an  examination  of  the  design  vas  made  with  respect  to  flexure  bending 
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stresses  and  column  strength.  Considering  the  flexure  system  diagram  given 
in  the  figure  belov^  the  maximum  bending  stress  of  the  flexure  occurs  at 
the  junction  point,  A,  with  the  strut.  Its  value  is  obtained  ns  follows: 


P 


(181) 

where 

M  ■  bending  moment  at  point  A 
I  =  flexure  moment  of  inertia 
6  =  angle  of  deflection  of  the 


system  =  y/h 


P  =  axial  load  on  flexure 
=  flexure  length 


but 

(182) 

where 

S  =  bending  stress  on  flexure  surface 
t  =  flexure  thickness 


M 


SI 

1/2 


Therefore,  by  equating  the  two  expressions  for  M  and  solving  for  the 
stress  S,  there  is  obtained: 


(183) 


and  for  y  =  0.006  inch,  h  =  14  inches,  P  =  275  pounds,  E  =  30  x  10^  psif 
w  =  16  inches,  t  =  l/s  inch,  and  =  l/2  inch,  the  result  is; 


S  «  15  pai 


(184) 
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With  regard  to  coliuon  buckling,  the  expressions  for  the  critical  loads, 
which  cause  buckling  of  long  columns  are  given  as  follows  for 
various  end  attachment  condi tions : 


1. 

2. 

3. 

4. 


where: 


E  =  Young's  modulus  of  elasticity 

I  =  moment  of  inertia  of  cross  section  about  a  transverse  axis 
=  column  length 

The  application  of  these  expressions  in  discussions  of  flexure  and 
strut  design  requires  a  judgment  not  only  of  the  applicable  expression, 
but  also  of  its  limit  of  applicability. 


Single  Flexure.  Considerations  of  the  end  conditions  of  an  elemental 
-  - 

flexure  lead  to  the  use  of  the  expression  P  =  ^ —  in  order  to  remain 

^  (21,)^ 

on  the  conservative  side  of  critical  loading.  The  limit  of  applicability 

in  calculations  can  be  established  by  dividing  the  expression  by  the 

cross-sectional  area,  A,  and  letting  r  =  Then  because  the 

critical  stress  is  S  =  P  /A,  there  is  obtained 

cr  cr 


S 

cr 


(185) 
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renains  wiLhin 


This  equation  is  applicable  as  lonj?  as  the  stress 
the  proportional  limit.  For  structural  steel  with  a  proportional  limit 
of  30,000  psi  and  E  =  30  x  10^  psi i 


(186) 


Consequently,  the  critical  load  for  a  steel  flexure  with 

> 

and  one  end  free,  having  —  =  50,  is  calculated  from  = 


one  end 
fl^EI 


f  ixed 


Nov,  because 


(187) 


For  tlic  flexure  dimensions  =  l/2  inch  and  t  =  l/8  inch,  the  slenderness 
ratio  is: 


_i  -  1-7 

r  TJs 

—  =  13.6  (188) 

s 

Because  —  ^  50,  the  flexure  will  not  fail  in  buckling. 


Strut.  Consideration  of 
elemental  flexures  leads 


the 

to 


end  condi tions 
the  use  of  the 


of  a  strut  separating 

n 

expression  P  =  — tt- 
^  c 


two 

to 
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remain  on  the  conservative  side.  Thus,  for  a  proportional  limit  of 
30,000  pai  for  steel,  it  follows  that: 

^  =  100  ( 189) 


Consequently,  the  critical  load  for  a  steel  strut  with  both  ends  pinned, 

hav'ing  —  ^  100,  is  calculated  from  P  =  — ^  -  , 
r  c  - 


For  the  strut  dimensions  L  =  13  inches  and  T  =  1 ;  inch,  the  slenderness 
ratio  is: 


7  =  177  (190) 

Hecnuse  —  >  100,  the  strut  will  fail  in  bucklinc  before  it  fails  in 
r 

compression.  However,  for  a  strut  with  the  stated  dimensions, 

-  37,000  pounds  or  133  times  the  normal  veiclit  supported  by  each 
strut. 


Pi  f  f  eren  t.i a  1  Thermal  Fxpans  i  on 


Tt  is  useful  to  investigate  the  ]>f>ssible  vi»rtical  differential  expansion 
of  the  concret  e  abutment  and  the  s  i  eol  .;i  1  uminiim  motor  mount  and  sus¬ 
pension  system.  With  reference  to  the  tiolov,  tlie  Vf'rtical  e\j)ansiiin 

are  ca 1 cul a  ted  and  compared . 


ALUM 


STEEL 


105 


If  the*  load  cell  attach  holts  are  removed  and  the  whole  system  is  heuti-u, 
the  ver t i rn  1  shift  of  the  motor  centerline  is: 

-  (!I  ot^  h  A  t.  ( iqi) 

where 

H  =  initial  height  of  steel  flexure  system 
h  =  initial  height  of  aluminum  motor  mount  assembly 
=  thermal  coefficient  of  expansion  of  steel 
=  thermal  coefficient  of  expansion  of  aluminum 
At  =  temperature  change 

Tlie  vertical  shift  of  a  point  of  the  concrete  originally  on  a  level 
with  the  motor  centerline  is: 

=  (H  ^  h)  At  (192) 

Using  the  following  values  (from  the  Handbook  of  Chemistry  and  Physics, 
lObO-bl) 

^s  ^  ^  in,  in.  degrer  ('  (l.b  \  10  i n . /i n .  degree  V) 

Or  =  2^.0  X  10  in.  in.  degrer  ('  (l“>.3  x  10  i n ./  in./ degree  F) 

“  10  -  X  10  in.  in.  degree  ('  (b.7  x  10  ^  in./ in. /degree  F) 

and  H  h  ~  12  inches,  H  =  18  inches,  and  H  =  1^  inches,  and  using  the 

average  value  of  ,  for  a  2“)  F  temperature  rise  there  is  obtained: 

F.  =  0,007  inch 
M  -  r 

K  =  0.0017  inch 
c 

Therefore,  the  indicated  differential  expansion  of  0.0011  inch  will 
result  in  a  slight  up-bending  of  the  local  cell,  producing  a  transverse 
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SYSTEM  AND  COMPONENT  DESCRIPTION 


The  thrust  total  impulse  measurement  system  is  generally  described  in  two 
ways,  functionally  and  constructionally  (Fig* 13  through  I5  ). 


FUNCTIONAL 

The  solid-propellant  motor,  when  fired,  produces  a  thrust  force  which 
compresses  a  strain-gage  load  cell,  producing  a  d-c  voltage  which  is 
directly  and  accurately  proportional  to  the  thrust  level.  This  d-c  volt- 
iige  equivalent  of  the  motor  thrust  controls  the  generation  of  a  train  of 
pulses  with  a  repetition  rate  directly  proportional  to  the  d-c  voltage, 
and,  therefore,  the  motor  thrust.  An  electronic  counter  records  the  num¬ 
ber  of  pulses  produced  during  the  propellant  burning  time,  and  the  total 
pulse  count  is  directly  proportional  to  the  thrust  total  impulse*  of  the 
burned  propellant.  The  calibration  of  the  entire  electromechanical 
system  is  accomplished  by  using  hydraulic  pressures  produced  by  accurate, 
remotely  operated,  dead  weights  and  supplied  to  a  precision  piston  to 
produce  accurate  static  forces.  The  transmission  of  these  forces  to  the 
motor  support  is  accomplished  by  two  long  tie  rods  extending  through  a 
concrete  retaining  abutment.  The  total  impulse  calibration  of  the  entire 
system  is  accomplished  by  counting  the  number  of  pulses  generated  in  one 
second  by  a  specific  calibration  force  on  the  motor  support  table.  The 
calibration  constant  is  obtained  by  dividing  the  number  of  pulses  gener¬ 
ated  in  one  second  by  the  value  of  calibration  force.  To  obtain  the  most 
accurate  use  of  the  measurement  system,  calibrations  are  made  immediately 
before  and  after  each  motor  firing. 


CONSTRUCTIONAL 

The  general  configuration  of  the  mechanical  system  and  components  of  the 
20-  by  3-*  by  5-foot,  22,000-pound,  two-piece  test  stand  is  shown  in 
Fig.  13  .  The  large,  13-foot,  18,000-pound  section  consists  of  the 
following. 
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Figure  14.  F.lectronic  Kquipoent  Cabinet 
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BALDWIN -LIMA -HAMILTON 
ULTRA  PRECISION  LOAD  CELL 
TYPE  C3P2B  special 
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Hlectricnl  Block  Diafirnm,  Sol  id-Propcllnnt  Toinl  Impulse  Measurement,  Syst  em 


Motor  Support  System 


The  solid-propellant  motor  is  directly  supported  by  two  adjustable  aluminum 
support  blocks  that  enable  alignment  with  the  load  ceil  axis*  The  support 
blocks  are  supported  by  a  welded  aluminum  I-beam  and  plate  assembly  which 
includes  a  braced  thrust  plate  in  direct  contact  with  the  head  end  of 
the  motor.  The  entire  assembly  is  supported  by  a  steel  compression-type 
para  1  lei ograra  flexure  system  w'hich  is  compliant  in  the  thrust  direction 
and  rigid  in  and  about  the  vertical  and  transverse  coordinate  axes.  Level¬ 
ing  jacks  at  the  bottom  of  the  flexures  enables  the  initial  alignment  of 
the  motor-mount  assembly  during  construction  of  the  system. 


Thrust  Dome  Assembly 

-  -  ■»  - 


Motor  thrust  forces  are  conducted  lo  and  away  from  the  load  cell  trans¬ 
ducer  by  a  conical-steel  load  collector  and  a  load  distributor,  respect¬ 
ively.  Symmetrical  mini mum- size  slots  enable  the  installation  and  re¬ 
moval  of  propellant  igniters  in  the  motor  head. 


Concrete  Abutment  and  1 -Deara  Frame 


Motor  thrust  forces  are  transferr«»d  to  a  w,»lded  I-beam  supporting  frame 
!)y  a  by  by  ^-foot  ,  1  1  , lOO-pounil  rnnforced-concrete  block.  Forces 
enter  the  block  through  111-  by  111-  by  2-1  2-inch  imbedded  steel  plates 
and  are  transferred  to  the  supporting  frame  through  a  2-inch-thick  steel 
base  plate.  Reinforcing  rods^"^  Vinch  diameter,  are  spaced  at  b-inch 
intervals  from  top  to  bottom  and  are  welded  to  the  base  plate.  The 
10- inch-high ,  I-beam  frame  is  constructed  of  three  lengthwise  sections 
welded  together  by  I-beam  cross  beams.  ?ix  tiedow^l  bolts  connect  this 
frame  to  the  lest  pod. 


113 


Cali  bra  Mon  fMston  A^Hemhly  and  T  n*  Uod^ 

The  hydraulic  piston  assembly  is  supported  from  the  concrete  abutment 
and  supplies  coaxial  calibration  forces  to  the  motor-mount  assembly  by 
means  of  two  rods  extending  through  clearance  holes  in  the  concrete 
ahu tment . 

The  smaller,  7- ^^oot-long ,  4000-pound,  section  of  the  two-piece  test 
stand  consists  of  the  following: 


Master  Pressure  Standard  (Dead  Weight  Tester) 

This  unit  produces  the  accurate  hydraulic  pressures  transmitted  to  the 
calibrating  piston.  The  production  of  the  pressures  is  accomplished  by 
the  placing  of  dead  weights,  with  remotely  controlled  lifters,  on  a 
piston  in  a  hydraulic  cylinder. 

A  hydraulic  pump,  operated  either  by  a  hand  wheel  or  an  electric  motor, 
maintains  the  proper  supply  of  oil  in  the  closed  hydraulic  system. 


Veather-Tight  Enclosure 

The  precision  calibration  system  is  protected  from  dirt  and  water  by  an 
air-  and  water-tight  enclosure  containing  two  sealed  access  panels  on 
each  side. 

The  genera  1  configuration  of  the  electrical  system  and  its  cotcponents 
is  shown  in  Fig,  14  ,  The  enclosure  contains  the  load-cell  pow'er  supply 

and  bridge  balance,  the  integrating  digital  voltmeter  which  performs  the 
electronic  integration  of  the  thrust-time  curv'e,  control  equipment  for 
the  operation  of  the  ca 1 ibrati on  system ,  and  equijKnent  which  enables 
satisfactory  tape  recording  and  playback  of  the  pulses  produced  by  the 
voltage- to- frequency  converter  section  of  the  integrator. 


REALISTIC  ERRORS 


A  realistic  error  breakdown  of  the  entire  measurement  system  is  presented 
in  Table  3  . 

TABLE  3 

REALISTIC  ERROR  BREAKDOWN 
(Entire  Measurement  System) 

Load  Cell 

Linearity,  ^  0.03 

Hysteresis,  ^  0.02 

Repeatability,  0.01 

Temperature  (lO  F) ,  ^  0.02 

Bridge  Balance 

Resolution  (ipv),  ^  0.02 

Power  Supply 

Stability  (2  C)  %  0.01 

Time-drift,  0.01 

Integrating  V^oltmeter 

Linearity,  ^  0.005 

Stability,  ^  0.01 

Calibrator 

Accuracy,  %  0.05 

Mechanical  Misalignment 

Total  Error,  ^  0.01 

D\'namic  Response  to  Rapid  Transients 

Total  Error,  ^  0.01 

Root-sum-square  error,  less  than  0.0° 
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LOAD  CVU. 


The  preciae  de  terra  i  na  ti  on  of  the  motor  thruBi^tiine  curve  depends  dvrectly 
on  the  quality  of  the  tranBdii<‘er  which  converts  motor  thrust  to  an  <‘qn  i - 
valent  electrical  voltage.  The  need  for  the  greatest  possible  accuracy 
and  linear! ty  of  th i a  transducer  is  evidenced  by  the  over- a  1 1  system 
error  limit  of  0.1^  in  the  measurement  of  the  area  under  the  HATES  motor 
thrv^st-time  curve  (total  impulse).  Efforts  to  fulfill  this  need  to  the 
greatest  possible  degree  led  to  direct  technical  discussions  wi th  the 
Baldvdn-Lima-Hami Iton  Corporation,  Waltham,  Massachusetts.  The  result 
of  thivse  technical  negotiations  was  a  proposed  design  of  a  load  cell 
with  an  accuracy  exceeding  that  of  the  Baldwin-Lima-Hami Iton  precision 
units.  The  performance  specifications  are  as  follows: 

Model:  Special  r“5P2B  (double  bridge) 

Full  scale:  10,0()0-pounds  compression 
Full  scale  output:  "J-mv/v  input  (±0.155^) 

Input  resistance:  3^0  ohms  nominal 

Nonlinearity  (terminal  method):  less  than  0.03^  full  scale 
Hysteresis:  Less  than  0.02i  full  scale 

Repeatability:  Less  than  0.01%  full  scale 

Tempera ture  effects; 

1.  On  zero  output:  less  than  0,\^%  full  scale/lOO  F 

2.  On  sensitivity:  less  than  0.08^  of  load  per  100  F 
Full-scale  deflection:  less  than  0.000  inch 

Important  additional  calculated  data  of  this  unit  were  obtained  con¬ 
cerning  mechanical  stiffness  to  transverse  toads,  bending  and  twisting 
moments,  and  are  as  follows: 

Mechanical 


1,  The  transverse  stiffness  is  in  the  order  of  1  x  10  Ib/in.  in 

any  direction  for  a  load  applied  at  the  top  of  the  ^ell  (3  inches 
f rora  the  base). 
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2.  The*  iorsienal  stiff  >e08  is  in  the  ordpr  of  20  x  10^  lb  in.,  rad 

3.  The  bendirift  stiffness  for  a  couple  applied  to  the  end  of  the 
cell,  tending  to  bend  the  axis  of  the  cell,  is  in  the  order  of 
25  X  10  lb  in.  rad  in  any  direction 


Electrical 


1.  The  cell  response  to  u  500-pound  transverse  load  is  a  maximum 
span  (sensitivity)  change  of  M),i5^  of  full  scale,  and  this 
change  increases  and  decreases  sinusoidally  as  the  direction 
is  rotated  about  the  cell, 

2.  Cell  response  to  torsion  will  be  negligible 

5.  A  moment  of  2500  lb  in.  applied  at  the  end  of  the  cell  has  the 
same  effect  as  item  1. 

These  mechanical  and  electrical  data  enabled  an  adequate  design  analysis 
to  be  made  of  a  suitable  motor-suspension  flexure  system.  Figure  16 
shows  the  dimensions  of  the  load  cell  and  the  detail  of  the  special 
heavy  end  mounting  plates  containing  alignment  center  holes  and  threaded 
holes  for  bolted  a ttcchmen^  s . 


POWER  SUPPLY 

The  low  over-all  system  error  limit  of  0.1^  in  the  measurement  of  thrust 
total  impulse  imposes  severe  requirements  on  the  load  cell  bridge  ex¬ 
citation  voltage.  Mea su rement -sys t em  accura<y  is  directly  related, 
through  the  load  cell  sensitivity,  to  the  power- supply  stability.  Ef¬ 
forts  to  obtain  a  suitable  strain  gage  f>ower  supply  led  to  technical 
negotiations  with  the  Systems  Research  ('orperat i on ,  Van  Nu^s,  California. 
The  result  of  these  discussions  was  a  proposed  new  design  of  power  supply 
vi tb  per f ormance  specifications  as  follows: 

Model :  5512 

Input  power;  95-155  vne,  'lO  cycle 


3/8  DIA  X  1/2  DEEP 

PILOT  HOLE - 

TYP  BOTH  ENDS 


Figure  16  .  load  fell 


Output;  tally  Ofo  IH  vd<* ,  0  to  UOO  rn  i  1  1  iamfjpr<*8 

Rp^ij  I  a  I  i  or) :  0,01*^  NL-FI.  and  O'!  to  Ml  vac 

Ripplo:  O.O*)  millivolt  rma 

dr  isolaiioii;  10,0(H1  nirgolims 

'■.1  a  1)1  i  1 1  y :  0 . 001^  C 

Noi  sp  at  bridg:r  (  '5"){)  ) :  2IJ\ 

Ovp r 1 oad  protection:  tuscd 

Overvoltage  protection:  limits  output  voltage  to  iH  vdr 
Current  regulation :  ova i lah le  by  jumper i ng  a t  rea r  rhass i h 
Size:  1-1  2  -  by  10-  by  1  0-  i  ru  be  s 

This  supply  has  a  solid-state  ronstrurtioii  incorporating  balanced  differ¬ 
ential  amplifiers  and  temperature  compensated  >.ener  circuitry  for  sta¬ 
bility.  Output  voltage  is  digitally  selected  in  1-millivolt  steps  by 
means  of  decade  switching  of  precision  voltage  dividers  connec  ted  in  a 
Kelvin-Varley  arrangement.  By  accurate  selection  of  close- tolerance 
resistors^  superior  accuracy  io  obtained  and  resolution  to  five  signi¬ 
ficant  f igures  or  one  pa rt  in  100,000  is  more  thiin  adequate  for  the  ap¬ 
plication.  Two  excitation  outputs  are  available  ai  the  rear  terminals. 
Isolation  to  ground  is  extremely  high  as  a  result  of  design  (’are  in  the 
layout  of  the  circuity  components  and  the  unit  io  mode  completely  float¬ 
ing  with  reference  to  ground.  I/esign  <are  is  also  taken  to  keep  capaci- 
M  VC  ‘.oupl  ing  to  a  minimum .  Fxre  1  lent  i  so  lut  i  cn  to  lino  is  achieved  h> 
using  transformers  with  a  (ompletrly  box-shielded  primary.  Special 
features  of  overload  protection  are  included  to  limit  th(»  supply  output 
to  18  volts,  regardless  of  (onditions,  as  a  result  of  zener  diode  effect¬ 
ing  the  feedback  loop  to  the  regulating  series  pass  transistors.  Low 
ripple  is  achieved  by  heavy  negative  feedback  coupled  with  preregulator 
filtering.  The  digital  programming  portion  of  tliese  units  provi^Jec  a 
precision  operation  and  excellent  application  where  exact  power  settings 
are  required.  Operation  is  simple  because  the  voltage  desired  appears 
immediatel  y  at  the  rear  termipais  after  selecting  the  five  front  panel 
knobs.  This  is  a  distinct  advantage  over  a  continuously  variable  control 
because  it  takes  time  to  reach  the  voltage  setting  desired  as  well  as 
expensive  additional  rc^adout  equipment.  The  additional  advantage  is 
that  any  specific  value  of  voltage  can  be  selected  giving  operational 


gimplifity  and  p  I  imi  no  1 1  on  of  prror.  This  fpaturo  of  exiellpnt  reprat- 
ah  1 1 1 ty  i g  hap f u 1  in  appl i ra  t i ons  of  re pen  ting  test  condi ti ons . 

Figure  17  shows  the  general  front  panel  configuration. 


BRIDGE  BALANCE 

The  over-all  system  error  limit  of  0.1^  also  itnpose^^  severe  ’Require¬ 
ments  on  the  circuitry  employed  for  the  load  cell  tput  zero  signa  1 
control.  Bridge  balance  must  be  accomplished  with  a  resolution  of 
approximately  ]  if  load  cell  signal  ,»rr^rff  of  0,01^  are  not 

to  be  exceeded  with  15t>0-pound  thrust  levels.  The  Systems  Research 
Corporation,  Van  Nuys,  California,  proposed  a  design  of  an  adequate  unit 
with  performance  and  description  os  follows; 

Model  4700 

Input:  Three  wires  to  bridge 

Resolution:  Balance  to  microvolt  with  350-ohra  bridge 
Size:  3-1/2  by  19-  by  6  inches 

Six  front  panel  mounted  controls  will  allow  nulling  bridge  unbalances 
to  ^55^-  Precision-tolerance  resistors  of  excellent  long-term  stability 
are  used  throughout.  Quality  switches  of  low-contact  resistance  ensure 
excellent  repeatability.  As  in  the  power  supply,  the  system  used  is 
decade  switching  of  precision  voltage  dividers  connected  in  a  Kelviii- 
Varley  arrangement.  Terminals  are  provided  for  a  limiting  resistor  if 
desired,  otherwise  these  terminals  are  jumpered. 

Figures  18  and  19  show  the  front  panel  configuration  and  circuit  de¬ 
sign  and  respectively. 

INTEGRATING  DIGITAL  VOLTMETER 

The  Beasureaent  of  the  are©  under  a  thrust-time  curve  is  readily  accomp¬ 
lished  by  the  use  of  a  voltage- to-frequency  converter  and  a  counter,  A 
thrust  (F)  on  the  load  cell  produces  a  dc  voltage  (V)  into  the  converter. 


Power  Supply 


EXCITATION 


Fi^ur<‘  10, 


lialariro  ('irciiif  I)<»Rign 


which  in  turn  produces  a  pulse  repetition  frequency  (f),  and  the  indivi¬ 
dual  prises  are  electronically  counted.  As  a  result,  the  total  number 
of  counts  (n)  is  proportional  to  the  time  integral  of  the  force  (total 
impulse).  Mathematically,  this  is  expressed  as 


N 


dt 


Cl  Cg 


F  dt 


therefore 


where 


and 


C  =  Cg 

i.‘b 

It  - 


To  comply  with  an  over-all  system  error  limit  of  0.1^  in  the  total  im¬ 
pulse  measurement,  the  linearity  and  stability  of  the  integrating  equip¬ 
ment  must  be  of  the  highest  possible  degree.  Also,  th«?  total  number  of 
counts  obtained  with  the  expected  minimum  thrust  (1500  pounds)  during 
an  expected  minimum  burning  time  interval  (2  seconds)  must  be  suffic- 
ently  great  so  that  one  count  represents  a  negligible  error  when  com¬ 
pared  to  O.ljC  of  the  total  number  of  counts  obtained.  If  one  count  is 
required  to  be  only  O.OlJt  of  the  total  number  obtained  (with  the  minimum 
thrust  and  time  conditions),  then  the  total  number  must  be  at  least 
10,000  or  5000  pps.  Consequently,  any  converter-counter  combination 
^ich  yields  5000  pulses  per  second,  with  the  minimum  thrust  levels  to 
be  encountered,  is  adequately  sensitive  for  the  application.  Linearity 
of  the  converter-counter  coTibination,  like  the  load  cell  transducer, 
must  be  the  greatest  possiM*  because  nonlinearities  occurring  anywhere 
in  the  system  produce  errors  ij  the  total  impulse  measurement. 
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The  I^ec  Model  inr*2401  A  integrating  digital  voltmeter,  manufactured 
by  the  Dymec  Division  of  Hewlett-Packard  Company,  fulfills  these  re¬ 
quirements  when  used  in  conjunction  with  the  special  Baldvin-Lima- 
Hamilton  load  cell.  Its  superior  linearity  (O.OO^jt  is  obtained  at  the 
expense  of  sensitivity  (O.l  vdc  yields  100,000  pps),  but  because  the 
special  Baldvin-Lima-IIamil ton  load  cell  produces  14  millivolts  with 
10,000-pounds  force,  its  sensitivity  is  sufficient  to  produce  the  de¬ 
scribed  negligible  error  (O.Ol^  with  minimum  thrust  and  time  conditions) 
with  the  one-count  resolution  common  to  high-quality  electronic  counters, 
as  the  following  statements  demonstrate: 

vco  sensitivity  is  100  mi  1  livolts— ^100,000  pps 

load  cell  10,000— ►  I'l  millivolts - ^54,000  pps 

1500— ►S.l  mill  ivolts  ■  ^  8100  pps 

O.OlJt  of  1500— ►8,1  •  10"^  millivolts - ►O.S  pps 

Therefore,  in  2  seconds  the  total  count  is  1.6.  Consequently,  the 
Baldwin-Liraa-Hamilton  load  cell  and  Dvmec  integrating  digital  voltmeter 
combination  sensitivity  and  resolution  are  consistent  with  the  assigned 
error  limitation  of  O.OljC. 

The  performance  features  and  sped  fi cut  ions  of  the  Djinec  2401  A  unit 
include  the  following; 

1.  Better  than  140-dec  ibel  c*ffe<  tive  lommoii  mode  rejection  at  all 
frequencies  ,  inc  hiding  d<‘ 

2.  Active  integration  minimi/es  error  caused  by  superimposed  noise 
and  provides  average  reading  of  input  over  sampling  period 

3.  Five-digit  readout  provides  up  to  one  part  in  100,000  reading 
resolution.  Minor  changes  accommodate  10^  counts. 

4.  0.005,^  linearity,  stability  ±0.05^  of  full  scale  per  day  on 
most  sensitive  range  (0, 1  vdc)  with  constant  temperature  and 
±10^  line  voltage  change 

5.  Temperature  sensitivity  ±0,002^^  of  reading  per  degree-centigrade 
on  the  most  sensitive  range  (O.l  vdc) 
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Figuii*  flhow.M  fb<*  front  fiunel  ronf  i  giira  1 1  on  . 


(^MJimATHlN  SYSTtlM 

From  a  mou Hiirpraoii t  arc  ura<  y  st.ainipoint  ,  the  most  important  Hinkle  (‘omi)Ofi- 
ent  of  the  entire  measurement  system  is  the  equipment  uaed  for  ealihra- 
t  1  on .  The  aerurary,  linearity,  repeatability,  stability,  and  all  other 
instrumentation  per f ormarue  parameters,  ran  be  determined  only  by  reference 
to  an  accepted  standard.  The  over-all  error  of  the  system  <an  be  no  less 
than  the  error  incurred  in  its  <  al  itiru  tion,  Efforts  to  obtain  the  most 
precise  force  standard  possible  for  this  purpose  led  to  direct  technical 
negotiations  with  the  Ruska  Instrument  Corporation,  Houston,  Texas*  The 
result  of  these  discussions  was  a  proposed  design  for  a  combined  master 
pressure  standard  and  a  precision  hydraulic  piston  for  the  production  of 
six  forces  in  the  0  to  10,000-pound  range.  To  obtain  the  most  accurate 
calibration  of  this  equipment  possible  it  was  also  proposed  that  this 
equipment  be  calibrated  by  the  National  Bureau  of  Standards,  Washington, 

I).  C. 

The  operating  principle  of  this  equipment  is  the  production  of  accurate 
hydraulic  pressures  by  means  of  prei ision  dead  weights  supported  by  a 
vertical  piston,  and  the  transmission  of  these  pressures  to  a  second 
precision  piston  oriented  to  as(*  the  accurate  output  forces  with  the 
equipment  to  be  calibrated.  The  application  of  this  principle  to  the 
ca 1 ibra ti on  of  a  horizontal  thrust  measuring  system  constitutes  the  most 
accurate  method  known  for  this  purpose,  because  it  eliminates  all  pivots, 
bell  cranks,  bearings,  flexures,  lever  arms,  etc,,  necessary  to  produce 
horizontal  forces  with  dead  weights  alone.  The  operational  flexibility 

of  this  equipment  in  the  Hates  motor  test  stand  enables  remotely  controlled 
calibrations  of  the  entire  measurement  system  to  be  made  immediately  be¬ 
fore  and  after  a  motor  firing,  with  11  controls  mounted  on  a  single  panel 
and  located  with  the  remaining  electronic  equipment  of  the  system.  Th' 
calibration  equipment  control  wiring  will  enable  operating  it  from  either 
che  test  stand  or  the  control  room. 
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&1  'oltmeter 


The  specifications  of  the  equipnent  include: 


1.  Six  calibration  points:  1000-,  11000  ,  ^*000-,  (>000-,  HOOO-,  and 
10 ,000-pound8  force  (dead  weights  raonufact ured  to  the  lo<al  g 
value  at  Kdwurds  Air  Force  liase  test  site). 

2 .  Accuracy  of  each  force  po ini:  tO. 0^^ 

Resolution:  0.02^  of  full  scale 

^1.  Float  time  of  dead  weight  gage  and  force  piston  cylinder  shall 
bp  at  least  1  minute  for  each  calibration  point 

'i.  Hydraulic  pump  capacity  shall  be  sufficient  to  provide  at  least 
one  complete  calibration  (ascending  and  descending),  including 
any  leakage ,  without  ref i lling 

().  Visual  means  shall  be  provided  at  the  control  console  to  indicate 
when  the  dead  weight  gage  is  floating  and  a  calibration  point 
may  be  taken 

7.  Weights  shall  be  applied  in  ascending  or  descending  order  select¬ 
able  through  a  multiposi tion  switch 

8.  The  piston  cylinder  of  both  the  dead  weight  gage  and  force  piston 
cylinder  shall  be  of  the  same  material. 

Figure  21  shows  this  oquipinent. 
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Figure  21^  Calibration  System  Control  Panel 


CALCUUTION  OF  COMPIiETF  SPRING  SYSTEM 


The  fflaihematicel  analyses  described  in  preceding  discussions  provided 
the  initial  justification  of  the  detaih‘d  design  of  the  entire  measure- 
aent  system.  However,  practical  consideration  of  system  and  component 
fabrication  and  assembly,  system  calibration,  and  motor  installation  and 
operation ,  require  some  deviation  f rom  the  simple,  i dea 1  image  of  the 
system  as  a  1'  or  2-degree-of-freedom  system.  These  deviations  were 
held  to  a  minimum,  however,  not  only  to  design  the  syst.em  as  close  as 
possible  to  the  systems  which  are  more  amenable  to  mathematicr4 1  analysis 
but,  more  important,  because  the  simpler  the  system  the  fewer  sources  of 
error . 

It  is  useful,  therefore,  to  estimate  the  stiffness  and  weights  of  the  de¬ 
signed  components  so  that  a  more  accurate  picture  may  be  obtained  of  the 
dynamic  system.  The  purpose  in  such  an  effort  is  to  (l)  determine  how 
close  the  designed  system  came  to  the  simpler  models  used  in  the  mathe¬ 
matical  analyses,  (2)  provide  a  more  accurate  basis  for  possible  addi¬ 
tional  analyses,  if  needed,  and,  (l)  identify  possible  sources  of  unde¬ 
sirable  dynamic  performance  which  ran  be  corrected  by  minor  changes  dur¬ 
ing  fabrication  and  assembly. 


calculation  of  component  STIFfWlSS 


Concrete  Abutment 


fh@  StflectioQ  of  th@  concrete  ebutnent,  fised  et  the  base,  is  obtained 
bj  the  use  of  the  folloving  equation  (Ref.  8  ): 


kl 

AE 
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^ere 

6  m  deflection  of  point  0  relative  to  base 

P  «  static  load  (10,000  pounds) 

E,Eg  concrete  nodnlus  of  elasticity  (3  x  10  psi)  and  shear 
(1.35  *  10^  psi) 

*X 

I  s  aoasnt  of  inertia  (l/l2  wh^) 

k  m  constant  (3/2  for  rectangular  cross  sections) 

a  length  of  bean  (32  inches) 

A  m  cross  section  parallel  to  base 


Substitution  ofthese  values  into  the  equation  yields  the  value  of  stiffness 
o2  the  concrete  abutment. 
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7.27  X  10® 
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Because  load  cell  stiffness  is  k 
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The  stiffness  ratio  is 
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Load  Diatributor  (Steel 


i 


Conical  Part. 


Area  of  elemental  croas  hatched  circular  alice  dA  -  frr 


y  P/dA  F  dx  F  ^ 
dO/dx  "  dFd®  ~  d?  ^^2 


»  ■ 


dx 

2 


The  total  axial  deflection  (6)  of  the  conical  part  associated  with 
force  (p)  is 
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Cylindrical  Part, 
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Cross-sectional  area  is 


=  Ft 
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TotsI  dsflsctioa  of  load  distributor  obtained  with  force  (F)  is 


a  ^ 


3F _  H  ^  p.  /  \ 

10  ^  EA  “  '10^  ■*’EA' 


Substituting  the  appropriate  values  yields  the  stiffness  value 

“u-  ■  rrs-  -  ‘ 

The  stiffness  ratio  is 


336  X  10^ 
2.04  X  10^ 
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Load  Collector  (Thrust  Doae) 

The  stiffness  of  this  conponent  was  estimated  to  be  approximately  l/4 
the  stiffness  of  the  load  distributor.  Therefore 
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A  40 


Rocket  Motor  Case 
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d<,=  12,  7S 


£  =20.2''- 
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>  XT  4  (20.2) 


=  25.2  X  10^  Ib/ft 


Stiffness  ratio  is 


k 


25. g  g  10^ 
2.04  %  10^ 


12.2 
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DIAGRAM  OF  C(MPIETE  SYSTEM 


On  the  basis  of  the  estimated  component  stif ,"*nes8es  and  weights,  the 
complete  spring-maes  system  is  given  below.  Damping  values  are  those 
regarded  as  nominal  in  the  computer  program. 


^2  > 


The  model  (nominal)  used  in  the  computer  program  was: 


FCt) 


The  two-degree-of-freedom  model  used  in  the  computer  analysis  is  a  simpli¬ 
fication  of  the  actual  system.  Lumped  masses  were  used  in  the  model 
which  were  comparable  to  the  distributed  masses  of  the  concrete  abutment 
and  solid  propellant  motor. 


The  load  distributor  and  collector  springs  were  omitted  from  the  computer 
model  altogether,  and  the  distributed  spring  and  masses  of  the  motor  were 
simply  lumped  as  a  mass.  Nevertheless,  it  is  felt  that  the  two-degree- 
of-freedom  computer  model  represented  the  essential  features  of  the  system, 
because  the  computer  results  showed  the  impulse  error  sensitivity  to 
abutment  stiffness  to  be  very  slight.  Therefore,  the  unification  of  the 
load  distributor  and  abutment  springs  is  entirely  reasonable.  The 
correctness  of  the  unification  of  load  collector  and  motor  springs  and 
masses  into  one  single  mass  (^00  pounds)  is  not  so  apparent.  To  further 
complicate  the  assessment  of  the  computer  model  correctness,  the  motor 
thrust  is  not  all  applied  at  one  location.  Approximately  of  it 
(Fl  (t)  in  the  detailed  diagram)  is  applied  at  the  head-closure  end  of 
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the  sotor  case,  and  the  Training  3^  (Fg  (t)  in  the  diagram)  is  applied 
at  the  noaale  asa^bly  and  it  is  traneisitted  to  the  head  end  through 
the  ct^paratively  eompliant  rocket  motor  case  ( 12  times  the  stiffness 
of  the  load  cell).  Consequently,  the  dynamic  behavior  of  the  part  of 
the  actual  system  beti^een  the  load  cell  and  the  motor  nossle  could  con^ 
eeivably  be  quite  different  from  that  of  the  limped  mass  representing 
it  in  the  computer  studies.  However,  most  of  the  misrepresentation  is 
associated  with  the  motor  hardware,  which  is  outside  the  test  stand 
design  and  therefore  beyond  Rocketdyne  design  control.  It  is  felt  that 
the  dynamic  behavior  of  the  motor  hardware  should  receive  additional 
attention  in  the  future  evaluations  of  the  constructed  measurement 
system. 
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ADDITIONAL  INSTOlIMflNTATrON 


VCO  PULSE  MODIFI01  FOR  TAPE  PLAYBACK 

Because  a  standard  tape  rt^corder  cannijt  record  and  playback  pulses  of 
short  duration  vith  good  fidelity^  it  is  neccHsary  fo  modify  the  pulse 
characteristic  at  the  input  and/or  output  of  the  recorder.  The  reason 
for  this  is  that  the  upper  frequency  I  irai  t  (j  f  the  recorder  is  near  t  he 
natural  ringing  frequency  of  the  head,  whicti  causes  the  system  to  ring 
electrically  when  activated  by  a  pulse,  Therefore,  when  a  pulse  is 
recorded,  the  playback  of  that  pulse  wi  1 1  fje  a  damped  oscillatory  wave 
with  a  frequency  equal  to  the  ringing  frequency  of  the  system,  and  the 
rate  of  decay  (damping)  determined  by  the  Q  of  the  system.  Tf  the  play- 
back  pulse  is  being  read  by  a  counter,  it  is  quite  possible  to  obtain 
more  than  one  count  per  initial  pulse  because  the  excursions  of  the  im¬ 
mediate  subsequent  ringing  can  be  of  sufficient  amplitude  to  activate  4i  e 
counter*  The  amplitude  resolution  of  the  counter  cannot  be  predictably 
set  to  discriminate  between  the  first  and  the  second  cycles  of  the  damped 
wave  because  a  slight  change  in  playback  pulse  amplitude  will  cause  the 
counter  to  either  miss  or  to  produce  multiple  counts. 

A  possible  method  of  solving  this  prohl<*m  has  been  tested  at  Rocketdync 
and  proved  practical.  The  pulses  can  be  first  modified  to  a  constant 
amplitude  and  duration  ofabout  (>  volts  am!  h  microseconds  by  use  of  a 
monostable  multivibrator  (one-shot).  These  pulses  can  he  recorded  on  a 
tape  recorder  with  a  flat  frequency  response  of  at  least  100  kilocycles 
such  as  an  Ampex  FR-100  used  on  direct  record  mode.  Figure  2UA  shows  a 
train  of  three  pulses  as  they  would  be  fed  to  the  recorder,  and  Fig.  22B 
shows  the  approximate  characteristic  of  the  playback.  Figure  214  shows 
one  of  the  playback  pulses  expanded  on  a  shorter  time  base.  This  signal 
can  be  fed  to  two  parallel  circuits  as  shown  iri  Fig,  2k  ( t.he  10  kilocycle 
hlpass  filter  is  present  to  remove  most  of  the  noise).  The  signal  of 
Fig.23A  can  be  fed  to  one  aide  of  the  summing  amplifier,  and  through  a 
delay  network  to  the  other  side.  In  this  way  the  two  signals  which  will 
be  added  together  look  like  the  wave  forms  ot  Fig.^TA  and  The  delay 
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FifKir&  22.  Tap®  Recorder  Tnpnt  and  Played^bsek  Signals 


can  be  adjusted  until  180  degree  phase  shift  is  produced  (for  100  kilo¬ 
cycles  this  will  be  5  fflicroseconds)  and  the  aaplitude  adjusted  until 
the  first  half-cycle  of  the  delayed  signal (23B)  is  equal  and  opposite 
to  the  second  half-=*cycle  of  the  original  signal.  V^hen  these  conditions 
exist,  all  succeeding  half-cycles  will  cancel  and  the  sum  of  the  tvo 
signals  (Fig,23C)  will  be  zero  except  for  the  first  half-cycle,  thus 
effectively  recreating  a  singl *  pulse.  (Some  small  amplitude  trash  will 
remain  as  a  result  of  harmonics,  but  will  be  negligible.)  The  delay  net¬ 
work  must  be  adjustable  near  5  microseconds  and  some  amplitude  control 
is  necessary  because  no  two  tape  recorders  will  produce  the  same  wave 
form  on  playback.  However,  once  set  for  any  particular  recorder,  no 
further  adjustment  should  be  needed  unless  the  system  is  to  be  used  with 
another  recorder.  Figure  25  shows  the  general  appearance  of  this  unit. 
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Nat'l  Aeronautics  it  Space  Admin.  1 
Washington,  D.  C.  20546 
Attn:  R.  W.  Ziem  (RPS) 

Nat'l  Aeronautics  it  Space  Admin.  1 
Goddard  Space  Flight  Center 
Greenbelt,  Maryland  20771 
Attn:  Library 

Defense  Documentation  Center  20 

Cameron  Station 
Alexandria,  Virginia  22314 

RTD (RTNP)  1 

Bolling  AFB 
Washington,  D.  C. 

Arnold  Eng.  Development  Center  1 
Attn:  AEOIM 

Air  Force  Systems  Command 
Tullahoma,  Tennessee  37389 

AFRPL  (RPR) 

Edwards,  California  92523 
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AFRPL  (RPM)  1 

Edwards,  California  92523 

AFFTC  (FTBAT-2)  1 

Edwards  AFB,  California  93523 

Office  of  Research  Analysis  (OAR)  1 
Attn:  RRRT 

Holloman  AFB,  New  Mexico  88330 


AFRPL  (RPFDE)  33 

Edwards,  California  92523 

Air  Force  Missile  Test  Center  1 

Attn:  MTBAT 

Patrick  Air  Force  Base,  Florida  32925 

Wright-Patterson  AFB,  Ohio  1 

45433 

Attn:  AFML(MAAE) 

Commanding  Officer  1 

Ballistic  Research  Laboratories 


Aberdeen  Proving  Ground,  Maryland 

Attn:  AMXBR-1 

21005 

Commanding  Officer  1 

Harry  Diamond  Laboratories 
Washington,  D.  C.  20438 
Attn:  AMXDO-TIB 

Commanding  Officer  1 

U.S.  Army  Research  Office  (Durham) 
Box  CM,  Duke  Station 
Durham,  North  Carolina  27706 

Commanding  Officer  1 

U.S.  Army  Chemical  Research  Sc 
Development  Laboratories 
Edgewcod  Arsenal,  Maryland  21010 
Attn:  Chief,  Physicochemical 
Research  Division 

Commanding  Officer  1 

Frankford  Arsenal 
Philadelphia,  Pennsylvania  19137 
Attn:  Propellant  Sc  Explosives 
Section,  1331 


Commanding  Officer  1 

Picatinny  Arsenal 
Dover,  New  Jersey  07801 
Attn:  Library 

Commanding  Officer  1 

Picatinny  Arsenal 

Liquid  Rocket  Propulsion  Laboratory 
Dover,  New  Jersey  07801 
Attn:  Technical  Library 

U.S.  Army  Missile  Command  4 

Redstone  Scientific  Info.  Center 
Redstone  Arsenal,  Alabama  35808 

Attn:  Chief,  Document  Section 

Commanding  General  1 

White  Sands  Missile  Range 
New  Mexico  88002 
Attn:  Technical  Library 

Bureau  of  Naval  Weapons  2 

Department  of  the  Navy 
Washington,  D.  C.  20360 
Attn.  DLI-  3 

Bureau  of  Naval  V/eapons  2 

Department  of  the  Navy 
Washington,  D.  C.  20360 
Attn:  RMMP-2 

Bureau  of  Naval  Weapons  1 

Department  of  the  Navy 
Washington,  D.  C.  20360 
Attn:  RMMP-3 

Bureau  of  Naval  Weapons  1 

Department  of  the  Navy 
Washington,  D.  C.  20360 
Attn;  RMMP-4 

Bureau  of  Naval  Weapons  1 

Department  of  the  Navy 
Washington,  D.  C.  20360 
Attn;  RRRE-6 


Distribution  List  (Cont'd)  p  3 


Commander  2 

U.  S.  Naval  Missile  Center 
Point  Mugu,  California  93041 
Attn:  Ti  chnical  Library- 

Commander  2 

U,  S.  Naval  Ordnance  Laboratory 
White  Oak 

Silver  Spring,  Maryland  20910 
Attn:  Library 

Comma.nder  5 

U.  S.  Naval  Ordnance  Test  Station 
China  Lake,  California  93557 
Attn:  Code  45 

Commander  (Code  753)  2 

U.  S.  Naval  Ordnance  Test  Station 
China  Lake,  California  93557 
Attn:  Technical  Library 

Superintendent  1 

U.  S,  Naval  Postgraduate  School 
Naval  Academy 
Monterey,  California 

Commanding  Officer  1 

U.  S.  Naval  Propellant  Plant 
Indian  Head,  Maryland  20640 
Attn:  Technical  Library 

Commanding  Officer  1 

Office  of  Naval  Research 
1030  E.  Green  Street 
Pasadena,  California  91101 

Director  (Code  6180)  1 

U.S,  Naval  Research  Lab. 
Washington,  D.  C.  20390 
Attn:  H.  W.  Carhart 

Director  1 

Special  Projects  Office 
Department  of  the  Navy 
Washington,  D.  C.  20360 


Commanding  Officer  1 

U.S.  Naval  Underwater  Ordnance 
Station 

Newport, Rhode  Island  02844 
Attn;  W.  W.  Bartlett 

Aerojet- General  Corporation  1 

P.  O.  Box  296 

Azusa,  California  91703 

Attn;  F.  M.  West,  Chief  Librarian 

Aerojet- General  Corporation  1 

11711  South  Woodruff  Avenue 
Downey,  California  90241 
Attn:  Florence  Walsh,  Librarian 

Aerojet-General  Corporation  3 

P.  O.  Box  1947 

Sacramento,  California  95809 
Attn:  Technical  Information  Office 

Aeronutronic  Div.  Philco  Corp.  1 
Ford  Road 

Newport  Beach,  California  92600 
Attn:  Dr.  L.  H.  Linder,  Manager 
Technical  Information  Dept. 

Aeroprojects,  Inc.  1 

310  East  Rosedale  Avenue 
West  Chester,  Pennsylvania  19380 
Attn:  C.  D.  McKinney 

Aerospace  Corporation  2 

P.  O.  Box  95085 
Los  Angeles,  California  90045 
Attn:  Library- Documents 

Allied  Chemical  Corporation  1 

General  Chemical  Division 
Research  Lab.  ,  P.  O.  Box  405 
Morristown,  New  Jersey  07960 
Attn:  L.  J.  Wiltrakis,  Security  Offic 

Amcel  Propulsion  Company 

Box  3049 

Asheville,  North  Carolina  28802 
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American  Cyanamid  Company  1 

1937  W.  Main  Street 
Stamford,  Connecticut  06902 
Attn:  Security  Officer 

IIT  Research  Institute  1 

Technology  Center 

Chicago,  Illinois  60616 

Attn:  C.  K.  Hersh,  Chemistry  Div. 

ARO,  Inc .  1 

Arnold  Engrg.  Dev.  Center 
Arnold  AF  Stati  1,  Tennessee  37389 
Attn:  Dr.  B.  H.  joethert 

Director  of  Engrg. 

Atlantic  Research  Corporation  2 

Shirley  Highway  &  Eds  all  Road 
Alexandria,  Virginia  22314 
Attn:  Library 

Atlantic  Research  Corporation  1 

Western  Division 

48-506  Easy  Street 

El  Monte,  California 

Attn:  H.  Niederman 

Battelle  Memorial  Institute  1 

505  King  Avenue 
Columbus,  Ohio  43201 
Attn:  Report  Library,  Room  6A 

Bell  Aerosystems  1 

Box  1 

Buffalo,  New  York  14205 
Attn:  T.  Reinhardt 

The  Boeing  Company  1 

Aero  Space  Division 

P.  O.  Box  3707 

Seattle,  Washington  98124 

Attn:  R.  R.  Barber,  Lib.  Ut.  Ch. 

Chemical  Propulsion  Information  3 
Agency 

Applied  Physics  Laboratory 

8621  Georgia  Avenue 

Silver  Spring,  Maryland  20910 


Defense  Metals  Information  Ctr.  1 
Battelle  Memorial  Institute 
505  King  Avenue 
Columbus,  Ohio  42301 

Douglas  Aircraft  Co.  ,  Inc.  i 

Santa  Monica  Division 
3000  Ocean  Park  Boule''’ard 
Santa  Monica,  California  90405 
Attn:  Mr.  J.  L.  Waisman 

The  Dow  Chemical  Company  1 

Security  Section 
Box  31 

Midland,  Michigan  48641 

Attn:  Dr.  R.  S.  Karpiuk,  1710  Bldg. 

E.  I.  duPont  deNemours  &  Co.  1 
Eastern  Laboratory 
Gibbstown,  New  Jersey  08027 
Attn:  Mrs.  Alice  R.  Steward 

Esso  Research  Sc  Engineering  Co.  1 
Special  Projects  Unit 
P.  O.  Box  8 

Linden,  New  Jersey  07036 
Attn:  Mr.  D.  L.  Baeder 

General  Dynamics/ Astronautics  1 
P.  O.  Box  1128 
San  Diego,  California  92112 
Attn:  Library  &  Information 
Ser'''’ices  (128-00) 

General  Electric  Company  1 

Flight  Propulsion  Division 
Evendale 

Cincinnati,  Ohio  45215 

Hercules  Powder  Company  1 

Allegany  Ballistics  Laboratory 
P.  O.  Box  210 

Cumberland,  Maryland  21501 
Attn:  Library 
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Hercules  Powder  Company  1 

Research  Center 
Wilmington,  Delaware  19899 
Attn;  Dr.  Herman  Skolnik,  Manager 
Technical  Information  Division 

Jet  Propulsion  Laboratory  1 

4800  Oak  Grove  Drive 
Pasadena,  California  91103 
Attn;  Library  I.TDS) 

Mr.  N.  E.  Devereux 

Arthur  D.  Little,  Inc.  1 

15  Acorn  Park 

Cambridge,  Massachusetts  02140 
Attn:  W.  H.  Varley 

Lockheed  Propulsion  Company  3 

P.  O.  Box  111 

Redlands,  California  92  374 

Attn;  Miss  Belle  Berlad,  Librarian 

Marquardt  Corp.  1 

16555  Staticoy  Street 
Box  2013  -  South  Annex 
Van  Nuys,  California  91401 

Martin  Company  1 

Baltimore,  Maryland  21203 
Attn;  Science  -  Technology 
Library  -  Mail  398 

Minnesota  Mining  fn  Manufacturing  2 

Co. 

900  Bush  Avenue 
St.  Paul,  Minnesota  55106 
Attn;  Code  001  3  R  &;  D 
VIA:  H.  C.  Zeman 

Security  Administrator 

Monsanto  Research  Corporation  1 

Boston  Labs.  ,  Everett  Station 
Chemical  Lane 

Boston,  Massachusetts  02149 
Attn;  Library 


New  ork  University  1 

Dept,  of  Chemical  Engineering 
New  York,  New  York  10053 
Attn;  P-  F.  Winternitz 

North  American  Aviation,  Inc.  1 
Space  &r  Information  Systems  Div. 
12214  Lakewood  Boulevard 
Downey,  California  90242 
Attn;  W.  H.  Morita 

Rocketdyne  3 

6633  Canoga  Avenue 
Canoga  Park,  California  91304 
Attn:  Library,  Dept-  596-306 

Rohm  &  Haas  Company  I 

Redstone  Arsenal  Research  Division 
Huntsville,  Alabama  36808 
Attn:  Librarian 

Space  Technology  Laboratory,  2 

Inc. 

1  Space  Park 

Redondo  Beach,  California  90200 
Attn;  STL  Tech.  Lib.  Doc- 
Acquisitions 

Texaco  Experiment  Incorporated  1 

P.  O.  Box  1-T 
Richmond,  Virginia  23202 
Attn:  Librarian 

United  Technology  Center  1 

P.  O.  Box  358 

Sunnyvale,  California  94088 
Attn:  Librarian 

Thiokol  Chemical  Corporation  1 

Alpha  Division,  Huntsville  Plant 
Huntsville,  Alabama  35800 
Attn:  Technical  Director 

Thiokol  Chemical  Corporation  1 

Alpha  Division 
Space  Booster  Plant 
Brunswick,  Georgia  11500 
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Thiokol  Chemical  Corporation  1 

Elkt(in  Division 
Eikton,  Maryland  21921 
Attn.  Librarian 

Thiokol  Chemical  Corporation  1 

Reaction  Motors  Division 
Denville,  New  Jersey  07834 
Attn:  Librarian 

Thiokol  Chemical  Corporation  1 

Rocket  Operations  Center 
P.  O.  Box  lo40 
Ogden,  Utah  84401 
Attn:  Librarian 

Thiokol  Chemical  Corporation  2 

Wasatch  Division 

P.  O.  Box  524 

Brigham  City,  Utah  84302 

Attn:  Library  Section 

Union  Carbide  Chemicals  Co.  1 

P.  O.  Box  8584  (Technical  Center) 
South  Charleston,  West  Virginia 
25303  Attn:  Dr.  H.  W.  Schulz 

United  Aircraft  Corporation  1 

Pratt  (a  Whitney  Fla.  Res.  &  Dev.  Ctr. 
P.  O.  Box  2691 

W.  Palm  Beach,  Florida  33402 
Attn:  Library 

United  Aircraft  Corporation  1 

Corporation  Library 
400  Main  Street 

East  Hartford,  Connecticut  06118 
Attn:  Dr.  David  Rix 

General  Electric  Company  1 

Apollo  Support  Department 
P.O.  Box  2600 

D.'^.yton  Beach,  Florida  32015 


British  Defence  Staff  4 

British  Embassy 

3100  Massachusetts  Avenue 

Washington,  D.  C.  20008 

Attn:  Scientific  Information  Officer 

Defence  Research  Member  4 

Canadian  Joint  Staff  (W) 

2450  Massachusetts  Avenue 
Washington,  D.  C.  20008 

Commanding  Officer  1 

Ammunition  Procurement  &  Supply 
Agency 

Joliet,  Illinois  60400 
Attn:  Engr.  Library 

AFRPL  (RPCS)  1 

Edwards,  California  93523 

Callery  Chemical  Company  1 

Research  &  Development 
Callery,  Pennsylvania  16024 
Attn:  Document  Control 

Ethyl  Corporation  1 

P.  O.  Box  3091 

Istrouma  Branch 

Baton  Rouge,  Louisiana  70805 

Hercules  Powder  Company  1 

Bacchus  Works 
Magna,  Utah  84044 
Attn:  Librarian 

Hynes  Research  Corporation  1 

308  Bon  Air  Avenue 

Durham,  North  Carolina  27704 

Walter  Kidde  Company  1 

675  Main  Street 
Belleville,  New  Jersey  07109 
Attn:  Security  Librarian 
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Olin  Mathieson  Chemical  Corp.  1 
Marion,  Illinois  62959 
Attn:  Research  Library 
Box  508 


Headquarters,  U.  S.  Air  Force  1 
Washington,  D.  C. 

20330 

Attn:  AFRSTD 


Olin  Mathieson  Chemical  Corp.  1 
Research  Library  l-K-3 
275  Winchester  Avenue 
New  Haven,  Connecticut  06511 
Attn:  Mail  Control  Room 

Miss  Laura  M.  Kajuti 

Penns  alt  Chemicals  Corporation  1 
Technological  Center 
900  First  Avenue 

King  of  Prussia,  Pennsylvania  19't06 

Purdue  University  1 

Lafayette,  Indianna  47907 
Attn;  M.  J.  Zucrow 

Rocketdyne,  A  Division  of  1 

North  American  Aviation,  Tnc. 

Solid  Propulsion  Operations 
P.  O.  Box  548,  McGregor,  Texas 
76657 

Attn:  Library 


AFRPL  (RPCL)  1 

Edwards,  California  93523 

Callery  Chemical  Company  1 

Research  &  Development 
Callery,  Pennsylvania  16024 
Attn:  Document  Control 

Chem.  Research  &  Dev.  Center  1 
FMC  Corporation 
P.  O.  Box  8 

Princeton,  New  Jersey  08540 
Attn:  Security  Officer 

Ethyl  Corporation  1 

P.  O.  Box  3091 

Istrouma  Branch 

Baton  Rouge,  Louisiana  70805 

Hynes  Research  Corporation  1 

308  Bon  Air  Avenue 

Durham,  North  Carolina  27704 


Shell  Development  Company 
1400  53rd  Street 
Emeryville,  California  94608 

Southwest  Research  Institute 
Department  of  Structural  Research 
8500  Culebra  Road 
San  Antonio,  Texas  78228 
Attn:  Dr.  Robert  C.  DeHart,  Director 

Texac  ),  Inc.  1 

P.  O.  Box  509 

Beacon,  New  York  12508 

Attn:  Dr.  R.  E.  Conary, Manager 

Los  Alamos  Scientific  Laboratory  1 
University  of  California 
P.  O.  Box  1663 

Los  Alamos,  New  Mexico  87544 


1  Walter  Kidde  Company  1 

675  Main  Street 
Belleville,  New  Jersey  07109 
Attn:  Security  Librarian 

1 

Olin  Mathieson  Chemical  Corp.  1 
Research  Library  l-K-3 
275  Winchester  Avenue 
New  Haven,  Connecticut  06511 
Attn:  Mail  Control  Room 

Miss  Laura  M.  Kajuti 

Pennsalt  Chemicals  Corporation  1 
Technological  Center 
900  First  Avenue 

K  ing  of  Prussia,  Pennsylvania  19406 
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Purdue  University  1 

Lafayette,  Indiana  47907 
Attn:  M.  J.  Zucrow 

Shell  Development  Company  1 

1400  53rd  Street 

Emeryville,  California  94608 

Southwest  Research  Institute  1 

Department  of  Structural  Research 

8500  Culebra  Road 

San  Antonio,  Texas  78228 

Attn:  Dr.  Robert  C.  DeHart,  Director 

Texaco,  Inc.  1 

P.  O.  Box  509 

Beacon,  New  York  12508 

Attn:  Dr.  R.  E.  Conary,  Manager 


Reports  prepared  under  a  Department  of  the  Air  Force  contract  or  project 
will  be  transmitted  VIA: 

Headquarters 

Air  Force  Systems  Command 
Attn:  SCS-41(3-1833) 

Andrews  Air  Force  Base 
Washington,  D.  C.  20331 


